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Abstract 


Reconstruction  of  the  fires  that  occurred  in  the  World  Trade  Center  (WTC)  1 ,  2,  and  7  on 
September  11,  2001,  relied  heavily  on  computer  simulations  because  examination  of  the  post-fire 
premises  was  not  possible  and  the  information  from  eyewitness  accounts  was  severely  limited  in  nature. 
These  simulations  were  perfonned  using  the  National  Institute  of  Standards  and  Technology  (NIST)  Fire 
Dynamics  Simulator  (FDS).  This  report  describes  a  series  of  six  fire  tests  of  single  workstations  (i.e., 
office  cubicles),  the  dominant  combustibles  in  the  WTC  buildings.  The  purpose  was  to  understand  how 
the  office  workstations  may  have  burned  and  to  provide  data  for  improvements  in  the  FDS  combustion 
algorithm  and  its  inputs  needed  to  approximate  the  burning  of  combustibles  as  complex  as  those  that 
comprise  an  office  workstation.  The  variables  in  the  tests  were  the  type  of  workstation,  the  presence  of 
jet  fuel  on  combustible  surfaces,  and  the  presence  of  inert  material  covering  part  of  the  combustible 
surfaces.  The  outcome  was  a  pragmatic,  empirical  approach  to  modeling  the  burning  behavior  that  led  to 
satisfactory  replication  of  the  heat  release  rate  from  the  workstations  and  the  dependence  of  that  rate  on 
the  test  variables. 

Keywords:  Fire,  Fire  Dynamics  Simulator,  fire  modeling,  full-scale  fire  tests,  workstations.  World  Trade 
Center. 
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degree  centigrade 
degree  Faiirenheit  (°F) 
degree  Fahrenheit  (°F) 
kelvin  (K) 


to 

kelvin  (K) 

degree  Celsius  (°C) 

degree  Celsius  (°C) 

kelvin  (K) 

degree  Celsius  (°C) 


Multiply  by 

T/K  =  t/°C  +  273.15 
t/°C  *  t  /deg.  cent, 
t/°C  =  (t/°F-32)/1.8 
T/K  =  (t/°F  + 459.67)/]. 8 
t/°C  =  T/K  2  273.15 


TEMPERATURE  INTERVAL 

degree  Celsius  (°C)  kelvin  (K)  1 .0  E+00 

degree  centigrade  -  degree  Celsius  (°C)  1 .0  E+00 

degree  Fahrenheit  (°F)  degree  Celsius  (°C)  5.555  556  E-01 

degree  Fahrenheit  (°F)  kelvin  (K)  5.555  556  E-01 

degree  Rankine  (°R)  kelvin  (K)  5.555  556  E-01 


VELOCITY  (includes  SPEED) 

foot  per  second  (ft/s)  meter  per  second  (m/s)  3.048  E-01 

inch  per  second  (in. /s)  meter  per  second  (m/s)  2.54  E-02 

kilometer  per  hour  (km/h)  meter  per  second  (m/s)  2.777  778  E-01 

mile  per  hour  (mi/h)  kilometer  per  hour  (km/h)  1 .609  344  E+00 

mile  per  minute  (mi/min)  meter  per  second  (m/s)  2.682  24  E+01 


VOLUME  (includes  CAPACITY) 

cubic  foot  (ft')  cubic  meter  (m') 

cubic  inch  (in.' )  cubic  meter  (m') 

cubic  yard  (yd')  cubic  meter  (m') 

gallon  (U.S.)  (gal)  cubic  meter  (m') 

gallon  (U.S.)  (gal)  liter  (L) 

liter  (L)  cubic  meter  (m') 

ounce  (U.S.  fluid)  (fl  oz)  cubic  meter  (m') 

ounce  (U.S.  fluid)  (fl  oz)  milliliter  (mL) 


2.831  685  E-02 
1.638  706  E-05 
7.645  549  E-01 
3.785  412  E-03 
3.785  412  E+00 
1.0  E-03 
2.957  353  E-05 
2.957  353  E+01 
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Genesis  of  This  Investigation 

Immediately  following  the  terrorist  attack  on  the  World  Trade  Center  (WTC)  on  September  1 1 ,  2001 ,  the 
Federal  Emergency  Management  Agency  (FEMA)  and  the  American  Society  of  Civil  Engineers  began 
planning  a  building  performance  study  of  the  disaster.  The  week  of  October  7,  as  soon  as  the  rescue  and 
search  efforts  ceased,  the  Building  Performance  Study  Team  went  to  the  site  and  began  its  assessment. 
This  was  to  be  a  brief  effort,  as  the  study  team  consisted  of  experts  who  largely  volunteered  their  time 
away  from  their  other  professional  commitments.  The  Building  Performance  Study  Team  issued  its 
report  in  May  2002,  fulfilling  its  goal  "to  determine  probable  failure  mechanisms  and  to  identify  areas  of 
fiiture  investigation  that  could  lead  to  practical  measures  for  improving  the  damage  resistance  of  buildings 
against  such  unforeseen  events." 

On  August  21,  2002,  with  funding  from  the  U.S.  Congress  through  FEMA,  the  National  Institute  of 
Standards  and  Technology  (NIST)  armounced  its  building  and  fire  safety  investigation  of  the  WTC 
disaster.  On  October  1,  2002.  the  National  Construction  Safety  Team  Act  (Public  Law  107-231),  was 
signed  into  law.  The  NIST  WTC  Investigation  was  conducted  under  the  authority  of  the  National 
Construction  Safety  Team  Act. 

The  goals  of  the  investigation  of  the  WTC  disaster  were: 

•  To  investigate  the  building  construction,  the  materials  used,  and  the  technical  conditions  that 
contributed  to  the  outcome  of  the  WTC  disaster. 

•  To  serve  as  the  basis  for: 

-  Improvements  in  the  way  buildings  are  designed,  constructed,  maintained,  and  used; 

-  Improved  tools  and  guidance  for  industry  and  safety  officials; 

-  Recommended  revisions  to  current  codes,  standards,  and  practices;  and 

-  Improved  public  safety. 

The  specific  objectives  were: 

•  Determine  why  and  how  WTC  1  and  WTC  2  collapsed  following  the  initial  impacts  of  the 
aircraft  and  why  and  how  WTC  7  collapsed; 

•  Determine  why  the  injuries  and  fatalities  were  so  high  or  low  depending  on  location, 
including  all  technical  aspects  of  fire  protection,  occupant  behavior,  evacuation,  and 
emergency  response; 

•  Determine  what  procedures  and  practices  were  used  in  the  design,  construction,  operation, 
and  maintenance  of  WTC  1,2,  and  7;  and 

•  Identify,  as  specifically  as  possible,  areas  in  current  building  and  fire  codes,  standards,  and 
practices  that  warrant  revision. 
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NIST  is  a  nonregulatory  agency  of  the  U.S.  Department  of  Commerce's  Technology  Administration.  The 
purpose  of  NIST  investigations  is  to  improve  the  safety  and  structural  integrity  of  buildings  in  the  United 
States,  and  the  focus  is  on  fact  finding.  NIST  investigative  teams  are  authorized  to  assess  building 
performance  and  emergency  response  and  evacuation  procedures  in  the  wake  of  any  building  failure  that 
has  resulted  in  substantial  loss  of  life  or  that  posed  significant  potential  of  substantial  loss  of  life.  NIST 
does  not  have  the  statutory  authority  to  make  findings  of  fault  nor  negligence  by  individuals  or 
organizations.  Further,  no  part  of  any  report  resulting  from  a  NIST  investigation  into  a  building  failure  or 
from  an  investigation  under  the  National  Construction  Safety  Team  Act  may  be  used  in  any  suit  or  action 
for  damages  arising  out  of  any  matter  mentioned  in  such  report  (15  USC  281a,  as  amended  by  Public 
Law  107-231). 

Organization  of  the  Investigation 

The  National  Construction  Safety  Team  for  this  Investigation,  appointed  by  the  then  NIST  Director, 
Dr.  Arden  L.  Bement,  Jr.,  was  led  by  Dr.  S.  Shyam  Sunder.  Dr.  William  L.  Grosshandler  served  as 
Associate  Lead  Investigator,  Mr.  Stephen  A.  Cauffman  served  as  Program  Manager  for  Administration, 
and  Mr.  Harold  E.  Nelson  served  on  the  team  as  a  private  sector  expert.  The  Investigation  included  eight 
interdependent  projects  whose  leaders  comprised  the  remainder  of  the  team.  A  detailed  description  of 
each  of  these  eight  projects  is  available  at  http://wtc.nist.gov.  The  puipose  of  each  project  is  summarized 
in  Table  P-1,  and  the  key  interdependencies  among  the  projects  are  illustrated  in  Fig.  P-1. 


Table  P-1.  Federal  building  and  fire  safety  investigation  of  the  WTC  disaster. 


Technical  Area  and  Project  Leader 

Project  Purpose 

Analysis  of  Building  and  Fire  Codes  and 
Practices;  Project  Leaders:  Dr.  H.  S.  Lew 
and  Mr.  Richard  W.  Bukowski 

Document  and  analyze  the  code  provisions,  procedures,  and 
practices  used  in  the  design,  construction,  operation,  and 
maintenance  of  the  structural,  passive  fire  protection,  and 
emergency  access  and  evacuation  systems  of  WTC  1,  2,  and  7. 

Baseline  Structural  Performance  and 
Aircraft  knpact  Damage  Analysis;  Project 
Leader:  Dr.  Fahim  H.  Sadek 

Analyze  the  baseline  performance  of  WTC  1  and  WTC  2  under 
design,  service,  and  abnormal  loads,  and  aircraft  impact  damage  on 
the  structural,  fire  protection,  and  egress  systems. 

Mechanical  and  Metallurgical  Analysis  of 
Structural  Steel;  Project  Leader:  Dr.  Frank 
W.  Gayle 

Detennine  and  analyze  the  mechanical  and  metallurgical  properties 
and  quality  of  steel,  weldments,  and  connections  from  steel 
recovered  from  WTC  1,  2,  and  7. 

Investigation  of  Active  Fire  Protection 
Systems;  Project  Leader:  Dr.  David 
D.  Evans;  Dr.  William  Grosshandler 

Investigate  the  perfomiance  of  the  active  fire  protection  systems  in 
WTC  1,  2,  and  7  and  their  role  in  fire  control,  emergency  response, 
and  fate  of  occupants  and  responders. 

Reconstruction  of  Thermal  and  Tenability 
Environment;  Project  Leader:  Dr.  Richard 
G.  Gann 

Reconstruct  the  time-evolving  temperature,  thermal  environment, 
and  smoke  movement  in  WTC  1 ,  2,  and  7  for  use  in  evaluating  the 
structural  performance  of  the  buildings  and  behavior  and  fate  of 
occupants  and  responders. 

Structural  Fire  Response  and  Collapse 
Analysis;  Project  Leaders:  Dr.  John 
L.  Gross  and  Dr.  Therese  P.  McAllister 

Analyze  the  response  of  the  WTC  towers  to  fires  with  and  without 
aircraft  damage,  the  response  of  WTC  7  in  fires,  the  performance 
of  composite  steel-trussed  floor  systems,  and  determine  the  most 
probable  structural  collapse  sequence  for  WTC  1,  2,  and  7. 

Occupant  Behavior,  Egress,  and  Emergency 
Communications;  Project  Leader:  Mr.  Jason 
D.  Averill 

Analyze  the  behavior  and  fate  of  occupants  and  responders,  both 
those  who  survived  and  those  who  did  not,  and  the  performance  of 
the  evacuation  system. 

Emergency  Response  Technologies  and 
Guidelines;  Project  Leader:  Mr.  J.  Randall 
Lawson 

Document  the  activities  of  the  emergency  responders  from  the  time 
of  the  terrorist  attacks  on  WTC  1  and  WTC  2  until  the  collapse  of 
WTC  7,  including  practices  followed  and  technologies  used. 
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NIST  WTC  Investigation  Projects 


Figure  P-1.  The  eight  projects  in  the  federal  building  and  fire  safety 
investigation  of  the  WTC  disaster. 


National  Construction  Safety  Team  Advisory  Committee 

The  NIST  Director  also  established  an  advisory  committee  as  mandated  under  the  National  Construction 
Safety  Team  Act.  The  initial  members  of  the  committee  were  appointed  following  a  public  solicitation. 
These  were: 

•  Paul  Fitzgerald,  Executive  Vice  President  (retired)  FM  Global,  National  Construction  Safety 
Team  Advisory  Committee  Chair 

•  John  Barsom,  President,  Barsom  Consulting,  Ltd. 

•  John  Bryan,  Professor  Emeritus,  University  of  Maryland 

•  David  Collins,  President,  The  Preview  Group,  Inc. 

•  Glenn  Corbett,  Professor,  John  Jay  College  of  Criminal  Justice 

•  Philip  DiNenno,  President,  Hughes  Associates,  Inc. 
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•  Robert  Hanson,  Professor  Emeritus,  University  of  Michigan 

•  Charles  Thornton,  Co-Chairman  and  Managing  Principal,  The  Thomton-Tomasetti  Group, 
Inc. 

•  Kathleen  Tiemey,  Director,  Natural  Hazards  Research  and  Applications  Information  Center, 
University  of  Colorado  at  Boulder 

•  Forman  Williams,  Director,  Center  for  Energy  Research,  University  of  California  at  San 
Diego 

This  National  Construction  Safety  Team  Advisory  Committee  provided  technical  advice  during  the 
Investigation  and  commentary  on  drafts  of  the  Investigation  reports  prior  to  their  public  release.  NIST 
has  benefited  from  the  work  of  many  people  in  the  preparation  of  these  reports,  including  the  National 
Construction  Safety  Team  Advisory  Committee.  The  content  of  the  reports  and  recommendations, 
however,  are  solely  the  responsibility  of  NIST. 

Public  Outreach 

During  the  course  of  this  Investigation,  NIST  held  public  briefings  and  meetings  (listed  in  Table  P-2)  to 
solicit  input  from  the  public,  present  preliminary  findings,  and  obtain  comments  on  the  direction  and 
progress  of  the  Investigation  from  the  public  and  the  Advisory  Committee. 

NIST  maintained  a  publicly  accessible  Web  site  during  this  Investigation  at  http://wtc.nist.gov.  The  site 
contained  extensive  information  on  the  background  and  progress  of  the  Investigation. 

NIST's  WTC  Public-Private  Response  Plan 

The  collapse  of  the  WTC  buildings  has  led  to  broad  reexamination  of  how  tall  buildings  are  designed, 
constructed,  maintained,  and  used,  especially  with  regard  to  major  events  such  as  fires,  natural  disasters, 
and  terrorist  attacks.  Reflecting  the  enhanced  interest  in  effecting  necessary  change,  NIST,  with  support 
from  Congress  and  the  Administration,  has  put  in  place  a  program,  the  goal  of  which  is  to  develop  and 
implement  the  standards,  technology,  and  practices  needed  for  cost-effective  improvements  to  the  safety 
and  security  of  buildings  and  building  occupants,  including  evacuation,  emergency  response  procedures, 
and  threat  mitigation. 

The  strategy  to  meet  this  goal  is  a  three-part  NIST-led  public-private  response  program  that  includes: 

•  A  federal  building  and  fire  safety  investigation  to  study  the  most  probable  factors  that 
contributed  to  post-aircraft  impact  collapse  of  the  WTC  towers  and  the  47-story  WTC  7 
building,  and  the  associated  evacuation  and  emergency  response  experience. 

•  A  research  and  development  (R&D)  program  to  (a)  facilitate  the  implementation  of 
recommendations  resulting  from  the  WTC  Investigation,  and  (b)  provide  the  technical  basis 
for  cost-effective  improvements  to  national  building  and  fire  codes,  standards,  and  practices 
that  enhance  the  safety  of  buildings,  their  occupants,  and  emergency  responders. 
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Table  P-2.  Public  meetings  and  briefings  of  the  WTC  Investigation. 


Date 

Location 

Principal  Aj^enda 

June  24. 2002 

New  York  City,  NY 

Public  meeting:  Public  comments  on  the  Druji  Plan  for  the 
pending  WTC  Investigation. 

August  21.  2002 

Gaithersburg,  MD 

Media  briefing  announcing  the  formal  start  of  the  Investigation. 

December  9,  2002 

Washington,  DC 

Media  briefing  on  release  of  the  Public  Update  and  NIST  request 
for  photographs  and  videos. 

April  8,  2003 

New  York  City,  NY 

Joint  public  forum  with  Columbia  University  on  first-person 
interviews. 

April  29-30.  2003 

Gaithersburg,  MD 

NCST  Advisory  Committee  meeting  on  plan  for  and  progress  on 
WTC  Investigation  with  a  public  comment  session. 

May  7,  2003 

New  York  City,  NY 

Media  briefing  on  release  of  May  2003  Progress  Report. 

Aueust  26-'' 7  ''003 

Gaithersburg,  MD 

NCST  Advmorv  Committpp  mpptinp  on  status  of  thp  WTC 
investigation  with  a  public  comment  session. 

September  17,2003 

New  York  City,  NY 

Media  and  public  briefing  on  initiation  of  first-person  data 
collection  nroiects 

V-VllW  11  \w/ll    L/l  \J  1  V  ^  lO  . 

December  2-3  ''003 

Gaithersburg,  MD 

NICEST  AdvTsorv  (^OTTimittpp  mpptinp  on  ^status  anH  initial  rpsiills 

-1.  ^  > —       A    1  \ vj  V  ijyjy  y    \_  w iiiiiiiii^^  iiiwwiiii ti.  v.'Ji  oicituiij  uii VJ  iiiiiiui  iwoi.J.ii^ 

and  release  of  the  Public  Update  with  a  public  comment  session. 

February  12,  2004 

New  York  City,  NY 

Public  meeting  on  progress  and  preliminary  findings  with  public 
cnmmpnts  on  i^ssiips  to  bp  ron'sidprpH  in  finrmiilatinp  finril 

W  V-fll  111        11  li3    Wll    1  k  J  i. )  LI      *  7    IV-/    KJ^    \-/\J  lltj             I              111    1  \.yl  1 1 1 1.1 1  CI  11 11  X111L4.1 

recommendations. 

June  18  ''004 

■JUlllVp'     lLJi|    J—\J  V7  T 

New  York  Citv  NY 

Media/public  briefing  on  release  of  June  2004  Progress  Report. 

June  22-23,  2004 

Gaithersburg,  MD 

NCST  Advisory  Committee  meeting  on  the  status  of  and 
nrplimmarv  findinfr*;  from  thp  WTC  Invpstipation  with  a  niiblir 

L/l  ^  lliillilCil  y    XlllUlllc^iJ  XI  will  lllW    VV    1         Ul  V  Voll^ClLlV.^11    Wllll  Cl  L./LIL/IIW 

comment  session. 

August  24,  2004 

Northbrook,  IL 

Public  viewing  of  standard  fire  resistance  test  of  WTC  floor 
svstpm  at  1  Jndprwritprs  I  abnratnrip*s  Inr 

ljV  01^111    Clt           llVJ^l  VVlllwlkJ    i  >  Li      V_'  1  CI  I V/1  1Wl>^  X11^> 

October  19-20,  2004 

Gaithersburg,  MD 

NCST  Advisory  Committee  meeting  on  status  and  near  complete 
set  of  preliminary  findings  with  a  public  comment  session. 

November  22,  z(jU4 

Gaithersburg,  MD 

NCST  Advisory  Committee  discussion  on  draft  annual  report  to 
Congress,  a  public  comment  session,  and  a  closed  session  to 
Hisniss  nrp-Hraft  rprommPTiHatinns  for  WTC  Investipation 

VJl  3^.1.1.31?    L.*!^    VJl  Ul  I   1      V.  VI  i  11 1 1^  1  lUU  11  V-f  liO   1  Vl      T  T    1  > —    111  v  ^Oli              V-Zll  • 

April  5,  2005 

New  York  City,  NY 

Media  and  public  briefing  on  release  of  the  probable  collapse 
sequence  for  the  WTC  towers  and  draft  reports  for  the  projects  on 
codes  and  practices,  evacuation,  and  emergency  response. 

June  23, 2005 

New  York  City,  NY 

Media  and  public  briefing  on  release  of  all  draft  reports  for  the 
WTC  towers  and  draft  recommendations  for  public  comment. 

September  12-13, 
2005 

Gaithersburg,  MD 

NCST  Advisory  Committee  meeting  on  disposition  of  public 
comments  and  update  to  draft  reports  for  the  WTC  towers. 

September  13-15, 
2005 

Gaithersburg,  MD 

WTC  Technical  Conference  for  stakeholders  and  technical 
community  for  dissemination  of  findings  and  recommendations 
and  opportunity  for  public  to  make  technical  comments. 

•    A  dissemination  and  technical  assistance  program  (DTAP)  to  (a)  engage  leaders  of  the 
construction  and  building  community  in  ensuring  timely  adoption  and  widespread  use  of 
proposed  changes  to  practices,  standards,  and  codes  resulting  from  the  WTC  Investigation 
and  the  R&D  program,  and  (b)  provide  practical  guidance  and  tools  to  better  prepare  faciUty 
owners,  contractors,  architects,  engineers,  emergency  responders,  and  regulatory  authorities 
to  respond  to  future  disasters. 

The  desired  outcomes  are  to  make  buildings,  occupants,  and  first  responders  safer  in  future  disaster 
events. 
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National  Construction  Safety  Team  Reports  on  the  WTC  Investigation 

A  final  report  on  the  collapse  of  the  WTC  towers  is  being  issued  as  NIST  NCSTAR  1.  A  companion 
report  on  the  collapse  of  WTC  7  is  being  issued  as  NIST  NCSTAR  1  A.  The  present  report  is  one  of  a  set 
that  provides  more  detailed  documentation  of  the  Investigation  findings  and  the  means  by  which  these 
technical  results  were  achieved.  As  such,  it  is  part  of  the  archival  record  of  this  Investigation.  The  titles 
of  the  full  set  of  Investigation  pubhcations  are: 

NIST  (National  Institute  of  Standards  and  Technology).  2005.  Federal  Building  and  Fire  Safety 
Investigation  of  the  World  Trade  Center  Disaster:  Final  Report  on  the  Collapse  of  the  World  Trade 
Center  Towers.  NIST  NCSTAR  1.  Gaithersburg,  MD,  September. 

NIST  (National  Institute  of  Standards  and  Technology).  2006.  Federal  Building  and  Fire  Safety 
Investigation  of  the  World  Trade  Center  Disaster:  Final  Report  on  the  Collapse  of  World  Trade  Center  7. 
NIST  NCSTAR  1  A.  Gaithersburg,  MD. 

Lew,  H.  S.,  R.  W.  Bukowski,  and  N.  J.  Carino.  2005.  Federal  Building  and  Fire  Safety  Investigation  of 
the  World  Trade  Center  Disaster:  Design,  Construction,  and  Maintenance  of  Structural  and  Life  Safety 
Systems.  NIST  NCSTAR  1-1.  National  Institute  of  Standards  and  Technology.  Gaithersburg,  MD, 
September. 

Fanella,  D.  A.,  A.  T.  Derecho,  and  S.  K.  Ghosh.  2005.  Federal  Building  and  Fire  Safety 
Investigation  of  the  World  Trade  Center  Disaster:  Design  and  Construction  of  Structural  Systems. 
NIST  NCSTAR  1-1  A.  National  Institute  of  Standards  and  Teclinology.  Gaithersburg,  MD, 
September. 

Ghosh,  S.  K.,  and  X.  Liang.  2005.  Federal  Building  and  Fire  Safety  Investigation  of  the  World 
Trade  Center  Disaster:  Comparison  of  Building  Code  Structural  Requirements.  NIST 
NCSTAR  1-lB.  National  Institute  of  Standards  and  Technology.  Gaithersburg,  MD,  September. 

Fanella,  D.  A.,  A.  T.  Derecho,  and  S.  K.  Ghosh.  2005.  Federal  Building  and  Fire  Safety 
Investigation  of  the  World  Trade  Center  Disaster:  Maintenance  and  Modifications  to  Structural 
Systems.  NIST  NCSTAR  1-lC.  National  Institute  of  Standards  and  Technology.  Gaithersburg, 
MD,  September. 

Grill,  R.  A.,  and  D.  A.  Johnson.  2005.  Federal  Building  and  Fire  Safety  Investigation  of  the  World 
Trade  Center  Disaster:  Fire  Protection  and  Life  Safety  Provisions  Applied  to  the  Design  and 
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E.1  INTRODUCTION 


Reconstruction  of  the  fires  that  occurred  in  the  World  Trade  Center  (WTC)  buildings  on 
September  11,  2001,  relied  heavily  on  computer  simulations  since  examination  of  the  post-fire  premises 
was  not  possible  and  the  infomiation  from  eyewitness  accounts  was  severely  limited  in  nature.  These 
simulations  were  performed  using  the  National  Institute  of  Standards  and  Technology  (NIST)  Fire 
Dynamics  Simulator  (FDS),  version  4.  FDS,  a  computational  fluid  dynamics  (CFD)  model,  solves 
numerically  a  form  of  the  Navier-Stokes  equations  with  subgrid-scale  modeling  appropriate  for  the  low- 
speed,  thermally  driven  flows  prevalent  in  fires.  The  prevalent  combustion  model  in  FDS  assumed  that 
the  generation  of  fuel  vapors  was  determined  by  a  heat  of  gasification  and  that  the  combustion  rate  was 
determined  by  the  rate  of  mixing  of  these  vapors  and  air. 

To  prepare  FDS  for  use  in  reconstructing  the  WTC  fires,  it  was  necessary  to  establish  its  accuracy  in 
simulating  the  burning  of  the  principal  combustibles  in  the  WTC  buildings.  The  most  common  floor 
layout  was  an  open  space  populated  by  a  large  array  of  workstations  or  cubicles.  These  approximately 
2.4  m  (8  ft)  square  units  were  geometrically  complex,  containing  a  desktop,  file  storage,  bookshelves, 
carpeting,  computer,  chair,  and  a  variety  of  amounts  and  locations  of  paper.  Each  cubicle  was  bounded 
on  all  four  sides  by  privacy  panels. 


E.2         DESCRIPTION  OF  FIRE  TESTS 


This  report  describes  a  series  of  six  fire  tests  of  individual  cubicles.  The  purpose  was  to  understand  how 
the  office  workstations  may  have  burned  and  to  provide  data  for  improvements  in  the  FDS  combustion 
algorithm  and  its  inputs  needed  to  approximate  the  burning  of  combustibles  as  complex  as  those  that 
comprise  an  office  workstation.  The  tests  involved: 

•  Instrumentation  to  measure  the  heat  release  rate  (HRR)  and  the  mass  loss  rate.  The  entire 
cubicle  was  placed  on  a  load  cell  under  the  NIST  9  m  by  12  m,  10  MW  collector  hood.  The 
hood  was  instrumented  to  measure  HRR  by  oxygen  consumption  calorimetry. 

•  A  skirted  ceiling  placed  over  the  workstation,  enabling  the  fire  effluent  to  form  a  hot  ceihng 
layer  such  as  would  have  occurred  in  a  room  fire. 

•  Placement  of  four  video  cameras  to  enable  following  the  progress  of  the  fires.  Strategically 
placed  heat  flux  gauges  and  thermocouples  provided  indicators  of  the  intensity  of  the  fire  and 
its  growth  pattern. 

•  Ignition  by  a  2  MW  spray  burner,  directed  downward  over  a  2  m  square  pan  simulating  an 
adjacent  burning  workstation. 
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•  Two  different  workstations  designs:  a  generic  workstation  analogous  to  those  used  in  the 
WTC  and  a  workstation  identical  to  one  used  on  floors  93  to  100  of  WTC  1.  The  visible 
differences  between  the  two  were  in  the  nature  of  the  privacy  panels,  the  file  cabinets,  and  the 
chair.  The  former  had  a  combustible  mass  of  250  kg;  the  latter  contained  about  15  kg  less 
wood  laminate  desk  surface  but  was  otherwise  similar. 

•  Inclusion  of  the  effects  of  inert  debris  from  fallen  ceiling  tiles  or  fractured  walls  and  jet  fuel 
(Jet  A)  similar  to  that  from  the  incident  aircraft  in  the  disaster. 

For  input  to  the  FDS  modeling,  the  thermophysical  properties  of  the  cubicle  materials  were  determined 
using  a  cone  calorimeter.  Values  of  HRR  and  ignition  delay  time  were  obtained  for  piloted  ignition  at 
heat  fluxes  of  20  kW/m",  25  kW/m"  (for  the  desk  surface  only  ),  45  kW/m",  and  70  kW/m".  None  of  the 
materials  yielded  a  constant  HRR.  However,  for  modeling  purposes,  the  samples  were  treated  as  if  they 
were  exhibiting  a  square-wave  HRR.  The  data  were  fit  to  a  linear  variation  of  HRR  with  incident  flux. 
The  ignition  data,  supplemented  by  the  heat  flux  values  at  which  the  ignition  delay  time  exceeded  1200  s 
(taken  as  non-ignition),  were  used  to  infer  an  ignition  temperature  and  effective  thermal  inertia  for  each 
material. 

E.3         FIRE  TEST  RESULTS 


An  initial  test  with  a  half  workstation  gauged  the  behavior  of  what  was  an  unknown  system.  The  results 
confirmed  that  the  chosen  ignition  source  would  provide  an  orderly  and  rapid  progression  of  fire  growth 
and  that  the  peak  fire  size  would  stay  within  the  capacity  of  the  calorimeter  hood. 

The  properties  determined  in  each  test  included  the  time  of  ignition  of  the  first  item;  the  total  mass  lost 
and  the  mass  loss  rate  over  the  duration  of  the  test;  the  total  heat  released  by  the  burning  workstation,  the 
HRR  over  the  duration  of  the  test,  and  the  net  HRR.  (the  total  HRR  minus  the  HRR  of  the  burner);  the 
burning  time;  and  the  effective  heat  of  combustion.  A  detailed  description  of  the  progress  of  each  fire 
was  also  prepared.  The  numerical  results  and  their  uncertainties  are  shown  in  Table  ES-1 .  These  data 
were  usable  in  modeling  the  open  burning  of  the  types  of  workstations  prevalent  in  the  WTC  buildings. 

The  following  observations  characterize  the  phenomena  that  controlled  the  spread  of  flame  on 
workstation  components  and  the  spread  from  component  to  component: 

•    In  Tests  1  through  Test  4  (no  Jet  A),  the  workstation  combustion  began  with  paper  on  the 
desktop  autoigniting  from  thermal  radiation  from  the  burner  flames.  The  flames  on  the  paper 
stacks  then  piloted  the  ignition  of  the  computer  monitor  shell,  the  desktop,  and  the  top  of  the 
chair  back.  The  chair  ignited  the  carpet  by  dropping  flaming  material  on  it  and  ignited  the 
underside  of  the  desk  surface  by  direct  flame  impingement.  The  desk  surfaces  tended  to  be 
slow  to  ignite  but  provided  most  of  the  heat  release  in  the  latter  parts  of  the  tests.  The  peak 
HRR  corresponded  to  all  accessible  combustible  surfaces  burning  simultaneously. 
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Executive  Summary 


•  The  jet  fuel  considerably  shortened  the  time  to  the  peak  HRR  value.  In  Test  5,  the  radiant 
flux  from  the  flames  vaporized  the  Jet  A  until  it  was  ignited  by  a  random  piece  of  flaming 
material  from  the  privacy  panel  adjacent  to  the  burner.  In  Test  6,  the  initial  event  was 
autoignition  of  the  paper  on  the  desktop,  followed  by  ignition  of  the  Jet  A.  Both  fires 
proceeded  rapidly  toward  involvement  of  all  the  combustibles.  Here  again,  the  peak  HRR 
corresponded  to  all  accessible  combustible  surfaces  burning  simultaneously. 

•  In  all  tests,  the  post-peak  behavior  included  both  the  burning  out  of  the  various  combustibles 
and  the  collapsing  of  three-dimensional  components  into  more  compact  piles  that  decreased 
their  burning  rates.  ,  - 

A  sensitivity  analysis  of  the  test  results  provided  the  following: 

•  Approximately  75  percent  of  the  heat  release  and  mass  loss  occurred  over  a  period  of 
approximately  20  inin,  insensitive  to  all  the  test  variables. 

•  The  peak  HRR  of  the  workstation  was  diminished  in  rough  proportion  to  the  fraction  of  total 
combustible  surface  area  covered  by  inert  materials.  It  was  estimated  that  a  moderate 
increase  in  the  surface  coverage  would  decrease  the  peak  HRR  proportionately  and  that  a 
very  large  increase  in  coverage  would  have  a  more  than  proportionate  effect  but  the  latter 
would  have  occurred  only  in  areas  that  were  highly  rubblized. 

•  The  presence  of  jet  fuel  had  only  a  small  effect  on  the  total  mass  loss  (approximately  200  kg), 
the  total  heat  released  (approximately  4  GJ)  and  the  peak  HRR  (approximately  7  MW)  from 
the  burning  workstation.  The  peak  HRR  value  was  estimated  to  be  insensitive  to  a  lesser 
amount  of  jet  fuel  present,  but  sensitive  to  a  larger  amount.  [The  Jet  A  in  the  tests  was 
consumed  shortly  after  speeding  the  ignition  of  the  cubicle.  More  Jet  A  would  have  resulted 
in  a  large  pool  of  a  highly  volatile  combustible.] 

•  The  presence  of  4  L  of  Jet  A  dramatically  decreased  the  times  to  peak  HRR,  from 
approximately  10  min  to  approximately  3  min.  Increasing  the  quantity  of  jet  fuel  was  not 
expected  to  affect  further  the  time  to  full  involvement  of  the  cubicle. 

•  The  tiles  and  Jet  A  were  nearly  independent  factors:  the  effect  of  the  tiles  on  the  burning  rate 
of  the  cubicle  was  almost  independent  of  the  presence  or  absence  of  the  Jet  A. 

•  The  overall  burning  behavior  was  not  sensitive  to  workstation  design,  with  the  only  (small) 
effect  being  due  to  the  difference  in  mass  and  the  exposed  flammable  area. 

•  Variation  in  the  mass  of  paper  on  the  cubicle  surfaces  would  not  have  affected  the  ignition  of 
the  heavier  workstation  components.  However,  the  total  absence  of  exposed  paper  would 
have  delayed  the  workstation  ignition  until  the  next  most  susceptible  component  autoignited 
or  random  flaming  debris  from  adjacent  fire  areas  served  as  a  pilot.  Additional  paper  in  the 
file  cabinets  would  have  prolonged  the  burning  time  of  the  workstation,  but  would  not  likely 
have  affected  the  peak  HRR  due  to  the  poor  air  supply  within  the  cabinets. 
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The  five  measures  of  fire  intensity  (HRR,  mass  loss  rate,  carbon  dioxide  volume  fraction  in  the  exhaust 
gases,  temperature  of  the  exhaust  gases  and  radiant  flux  seen  by  the  gage  viewing  the  entire  fire  from  the 
south  side  of  the  cubicle)  presented  a  consistent  representation  of  the  progress  of  the  fire,  given  their 
individual  measurement  characteristics. 


E.4         SIMULATION  OF  FIRE  TESTS 


Preliminary  simulations  indicated  that  the  traditional  representation  of  materials  combustion  was  not 
sufficient  to  replicate  the  experimental  results  with  sufficient  accuracy.  This  was  not  surprising  since 
during  the  fire  tests  some  important  phenomena  had  been  identified  as  not  being  capable  of  inclusion  in  a 
computational  model  of  current  sophistication.  These  included:  changes  in  the  shape  (e.g.,  by  collapse) 
and  surface  area  of  combustibles  that  affected  their  burning  rate;  three-dimensional  burning  behavior, 
especially  of  the  chair;  fabric  peeling  on  the  privacy  panels;  and  anticipation  of  whether  a  combustible 
would  ignite  from  direct  flame  impingement  or  by  radiation-driven  autoignition. 

A  pragmatic,  empirical  approach  was  taken  to  the  representation  of  the  various  combustibles. 

•  The  carpet,  desk,  and  privacy  panel  data  from  the  cone  calorimeter  were  used  directly.  The 
carpet  and  privacy  panel  were  modeled  assuming  a  burning  rate  proportional  to  the  heat  flux 
from  the  surrounding  environment.  The  desk  was  modeled  as  a  charring  solid,  in  which  a 
pyrolysis  front  propagates  through  the  material  leaving  a  layer  of  char  behind  that  insulates 
the  material  and  reduces  the  burning  rate. 

•  All  the  other  items  within  the  workstation  were  modeled  as  a  single  entity  by  lumping  their 
masses  together  into  three  homogeneous  "boxes"  and  distributing  these  throughout  the 
workstation.  The  combustion  properties  were  adjusted  to  match  the  results  of  the 
experiments.  Thus,  the  results  from  the  single  workstation  bums  served  to  calibrate  the 
model.  The  properties  of  each  box  were: 

-  Dimensions  (length,  width,  height):  0.5  m  by  0.5  m  by  1.0  m  (20  in.  by  20  in.  by  39  in.) 

-  Ignition  temperature:  370  °C 

-  HRR:  2.25  MW  after  a  100  s,  linear  ramp  from  the  point  of  ignition.  The  heat  release 
then  decreased  in  time  resulting  from  a  decrease  in  surface  area  as  the  boxes  burned. 

Simulations  of  Tests  2,  3,  5,  and  6  were  then  performed  to  determine  how  well  the  above  approximations 
led  to  agreement  between  the  measured  and  computed  HRR  profiles.  The  magnitudes  and  shapes  of  the 
HRR  versus  time  curves  were  similar  in  form  and  magnitude,  with  the  following  notes: 

•  For  Tests  2  and  3  (no  Jet  A),  the  simulations  underestimated  the  times  to  the  HRR  peak  by 
about  2  min.  This  was  due  to  the  strong  influence  of  the  melting  of  the  chair  plastic  onto  the 
carpet,  a  feature  not  captured  in  the  model. 

•  The  predicted  decrease  in  the  HRR  peak  resulting  from  tiles  and  the  increase  in  the  peak  from 
the  Jet  A  were  similar  in  magnitude  to  the  measured  results. 
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•    The  simulations  correctly  predicted  the  small  increase  in  the  time  to  the  peak  with  the 

addition  of  the  tiles,  but  they  did  not  predict  the  more  than  a  factor  of  three  reduction  in  the 
time  to  the  peak  for  the  addition  of  Jet  A.  FDS  consumed  the  Jet  A  immediately  with 
relatively  small  effect  on  the  growth  of  the  overall  HRJl.  For  the  actual  fire,  the  ignition  of 
the  Jet  A  occurred  a  bit  later,  but  then  quickly  caused  the  entire  cubicle  to  become  fully 
involved.  The  error  in  the  time  to  the  HRR  peak  was  less  than  2  min. 

These  effects  were  small  enough  to  be  averaged  out  when  multiple  workstations  were  burning 
simultaneously  at  various  stages  of  their  burning.  It  was  estimated  that  on  the  order  of  40  workstations 
were  burning  on  a  given  floor  at  the  height  of  the  fires  in  WTC  1 .  Thus,  the  chosen  set  of  parameters  and 
approximated  component  burning  descriptions  gave  a  reasonable  description  of  the  actual  workstation 
heat  release  rate  behavior  and  its  dependence  on  the  factors  varied  in  the  tests  (inert  tile  coverage  and 
presence  of  Jet  A). 
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Introduction 


1.1  CONTEXT  OF  TESTS 


Reconstruction  of  the  fires  that  occurred  in  World  Trade  Center  (WTC)  1,  2,  and  7  on 
September  11,  2001,  relies  heavily  on  computer  simulations  since  examination  of  the  post-fire  premises 
was  not  possible  due  to  the  collapse  of  the  buildings,  and  the  information  from  eyewitness  accounts  is 
severely  limited  in  nature.  These  simulations  were  to  be  performed  using  the  National  Institute  of 
Standards  and  Technology  (NIST)  Fire  Dynamics  Simulator  (FDS),  version  4  (NIST  NCSTAR  1-5F'). 
FDS,  a  computational  fluid  dynamics  (CFD)  model,  solves  numerically  a  form  of  the  Navier-Stokes 
equations  with  a  subgrid-scale  modeling  appropriate  for  the  low-speed,  thermally  driven  flows  prevalent 
in  fires. 

The  traditional  combustion  model  in  FDS  is  relatively  simple.  Enthalpy  input  (from  an  ignition  source  or 
an  existing  fire)  pyrolyzes  the  exposed  combustibles  or  fuels.  The  rate  of  gasification  and  the  ignition 
condition  for  a  combustible  are  determined  by  its  thennophysical  properties  and  ignition  temperature. 
The  air  space  in  the  compartment  is  divided  into  rectilinear  volume  elements.  Each  of  these  elements  is 
characterized  by  its  mixture  fraction,  a  conserved  scalar  quantity  defined  as  the  mass  fraction  of  the 
contained  gas  that  originated  as  fuel.  If  an  element  is  bounded  by  elements  with  mixture  fractions  that  are 
rich  (fijel-air  equivalence  ratio  greater  than  unity)  and  lean  (fuel-air  equivalence  ratio  less  than  unity), 
then  combustion  is  deemed  to  occur  within  that  element.  FDS  assumes  that  the  combustion  rate  is 
mixing-controlled,  and  that  the  reaction  of  fuel  and  oxygen  is  infinitely  fast.  The  heat  released  by  the 
combustion  is  dispersed  according  to  the  calculated  fluid  dynamics,  the  radiant  flux  to  the  surfaces  of 
combustibles  are  recalculated,  and  the  process  is  repeated  until  all  the  combustibles  have  been  consumed. 
The  "flames"  spread  as  successive  volume  elements  reach  combustible  mixture  fractions.  The  mass 
fractions  of  the  major  combustion  products  are  derived  from  the  mixture  fraction  by  means  of  "state 
relations,"  empirical  expressions  arrived  at  by  a  combination  of  simplified  analysis  and  measurement. 

To  prepare  FDS  for  use  in  reconstructing  the  WTC  fires,  it  was  necessary  to  establish  its  accuracy  in 
simulating  the  combustion  of  the  principal  combustibles  in  the  three  buildings.  Conversations  with 
building  occupants  and  managers  indicated  some  variety  of  furnishings  and  arrangements  in  the  various 
tenant  areas  (NIST  NCSTAR  1-5).  The  most  common  floor  layout  was  a  continuous  open  space 
populated  by  a  large  array  of  workstations  or  cubicles.  These  units  were  geometrically  complex  and 
contained  a  mixture  of  combustible  materials.  Each  approximately  2.44  m  (8  ft)  square  cubicle  typically 
was  bounded  on  all  four  sides  by  privacy  panels,  with  an  entrance  opening  in  one  side  only.  Within  the 
area  defined  by  the  panels  was  a  self-contained  workspace:  desktop,  file  storage,  bookshelves,  carpeting, 
chair,  etc.  Presumably,  there  were  a  variety  of  amounts  and  locations  of  paper,  both  exposed  on  the  work 
surfaces  and  contained  within  the  file  cabinets  and  bookshelves.  These  cubicles  were  grouped  in  clusters 


'  This  reference  is  to  one  of  the  companion  documents  from  this  Investigation.  A  list  of  these  documents  appears  in  the  Preface 
to  this  report. 
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or  rows,  with  up  to  210  units  on  a  given  floor.  On  a  mass  basis,  the  workstations  were  generally  the 
dominant  combustibles  on  the  floors  where  extensive  fires  were  observed. 

This  report  describes  a  series  of  six  fire  tests  of  individual  cubicles  performed  during  July  and 
August  2003.  The  most  important  characteristic  of  the  combustion  is  the  heat  release  rate  (HRR)  of  the 
fire.  It  dominates  the  growth  of  the  fire  and  the  concomitant  hazards.  The  mass  loss  rate  (MLR),  also 
measured,  provides  parallel  information,  since  for  a  constant  heat  of  combustion,  the  MLR  and  HRR  are 
proportional.  Madrzykowski  (1996,  1998)  had  measured  the  heat  release  rate  behavior  of  workstations 
that  were  somewhat  different  from  those  studied  here.  He  found  the  peak  heat  release  rate  values  varied 
with  such  details  as  the  number  of  privacy  panels  and  the  failure  of  a  thermoplastic  shelf  support  leading 
to  a  shelf  collapse  which  spewed  loose  papers  onto  the  fire,  immediately  enhancing  it  substantially.  Here, 
two  different  cubicle  designs  were  combusted.  As  a  result  of  the  airplane  impacts  into  the  two  towers, 
workstations  in  some  locales  would  have  been  showered  with  inert  debris  from  the  ceiling  tile  system  or 
fractured  walls  and  also  may  have  been  sprayed  with  Jet  A  (jet  fuel).  Thus,  these  two  effects  were  also 
included  in  the  test  series.  The  tests  are  documented  in  Chapter  2.  The  results  are  presented  and 
discussed  in  Chapter  3. 

The  thermophysical  properties  of  the  cubicle  materials  (carpet,  desk  [wood],  computer  monitor,  chair, 
privacy  panel,  and  stacked  paper)  were  determined  using  the  cone  calorimeter  (ASTM  E  1354)  for  input 
to  the  FDS  simulations.  The  results  of  these  measurements  are  included  in  Appendix  A. 

Preliminary  simulations  indicated  that  the  traditional  representation  of  materials  combustion  was  not 
sufficient  to  replicate  the  experimental  results  with  sufficient  accuracy.  An  empirically  derived  revision 
was  incorporated  and  is  described  in  Chapter  4,  along  with  the  resulting  successful  simulations. 
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Description  of  Experiments 


2.1  WORKSTATIONS 
2.1.1        Generic  Workstations 

The  number  of  different  types  of  workstations  in  the  World  Trade  Center  (WTC)  buildings  is  unknown 
and  probably  large.  Marsh  &  McLennan,  the  tenant  in  floors  93  to  100  of  WTC  1,  had  at  least  four 
different  workstation  layouts  (Fleck  2003).  While  there  were  likely  to  have  been  some  differences  in  fire 
behavior  among  these,  their  general  complexity  and  variety  of  combustible  materials  make  any  one  of 
them  a  reasonable  subject  for  examination  of  the  Fire  Dynamics  Simulator  (FDS)  predictive  capability. 

Five  of  the  six  experiments  focused  on  the  burning  behavior  of  a  type  of  single  workstation  that  was  a 
generic  analog  (not  a  replica)  of  one  type  of  station  used  on  floors  93  to  100  of  WTC  1.  Figures  2-1  and 
2-2  show  the  generic  workstation  layout,  as  tested. 


Source:  NISI. 

Figure  2-1.  View  of  the  north  and  east  sides  of  the  generic  work  station. 
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Source:  NIST. 

Figure  2-2.  View  from  above  of  the  west  side  of  the  generic  workstation.  The 
horizontal  insulated  pipe  above  the  west  wall  is  the  feed  line  to  the  four 

nozzles  of  the  spray  igniter. 

Figure  2-1  is  a  view  through  the  cubicle  opening  at  the  right  hand  (east)  side  of  the  layout,  showing  the 
bookshelf  on  the  rear  (north)  wall,  and  the  two-drawer  file  cabinets  beneath  the  east  desk  surface. 
Figure  2-2  is  an  elevated  view  looking  down  at  the  left  hand  (west)  side  of  the  cubicle.  The  computer 
processor  and  the  wastebasket,  both  beneath  the  desk  surface,  are  visible  on  the  right.  The  insulated  feed 
line  containing  the  four  fuel  spray  nozzles  is  visible  above  the  top  of  the  west  cubicle  wall;  the  fire  these 
produce  is  actually  on  the  outer  side  of  the  west  cubicle  wall.  Note  the  positioning  of  the  chair, 
immediately  adjacent  to  the  west  desk  surface  and  with  a  portion  of  its  seat  below  the  desk  surface  section 
on  which  the  keyboard  rests.  Note  also  that  paper,  in  various  sheet  counts,  is  distributed  around  the  three 
closed  sides  of  the  cubicle  with  the  largest  concentration  being  in  the  bookcase. 

The  layout,  including  the  placement  of  the  various  non-stationary  items,  was  suggested  by  personnel  from 
the  company  that  supplied  office  furnishings  to  the  occupants  of  the  noted  floors  of  WTC  1  (Fleck  2003). 
The  mass  and  distribution  of  papers  and  other  office  clutter  is  based  on  information  also  provided  by  this 
source,  who  was  a  frequent  visitor  to  these  offices. 
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The  workstation  was  a  nominally  2.44  m  by  2.44  m  (8  ft  by  8  ft)  cubicle  with  privacy  wall  panels  1 .22  m 
(4  ft)  high,  except  for  one  1.52  m  (5  ft)  tall  panel  section  supporting  the  bookcase.  The  entry  opening  was 
1.22  m  (4  ft)  wide. 

•  The  panels  were  covered  on  both  sides  with  a  thermoplastic  fabric.  Their  interior  structure 
was  composed  of  layers  of  fiberglass  and  perforated  steel.  Each  panel  section  had  a  38  mm 
(1.50  in.)  softwood  frame  around  its  periphery  that  was  covered  in  part  by  fabric  and  in  part 
by  a  steel  outer  frame. 

•  The  work  surfaces  were  formed  from  four  unequal  sections  of  laminated  medium  density 
fiberboard  supported  by  brackets  from  the  wall  panels.  The  tops  of  the  surfaces  were  0.75  m 
(29.5  in.)  above  the  floor.  In  the  assemblies  used  for  all  but  one  of  the  tests,  the  brackets 
supplied  were  composed  of  an  aluminum  alloy  which  could  be  expected  to  melt  relatively 
early  in  the  fire.  Since  the  actual  WTC  stations  had  steel  support  brackets,  steel  brackets 
were  added  here  as  well  to  prevent  premature  collapse  of  the  work  surface.  The  tops  of  the 
panel  sections  were  wired  together  for  the  same  reason;  aluminum  clips,  not  present  in  the 
actual  WTC  work  stations,  were  the  normal  link  holding  the  panels  together  in  the  generic 
work  stations. 

•  The  "bulldog"  office  chair  was  upholstered  with  a  non-thermoplastic  fabric  over  polyurethane 
foam  (seat  and  seat  back)  supported  by  a  one-piece  thermoplastic  shell.  Its  five-legged  base 
was  also  thermoplastic  with  steel  framing  and  support  elements. 

•  The  bookcase  was  1.22  m  (48  in.)  long  and  had  a  steel  shelf  and  top.  These  were  supported 
on  their  ends  by  combustible  end  panels.  The  steel  front  closure  panel  was  fabric-covered 
steel,  and  it  was  open  (on  top  of  the  bookcase). 

•  The  carpet  tiles  were  nylon  fiber-faced  over  a  dense  foam  rubber  backing.  A  square  area 
2.74  m  by  2.74  m  (9  ft  by  9  ft)  was  covered  with  36  carpet  tiles. 

•  The  computer  monitor  was  a  nominally  17  in.  CRT-based  unit.  Its  front  face  was  taped  with 
fiberglass  tape,  and  it  was  pointed  toward  the  wall  panel  opposite  the  cubicle  opening  for 
safety  in  the  event  of  an  implosion.  The  keyboard  was  placed  in  its  normal  location,  parallel 
to  the  angled  segment  of  the  work  surface.  The  computer  processor  (tower-type  container 
with  plastic  only  on  the  front  face  of  the  container)  was  placed  on  the  floor  next  to  a  waste 
paper  basket  (both  on  the  north  side  of  the  cubicle,  opposite  the  opening). 

•  The  wastebasket  was  thermoplastic  and  it  contained  one  ream  of  copier  paper  atop  five 
balled-up  paper  ream  wrappers. 

•  Copier  paper  was  chosen  to  substitute  for  all  types  of  paper-based  products  expected  in  a 
work  station  (loose  leaf,  catalogs,  books,  etc.)  since  it  is  reproducible  and  similar  in  general 
fire  behavior  to  all  such  products.  The  paper  was  20  lb,  84  brightness.  Besides  being  in  the 
wastebasket,  copier  paper  was  distributed  elsewhere  as  indicated  in  Figs.  2-1  and  2-2.  It  was 
used  both  in  horizontal  and  in  vertical  arrays,  arrangements  that  potentially  might  have 
behaved  quite  differently,  but  did  not  appear  to  do  so  in  any  of  the  tests  here.  The  paper  in 
the  bookcase  and  the  paper  in  the  southwest  comer  (to  the  left  of  the  cubicle  opening  as 
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indicated  in  Fig.  2-2)  were  placed  vertically  in  open,  thin  cardboard  "document  boxes,"  each 
with  an  89  mm  (3.5  in.)  wide  interior  dimension.  These  boxes  had  a  sloped  top  edge  such 
that  they  equaled  the  long  dimension  of  the  paper  in  the  rear  but  were  only  100  mm  (4  in.)  tall 
in  front.  In  the  bookcase,  each  document  box  contained  1 .33  reams  of  copier  paper  (a  rather 
loose  degree  of  packing).  The  10  boxes  were  grouped  in  sets  of  two  with  equal, 
approximately  3  cm  ( 1 .2  in.)  spacing  between  each  pair.  The  four  boxes  on  the  desk  surface 
contained  one  ream  each  of  copier  paper  (an  even  looser  degree  of  packing)  and  were 
grouped  as  a  set  of  four  units  in  a  line,  in  contact  with  each  other.  Other  paper  was  laid  as 
open  horizontal  stacks  on  the  desk  surface  at  locations  seen  in  Figs.  2-1  and  2-2.  A  few 
sheets  were  tacked  (as  groups  of  three  to  five  sheets)  to  the  cubicle  walls  on  three  sides. 

•    The  three  lateral  file  cabinets  were  0.91  m  wide,  0.51  m  deep,  and  0.68  m  high  (36  in.  by 
20  in.  by  27  in.).  Each  had  two  painted  steel  drawers  and  rested  directly  on  the  carpet  tiles. 
Each  cabinet  was  capable  of  holding  up  to  about  150  kg  (330  lb)  of  paper.  [Such  paper 
content  is  normally  included  in  the  fire  load  estimates  of  an  office  occupancy  but,  when 
enclosed  in  steel,  it  is,  in  effect,  partially  discounted  (AISI  1981;  NFPA  2003).]  Here,  two  of 
the  cabinets  contained  a  very  limited  amount  of  copier  paper  as  a  rough  means  of  assessing 
the  extent  to  which  paper  in  file  drawers  might  contribute  to  a  fire.  It  was  anticipated  that  this 
contribution  would  be  small  and  occur  well  past  the  fire  peak,  thus  having  no  significant 
influence  on  important  portions  of  the  HRR  curves. 

Two  piles  were  laid  horizontally  in  two  file  drawers.  Each  pile  consisted  of  one  ream  laid  in  the  drawer 
comer  with  the  long  dimension  of  the  paper  aligned  with  the  depth  dimension  of  the  drawer  topped  by  a 
second  ream  aligned  at  90  degrees  to  the  lower  ream.  Both  placements  were  in  the  position  most 
vulnerable  to  heat  exposure,  i.e.,  they  were  against  the  north  wall  of  the  file  cabinets  and  thus 
immediately  adjacent  to  an  area  subject  to  intense  burning,  the  open  area  under  the  desk  surface. 

Table  2-1  lists  the  mass  of  each  type  of  workstation  component  and  an  estimate  of  the  combustible  mass 
that  each  represented.  In  Table  2-2,  the  materials  are  grouped  into  four  categories.  Wood-based 
materials  accounted  for  nearly  half  of  the  combustible  mass;  paper,  one  quarter;  and  synthetic  materials, 
the  remainder.  The  heats  of  combustion  of  these  materials,  based  on  the  cone  calorimeter  results  reported 
in  Appendix  A,  varied  over  a  range  of  nearly  a  factor  of  three. 
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Table  2-1.  Contents  of  generic  workstation. 


Component 

Mass  (kg) 

Combustible 
Fraction 

Combustible 
Mass  (kj^) 

Work  surface  -  melamine  laminate  over  medium  density 
fiber-board  4  pieces,  6.1  m  total  length  x  0.61  m 
wide  X  28  mm  thick 

82.8 

1.0 

82.8 

27  reams  of  paper  and  14  document  boxes 

63.7 

1.0 

63.7 

Plastic  kick  plates  and  trim  (base  of  walls,  inside  and 
outside) 

7.1 

1.0 

7.1 

Computer  keyboard 

1.2 

1.0 

1.2 

Plastic  waste  basket 

0.7 

1.0 

0.7 

36  nylon  carpet  tiles  with  rubber  backing 

38.0 

0.9 

34.2 

Shelf  ends  -  particle  board  or  dense  foam 

3.8 

0.9 

3.4' 

Bulldog  chair  -  fabric,  foam,  thermoplastic  shell  and  base 

19.4 

0.8 

15.5 

Computer  monitor  -  ABS 

17.6 

0.3 

5.3 

Cnmniitpr  nrocessor 

12.3 

0.3 

3.7 

9  wall  panels  with  aluminum  angle,  wood  frames, 
fiberglass,  and  metal  mesh,  3  sizes 

168.2 

0.18 

30.3 
S  ^  W  wooH 
5.0  kg  fabric 

Book  shelf 

8.3 

0.1 

0.8 

3  two-drawer  steel  filing  cabinets,  0.91  m  long  x  0.51  m 
deep  X  0.76  m  high 

142.5 

0.0 

0.0 

Total 

557.1 

0.45 

248.7 

a.  Table  2-2  assumes  that  the  3.4  kg  is  plastic. 


Table  2-2.  Categories  of  mat 

erials  in  the  generic  workstations. 

Material 

Mass  (kg) 

Fraction  of  Total  Mass 

Effective  Heat  of 
Combustion  (MJ/kg) 

Wood/laminate 

1 11.5 

0.45 

14 

Paper 

63.7 

0.25 

14 

Plastics" 

39.3 

0.16 

16  to  38 

Carpet 

34.2 

0.14 

22 

a.  Includes  computer  monitor  shell  (16  MJ/kg),  wall  fabric  (30  MJ/kg),  and  chair  composite  (38  MJ/kg). 


2.1.2        WTC  Workstation 

One  of  the  six  tests  was  of  a  workstation  that  was  identical  in  design  to  one  of  those  on  floors  93  to  100  of 
WTC  1 .  The  principal  differences  between  it  and  the  generic  workstations  used  in  the  other  five  tests 
were  as  follows: 

•    The  wall  panel  had  a  3  mm  (0. 125  in.)  layer  of  flame-retarded  polyester  fiber  beneath  the 
outer  fabric,  a  more  open  steel  panel  beneath  this,  a  central  fiberboard  layer  (3  mm  thick)  and 
an  all  steel  (i.e.,  no  wood)  peripheral  frame.  The  fiberboard  roughly  cut  to  one  half  the 
amount  of  woody  fuel  within  the  wall  panels  of  the  enclosure.  This  decrease  in  woody  fuel 
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was  only  about  5  percent  of  the  total  available  in  the  desk  surfaces.  Also,  its  enclosure  deep 
within  the  wall  panels  delayed  its  burning. 

•  The  four  file  cabinets  had  a  total  face  length  of  2.67  m  ( 1 05  in. ),  rather  than  three  with  a  total 
face  length  of  2.73  m  (107  in.).  They  had  a  flammable,  charring  plastic  surface  on  the  drawer 
fronts  that  added  about  10  percent  to  the  exposed  flammable  area. 

•  The  chair  was  constructed  somewhat  differently.  Its  seat  and  back  were  separate  pieces,  it 
had  an  additional  3  kg  (6.6  lb)  of  mass,  and  it  behaved  as  if  its  upholstered  surfaces  were 
flame-retarded.  The  chair  was  certified  to  pass  California  Technical  Bulletin  117,  which 
usually  requires  somewhat  flame-retarded  polyurethane  foam. 

Table  2-3  summarizes  the  mass  of  the  components  and  their  estimated  combustible  mass  content  for  the 
WTC  workstation;  Table  2-4  groups  the  materials  into  the  same  categories  as  in  Table  2-2.  Items  in  bold 
type  differed  from  the  generic  workstation  summarized  in  Table  2-1.  Note  that  the  WTC  workstation  had 
a  slightly  lower  total  combustible  mass,  but  the  plastic  fronts  on  the  file  cabinets  meant  there  was  more 
exposed  area  of  flammable  material  and  the  total  mass  of  plastic  was  greater. 


Table  2-3.  Contents  of  WTC  workstation. 


Component 

Mass  (kg) 

Combustible 
Fraction 

Combustible 
Mass  (kg) 

Work  surface — melamine  laminate  over  medium  density 
fiberboard,4  pieces,  6.1  m  total  length  x  0.61  m  wide  x  30  mm 
thick 

79.4 

1.0 

79.4 

27  reams  of  paper  and  14  document  boxes 

63.7 

1.0 

63.7 

Computer  keyboard 

1.2 

1.0 

1.2 

Plastic  waste  basket 

0.7 

1.0 

0.7 

36  nylon  carpet  tiles  with  rubber  backing 

38.0 

0.9 

34.2 

Plastic  doors  for  2  cm  to  76  cm  wide  shelves 

4.4  (est.) 

1.0 

4.4 

High-end  bulldog  chair — fabric,  polyurethane  foam, 
thermoplastic  shell  and  base 

22.4 

0.8 

17.9 

Computer  terminal,  ABS 

17.6 

0.3 

5.3 

Computer  processor 

12.3 

0.3 

3.7 

6  w  all  panels  of  steel  construction  with  an  0.32  cm  thick 
hardboard  septum,  fiberglass  insulation,  and  fabric 

89.2 

0.23 

20.5 
15.4  kg  wood 
5.1  kg  fabric 

Three  2-door  filing  cabinets  similar  to  those  above  plus  a 
smaller  portable  filing  cabinet 

136.9 

0.02  (est.) 

2.7  (est.) 

Total 

465.8 

0.50 

233.7 

Note:  Bold  items  have  values  different  from  those  in  the  generic  workstation. 
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Table  2-4.  Categories  of  materials  in  the  WTC  workstation. 


Material 

Mass  (kg) 

Fraction  of  Total 

h  IfV^rtiv'o  ilt'iit  f if  r '<kinhiist ifiri 
(M.l/k«) 

Wood/laminate 

94.8 

0.40 

14 

Paper 

63.7 

0.27 

14 

Plastics" 

41.0 

0.18 

16  to  38 

Carpet 

34.2 

0.15 

22 

a.  Includes  computer  monitor  shell  (16  MJ/kg),  wall  fabric  (30  MJ/kg),  and  chair  composite  (38  MJ/kg). 
Note:  Bold  items  have  values  different  from  those  in  the  generic  workstation. 

2.2  TEST  ASSEMBLY 

Each  cubicle  rested  on  top  of  a  double  layer  of  13  mm  (VS  in.)  thick  calcium  silicate  sheets.  The  support 
structure  beneath  this  was  in  turn  placed  atop  a  set  of  four  weighing  cells,  one  at  each  comer. 


Since  a  hot,  smoky  upper  gas  layer  in  a  room  can  play  a  significant  role  in  fire  spread,  a  3.66  m  by  3.66  m 
(12  ft  by  12  ft)  ceiling  was  created  here.  It  consisted  of  a  single  layer  of  13  mm  (1.2  in.)  thick  calcium 
silicate  board  arrayed  2.74  m  (9  ft)  above  the  floor  of  the  workstation.  The  calcium  silicate  board  was 
supported  by  a  frame  formed  from  steel  pipe.  Water  flowing  through  the  pipe  precluded  heat  damage  to 
the  frame.  There  were  only  four  water  pipes  flowing  beneath  the  ceiling  to  minimize  heat  loss  from  the 
hot  gases  trapped  beneath  the  ceiling.  The  calcium  silicate  layer  was  surrounded  on  all  four  sides  by  a 
0.61  m  (2  ft)  high  steel  skirt.  This  allowed  the  formation  of  a  hot  smoke  layer  above  the  workstation, 
providing  for  radiative  heating  by  the  accumulated  hot  gases. 

The  entire  assembly  was  placed  beneath  the  hood  of  the  NIST  10  MW  calorimeter  hood  (9  m  by  12  m)  so 
as  to  capture  the  full  smoke  plume  and  allow  heat  release  rate  measurements  during  the  burning  of  the 
workstation.  The  air  supply  to  the  fire  was  thus  not  restricted,  and  the  burning  behavior  was  that  of  a 
fully  open  condition. 

2.3  INSTRUMENTATION 

The  types  and  location  of  the  various  instruments  used  in  this  test  series  are  listed  in  Table  2-5.  The  key 
measurement  for  this  study  is  the  heat  release  rate  (HRR)  of  the  fire  throughout  a  test.  The  HRR 
measurement  instrumentation  and  calculation  process  is  essentially  the  same  as  for  the  NIST  3  MW 
calorimeter  (Bryant  et  al.  2003)  and  is  briefly  summarized  here. 
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Table  2-5.  Test  instrumentation. 

Type  Instrument 

Location 

Heat  release  rate 

8  bi-directional  probes,  each  with  associated  thermocouple,  for  mass  flow 
detennination  m  exhaust  duct;  associated  pressure  transducers  in  enclosure  on 
roof;  crossed  Tee  gas  sampling  probe  in  exhaust  duct; 
oxygen,  COt  and  CO  analyzers  in  calorimeter  control  room 

4  load  cells 

Positioned  at  corners  of  cubicle  support  structure 

Water-cooled  Schmidt-Boelter 
total  heat  flux  gauge 

Near  the  center  of  the  cubicle  floor  pointing  upward 

Garden  total  heat  flux  gauge 

2  m  height  and  about  3  m  from  the  cubicle  side  with  the  entrance,  to  assist 
with  the  timing  of  certain  events  in  the  fire  growth. 

7  themiocouples  with  high 
temperature  insulation  for  lead 
wires 

Positioned  on  the  underside  of  the  desk  to  monitor  flame  spread  in  some  tests 

4  video  cameras 

2  positioned  to  view  from  the  east,  one  from  the  south,  and  one  from  the 
southwest. 

As  shown  by  Huggett  (1980),  for  viitually  all  materials  found  in  buildings,  the  average  heat  of 
combustion  per  gram  of  oxygen  consumed  is  13.1  kcal.  The  determination  of  oxygen  consumption 
begins  with  the  collection  of  all  the  fire  effluent  in  the  hood  and  includes  several  continuous 
measurements  in  the  hood  exhaust:  total  mass  flow  and  the  concentrations  of  oxygen,  CO2  and  CO.  The 
latter  two  provide  a  measure  of  the  completeness  of  combustion,  which  affects  the  conversion  of  oxygen 
consumption  rate  to  HRR. 

The  total  exhaust  mass  flow  was  derived  from  measurement  of  the  exhaust  flow  velocity  and  temperature, 
combined  with  the  composition  of  the  exhaust  stream.  The  flow  velocity  was  measured  at  eight  locations 
in  the  2  m  (6.5  ft)  diameter  exhaust  duct  using  bi-directional  probes  arrayed  in  the  vertical  direction.  The 
temperature  was  also  measured  at  each  of  these  locations  using  1.6  mm  diameter  sheathed 
chromel/alumel  thermocouples. 

The  exhaust  gas  was  sampled  through  a  four-armed  cross  that  spanned  the  duct  diameter  vertically  and 
horizontally,  with  each  arm  having  four  6  mm  (%  in.)  diameter  holes  spaced  at  equal  duct  area  points 
along  the  arm  span.  The  sampling  lines  and  diaphragm  pump  were  all  heated  above  350  K  to  avoid 
condensation  in  the  sampling  line.  A  dry  ice  trap  followed  by  a  Drierite  dessicant  trap  removed  moisture 
before  the  flow  entered  a  paramagnetic  type  oxygen  monitor  and  a  nondispersive  infrared  analyzer  for  CO 
andCO.. 

The  total  mass  of  the  cubicle  was  monitored  continuously  during  the  tests.  The  decline  in  total  mass 
provided  a  measure  of  the  fuel  gasification  rate  (or  mass  loss  rate,  MLR)  and  enabled  computation  of  the 
effective  heat  of  combustion  of  the  cubicles  and  their  contents.  The  cubicle  (including  its  contents  and 
flooring)  rested  on  a  platform  that  was  supported  on  a  network  of  four  load  cells.  Each  load  cell  had  a 
total  weight  capacity  of  900  kg  (2000  lb)  and  was  thermally  protected  from  the  fire  by  calcium  silicate 
board  enclosures.  The  signals  generated  by  the  individual  load  cells  were  measured  by  a  load  cell 
terminal  for  intermediate  display  and  retransmitted  to  the  data  acquisition  system.  During  test 
preparation,  the  output  from  each  load  cell  terminal  was  monitored  to  insure  the  load  cells  were  not 
overloaded. 
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In  Tests  1,  2,  and  3,  a  Schmidt-Boelter  total  heat  flux  gage  was  placed  on  the  cubicle  floor  on  its  cast-west 
centerline  and  0.8  m  from  the  cubicle  opening.  The  sensor  face  was  127  mm  (5  in.)  above  the  carpet  and 
it  measured  the  local  heat  flux  to  the  carpet.  The  gage  was  cooled  by  hot  (68      to  78  ''C)  water  to 
preclude  condensation  of  water  vapor  on  the  gage  surface.  An  external  Gardon-type  total  heat  flux  gage 
was  also  mounted  externally  approximately  1.2  m  (4  ft)  from  the  cubicle  and  approximately  1  m  (3  ft) 
above  the  cubicle  floor  level  so  as  to  view  the  entire  fire  from  the  cubicle  opening  side.  This  provided  a 
signal  that  was  roughly  proportional  to  the  instantaneous  HRR  of  the  fire,  useful  for  checking  the  timing 
of  fire  events  and  assessing  whether  the  calorimeter  properly  estimated  the  peak  HRR. 

In  Tests  2,  3,  and  4,  thermocouples  were  placed  on  the  underside  of  the  desk  surface  in  an  attempt  to 
follow  the  progress  of  the  fire  in  that  important  location.  However,  despite  the  use  of  high  temperature 
glass  fiber  insulation,  the  thennocouple  leads  short  circuited,  and  no  usable  data  were  obtained. 

The  signals  from  all  the  measurement  devices  were  collected  by  a  high  speed  data  acquisition  system. 
The  signal  from  each  device  was  recorded  every  1  s.  During  this  recording  cycle,  each  data  point  was 
averaged  over  a  period  of  1/200  s. 

Each  test  was  also  monitored  by  four  video  cameras  located  as  shown  in  Fig.  2-3.  These  enabled 
observation  of  the  progression  of  flame  spread  over  the  various  surfaces  in  the  workstation  assembly. 
This  was  useful  in  explaining  the  nature  of  the  measured  HRR  curve  over  time  and  any  differences  from 
the  predicted  HRR  behavior.  A  complete  documentation  of  the  fire  growth  was  not  possible  due  to 
blockage  of  the  view  by  the  cubicle  partitions  and  interference  by  the  flames  themselves  as  the  fire  grew 
in  intensity.  An  observer  narrated  the  fire  growth  to  one  of  the  cameras  to  supplement  what  the  cameras 
recorded  directly. 
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Figure  2-3.  Location  of  video  cameras. 
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2.4         TEST  PARAMETERS 

After  extensive  discussion,  two  factors  were  identified  that  might  affect  the  fire  behavior  in  the  tests  in  a 
manner  that  would  be  a  good  test  of  the  FDS  simulation  capability  and  that  at  the  same  time  were 
germane  to  fire  behavior  in  the  WTC  buildings: 

•  Presence  and  Quantity  of  Aviation  Fuel  (Jet  A).  The  impact  of  the  airplanes  could  be 
expected  to  leave  varying  amounts  of  Jet  A  deposited  on  various  surfaces,  with  some  areas 
being  heavily  engulfed  and  others,  further  from  the  impact  zone,  being  dry.  Because  of  its 
high  volatility  and  heat  of  combustion,  even  a  limited  amount  of  Jet  A  has  the  potential  to 
alter  the  burning  and  spread  rates  of  the  fire.  The  quantity  chosen,  4  L  per  workstation,  was 
sufficient  to  fonn  a  continuous  layer  on  the  major  horizontal  surfaces  and  thus  was  sufficient 
to  alter  the  fire  growth  process  substantially. 

•  Presence  and  Quantity  of  Inert  Rubble  on  Surfaces.  Survivors  of  the  WTC  disaster  reported 
significant  falling  of  ceiling  tiles,  even  on  floors  some  distance  from  the  airplane  impacts. 
These  tiles  are  both  low  in  flammability  (up  to  1 5  percent  by  mass  combustible  content)  and 
good  thermal  insulators.  As  such,  they  are  capable  of  retarding  or  even  preventing  the  fire 
involvement  of  surfaces  onto  which  they  fall.  In  two  of  the  tests,  the  tiles  were  simulated  by 
40  13  mm  (1/2  in.)  thick  slabs  of  calcium  silicate  board  cut  into  15  cm  by  30  cm  (6  in.  by 

12  in.)  rectangles.  These  were  distributed  on  the  desk,  chair  and  floor  surfaces  as  shown  in 
Fig.  2-4.  The  tiles  were  placed  only  on  areas  open  to  the  ceiling  (including  part  of  the  chair 
seat  and  part  of  the  tops  of  some  of  the  paper  piles),  but  not  under  the  desk  or  chair.  To 
obtain  a  measurable  but  not  overwhelming  effect,  the  tiles  covered  approximately  30  percent 
of  this  upward-facing  area. 


Source:  NISI. 

Figure  2-4.  View  from  the  south  of  the  generic  work  station 
with  calcium  silicate  "tiles"  in  place  on  exposed,  upward 
facing  surfaces  (floor,  desk,  chair  seat). 
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2.5  IGNITION  SCENARIO 


As  with  any  complex  assemblage  of  fuels,  the  heat  release  rate  history  for  a  workstation  subsequent  to 
ignition  can  be  expected  to  depend  significantly  on  the  mode  and  location  of  that  ignition  event.  In  the 
WTC  towers,  the  photographic  evidence  shows  an  overwhelming  ignition  event  in  the  area  of  impact. 
This  fire  receded,  then  returned  and  spread  to  other,  less  disrupted  areas,  providing  a  continuing  scries  of 
large  ignition  sources.  This,  combined  with  the  rationale  that  only  the  larger  fires  could  heat  enough 
structural  members  sufficiently  to  influence  the  eventual  building  collapse,  suggested  that  a  high  intensity 
ignition  source  for  these  tests  would  be  a  proper  choice. 

The  ignition  source  chosen  was  a  2  MW,  four-nozzle  spray  burner.  The  nozzles  were  directed  downward 
over  a  1  m  by  2  m  pan,  creating  a  fire  that  simulated  the  burning  of  an  adjacent  workstation.  As  shown  in 
Fig.  2-2,  the  burner  and  pan  were  located  immediately  adjacent  to  the  exterior  of  the  west  panel  of  the 
test  station,  with  the  pan  bottom  being  0.81  m  (2  ft  8  in.)  above  the  floor.  The  fuel  spray  was  somewhat 
coarse,  and  some  droplets  made  it  to  the  pan,  although  there  was  no  significant  liquid  accumulation  there. 
The  size  and  placement  of  this  ignition  source  were  guided  by  preliminary  FDS  predictions.  While  a  real 
workstation  would  yield  a  transient  HRR  and  thus  a  time-varying  heat  input  to  its  neighbor,  this 
complexity  was  deemed  undesirable  for  the  task  at  hand. 

The  igniter  fire  was  supported  by  a  flow  of  commercial-grade  liquid  heptanes  (mixed  isomers,  but 
referred  to  as  "heptane"  for  the  remainder  of  this  report)  sprayed  at  a  nominally  steady  rate  from  the  four 
nozzles.  The  heptane  was  supplied  from  a  reservoir  by  a  variable-speed  pump.  The  desired  flow  was 
preset  before  the  test  and  measured  in  triplicate  by  catching  the  flow  from  each  nozzle  in  a  volumetric 
cylinder  for  a  period,  typically,  of  20  s.  The  nominally  2  MW  fire  plume  supported  by  the  heptane  flow 
typically  impinged  almost  continually  on  the  ceiling  above  the  igniter.  There  was  thus  a  near-continuous 
wall  of  flames  radiating  toward  the  workstation  along  the  central  three-fourths  of  the  length  of  one  wall 
panel.  In  addition,  the  workstation  was  subject  to  radiation  from  the  hot  ceiling  and  the  hot  smoke 
captured  below  the  ceiling.  Typically  (with  one  exception)  the  spray  burner  was  kept  on  until  late  in 
the  fire. 


2.6         EXPERIMENTAL  DESIGN 


Since  the  test  series  was  to  examine  the  effects  of  two  factors,  with  two  values  of  each,  the  basic  test 
series  utilized  a  full  factorial  design.  These  four  tests  are  depicted  in  Table  2-6  as  Tests  2,  3,  5,  and  6. 
This  is  the  minimum  number  of  tests  required  to  provide  a  first  order  statistical  assessment  of  the  impact 
of  two  parameters.  These  tests  were  preceded  by  a  test  that  utilized  nominally  one  half  of  a  generic 
workstation,  though  it  included  both  a  full  chair  and  full  computer.  This  was  done  to  gage  the  behavior  of 
what  was  an  unknown  system.  A  sixth  test  compared  the  results  of  Test  2  with  the  behavior  of  a  different 
style  of  workstation. 
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Table  2-6.  Experimental  test  plan. 


Test 

Specimen 

Tiles 

Jet  Fuel 

1 

Half  generic  workstation 

None 

None 

2 

Generic  workstation 

None 

None 

3 

Generic  workstation 

40 

None 

4 

WTC  workstation 

None 

None 

5 

Generic  workstation 

None 

4L 

6 

Generic  workstation 

40 

4L 

The  results  of  Test  1,  using  the  half  workstation,  confirmed  that  the  chosen  ignition  source  would  provide 
an  orderly  and  rapid  progression  of  fire  growth  over  the  various  fuel  surfaces  in  the  workstation.  It  also 
demonstrated  that  the  peak  fire  size,  while  large,  would  most  likely  stay  within  the  capacity  of  the 
calorimeter  hood. 

2.7         TEST  PROCEDURE 

The  following  was  the  sequence  followed  in  performing  each  of  the  tests: 

•  The  workstation  was  typically  assembled  shortly  before  a  test.  The  humidity  in  the  Large 
Fire  Laboratory  at  the  time  of  the  testing  was  approximately  70  percent  for  most  of  the  tests. 
Cellulosic  materials  are  hygroscopic,  and  significantly  increased  moisture  content  could  slow 
the  ignition  process  compared  to  a  more  nornial  relative  humidity  of  50  percent.  Thus,  the 
most  susceptible  such  material,  the  copier  paper,  was  covered  with  plastic  sheeting  if  it  were 
to  be  exposed  for  more  than  2  h.  This  probably  did  not  preclude  significant  moisture  pick-up 
in  the  outer  portions  of  paper  piles. 

•  A  22.7  kg  (50  lb)  weight  was  added  to  the  platform  (and  then  removed)  to  verify  the  proper 
differential  output  of  the  load  cell  electronics. 

•  The  data  channels  to  be  scanned  were  verified,  the  appropriate  delay  times  for  the  hood 
measurements  were  selected,  span  values  (mainly  calibration  gas  cylinder  volume  fractions) 
were  entered,  as  were  the  calorimetry  settings  for  the  combustion  products  expansion  factor 
and  the  generic  heat  of  combustion  per  mass  of  oxygen. 

•  The  gas  analysis  instruments  in  the  calorimeter  were  zeroed  and  spanned  before  each  test, 
using  known  gas  mixtures  with  oxygen,  nitrogen,  COt  and  CO.  Procedures  generally 
followed  those  described  by  Bryant  et  al.  (2003). 

•  The  background  heat  release  rate  was  measured  for  a  few  minutes  immediately  before  the  test 
was  to  start  in  order  to  establish  the  calorimeter  basehne.  Similarly,  the  baseline  values  of  the 
thermocouples  and  heat  flux  gages  were  checked. 

•  The  heptane  flow  was  measured  in  triplicate  during  the  zero/spanning  process. 

•  In  Tests  4  and  6,  the  panels  representing  the  fallen  tiles  were  laid  in  place. 
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•  In  Tests  5  and  6,  the  jet  fuel  was  placed  on  the  horizontal  surfaces  immediately  prior  to  the 
start  of  the  test.  A  sprinkling  can  was  used  in  two  separate  operations,  each  involving  2  L  of 
the  liquid  fuel.  First  the  liquid  was  sprinkled  on  the  horizontal  desk  surfaces  and  the  objects 
on  them  (i.e.,  the  various  paper  stacks  or  document  boxes,  the  computer  monitor,  and  the 
keyboard)  using  a  timed  movement  that  attempted  to  allot  an  equal  amount  of  liquid  to  each 
one-third  of  the  total  surface.  The  desk  surface  had  been  leveled  after  installation  so  that  the 
liquid  would  not  run  preferentially  in  one  direction.  Next,  an  equal  amount  of  the  liquid  was 
sprinkled,  in  a  similar  timed  manner,  on  the  central  open  section  of  the  carpet  (not  under  the 
desk  surfaces).  Since  the  chair  occupied  a  portion  of  this  space,  the  allotment  for  that  portion 
went  onto  the  chair  seat  and  back  surfaces.  For  the  objects  on  or  in  contact  with  the  desk 
surface,  there  was  some  tendency  for  the  Jet  A  to  wick  into  them  if  they  were  porous.  This 
was  true  of  the  paper,  the  inert  tiles,  and  the  wall  panel  fabric  just  above  the  desk  surface. 

•  With  all  lights  and  cameras  on,  the  test  was  initiated  by  starting  the  heptane  flow  and 
immediately  igniting  it  with  a  torch.  This  defined  time  zero. 

•  The  heptane  flow  to  the  igniter  was  left  on  until  well  past  the  peak  in  heat  release  unless  there 
was  an  indication  that  the  fire  was  going  to  exceed  the  10  MW  calorimeter  capacity.  This 
happened  only  in  Test  5.  In  that  case,  the  flow  was  turned  off  approximately  575  s  into  the 
test. 

•  The  background  HRR  was  measured  for  a  few  minutes  immediately  after  the  test  in  order  to 
determine  the  degree,  if  any,  of  calorimeter  baseline  shift. 

•  The  residual  weights  of  the  paper  piles  in  the  two  file  drawers  and  of  the  paper  piles  on  the 
desk  were  obtained  as  a  measure  of  the  participation  of  this  paper  in  the  overall  fire. 
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3.1  DATA  REDUCTION 


The  standard  data  analysis  program  used  in  these  tests  computed  the  heat  release  rate  (HRR)  as  a  function 
of  time  during  the  tests.  These  results  allowed  the  preparation  of  the  HRR  plots  (Figs.  3-1  through  3-6) 
using  a  standard  plotting  package.  Each  graph  shows  both  the  full  heat  release  rate  curves  of  a  fire  test 
and  an  annotated  portion  that  omits  most  of  the  tail-off  of  the  fire.  Figure  3-7  is  a  sequence  of 
photographs  taken  during  Test  2. 

Additional  analysis  was  required  for  computing  the  net  heat  release  rate,  the  total  heat  released  by  the 
workstation,  the  mass  burning  rate,  and  the  heat  of  combustion.  The  methods  for  doing  these 
supplemental  calculations  are  given  in  succeeding  sections.  A  comprehensive  discussion  of  the  derivation 
of  the  heat  release  equation,  the  basis  for  the  selection  of  the  constant  parameters,  and  the  data  analysis 
program,  including  the  real  time  output  corrected  for  the  flow  time  delay  and  the  instrument  response 
time,  can  be  found  in  Bryant  et  al.  (2003). 

The  results  of  both  these  calculations  and  observations  of  the  video  tapes  are  included  in  Table  3-1. 
3.1.1        Burning  Time  and  Ignition  Time 

Two  measures  were  used  for  the  burning  time.  The  full-width  at  the  half-height  of  the  net  heat  release 
rate  (FWHH)  was  used  as  a  measure  of  the  intense  period  of  burning.  The  time  at  which  the  mass  loss  of 
the  workstation  including  Jet  A  decreased  by  75  percent  the  total  mass  loss,  t(75  %),  was  also  used  as  a 
measure  of  the  time  interval  for  intense  burning.  These  two  quantities  were  obtained  from  the  recorded 
mass  loss  data  and  firom  the  reduced  heat  release  rate  data.  The  net  heat  release  data  had  pronounced 
oscillations  past  the  peak  and  there  were  two  or  more  times  corresponding  to  the  half-height  of  the  net 
heat  release  weight.  The  average  of  the  shortest  and  longest  times  was  used  in  computing  the  FWHH. 

The  ignition  time  of  the  first  item,  t,o,  was  determined  from  analyzing  the  video  from  the  tests. 
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Time  from  Ignition  (s) 

Figure  3-1  a.  Complete  heat  release  rate  curve,  Test  1,  one  half  of  generic  workstation, 

no  Jet  A,  no  tiles. 
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Figure  3-1  b.  Annotated  partial  heat  release  rate  curve,  Test  1,  one  half  of  generic 

workstation,  no  Jet  A,  no  tiles. 
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Figure  3-2a.  Complete  HRR  curve,  Test  2,  complete  generic  workstation,  no  Jet  A, 

no  tiles. 
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Figure  3-2b.  Annotated  partial  HRR  curve.  Test  2,  complete  generic  workstation, 

no  Jet  A,  no  tiles. 
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Figure  3-3a.  Complete  HRR  curve,  Test  3,  generic  workstation,  no  Jet  A,  30  percent  tile 
coverage  on  exposed,  upward-facing  surfaces. 
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Figure  3-3b.  Annotated  partial  HRR  curve,  Test  3,  complete  generic  workstation, 
no  Jet  A,  30  percent  tile  coverage  on  exposed,  upward-facing  surfaces. 
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Time  from  Ignition  (s) 
Figure  3-4a.  Complete  HRR  curve,  Test  5,  generic  workstation,  4  L  Jet  A,  no  tiles. 
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Figure  3-4b.  Annotated  partial  HRR  curve,  Test  5,  complete  generic  workstation, 

4  L  Jet  A,  no  tiles. 
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Figure  3-5a.  Complete  HRR  curve,  Test  6,  generic  workstation,  4  L  Jet  A,  30  percent 
tile  coverage  on  exposed,  upward  facing  surfaces. 
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Figure  3-5b.  Annotated  partial  HRR  curve,  Test  6,  complete  generic  workstation, 
4  L  Jet  A,  30  percent  tile  coverage  on  exposed,  upward  facing  surfaces. 


22 


NISTNCSTAR  1-5C,  WTC  Investigation 


Test  Results  and  Interpretation 


Time  from  Ignition  (s) 
Figure  3-6b.  Annotated  partial  HRR  curve,  Test  4,  WTC  workstation,  no  Jet  A,  no  tiles. 
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Source:  NIST. 

Figure  3-7a.  2  MW  igniter  flame  at  time  zero. 


Source:  NIST. 

Figure  3-7b.  First  ignited  object  within  cubicle,  paper  stack  on  desk, 
at  approximately  60  s  after  start  of  2  MW  igniter  flame. 
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Source:  NIST. 


Figure  3-7d.  Collapsing  chair  fire  near  its  peak  at  512  s. 
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Source:  NIST. 

Figure  3-7e.  Work  station  fire  near  peak  HRR  at  533  s. 
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Table  3-1.  Key  results  from  the  workstation  fire  test  burns. 


Quantity 

Test 

1 

2 

3 

4 

5 

6 

Workstation 

'/2  Generic 

Generic 

Generic 

WTC 

Generic 

Generic 

Tiles 

N 

N 

Y 

N 

N 

Y 

Jet  fuel 

N 

N 

N 

N 

Y 

Y 

Peak  HRJl'(MW) 

5.92/5.77 

8.70/8.48 

7.56/7.30 

9.89/9.66 

9.12/8.91 

7.960/7.60 

Time  to  peak  (s) 

490 

530 

590 

510 

160 

200 

Net  peak  HRR'  (MW) 

3.82/3.67 

6.95/6.73 

5.53/5.27 

7.72/7.46 

7.38/7.17 

6.17/5.95 

Peak  MLR  (kg/s) 

0.197 

0.308 

0.263 

0.420 

0.336 

0.293 

Time  to  peak  (s) 

J  JU 

JOU 

/ion 

1  Af> 

1  oU 

i  .zu 

A  1  1 

1  An 

1  lA 

X                illlCI  Veil        S  J 

1  JU  10  ZDJ 

jU  to  jZUU 

1  An  in  Ann 
1  ou  to  ouu 

on  t^-v  0  1  nn 
JU  to  L  J  uu 

n  (<->  o'^fin 
U  to  ZjUU 

Tn  in  o'son 

ZU  lO  ZjZU 

69.1 

205.0 

213.6 

173.6 

200.2 

205.3 

Effective  heat  of 
combustion  (MJ/kg  ) 

17.4 

19.8 

19.3 

16.9 

1 8.0 

1 8.Z 

FWHH'^  (s) 

244 

445 

318 

451 

t(75  %y 

1311 

1453 

833 

1009 

t^g  (s  )  (item  ignited) 

39 
(paper) 

67 
(paper) 

56 
(paper) 

50 
(paper) 

90 
(Jet  A) 

114 
(paper) 

a.  The  first  number  is  the  calorimeter  output;  the  second  is  a  10s  average  about  the  absolute  peak. 

b.  The  time  interval  applies  to  both  the  net  heat  released  and  to  the  total  mass  loss. 

c.  There  was  some  spillage  of  smoke  in  Test  4,  which  may  partly  account  for  the  lower  heat  of  combustion. 

d.  Full  width  half  height  of  net  HRR  curve. 

e.  Time  at  which  75  percent  heat  had  been  released  and  75  percent  of  mass  had  been  lost, 
f  Time  of  ignition  of  first  object  within  workstation. 


3.1.2        Heat  Release  Rate 


The  HRR  was  computed  from  the  measured  species  concentrations  and  flow  measurements  using  the 
following  equation: 


L  An 


(1  ~  ^//,o 


Mo. 


(3-1) 


l  +  ^(a-l)  -        -  M 


air 


The  oxygen  depletion  factor,  (j) ,  was  expressed  in  terms  of  the  volume  fractions  of  the  species  in  the  stack 
and  the  volume  fraction  in  the  ambient  air  (indicated  by  a  superscript ")  as  follows: 

,  (1  ~  ^co.  ~  ^co )  ~  ^o,  (1  ~  ^co. ) 

H>  =  —  (3-2) 
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The  terms  in  the  equation  are  defined  below.  For  the  constant  parameters,  the  numerical  values  used  in 
the  calculations  were: 

'    ( A//^,         ^  =  heat  of  combustion  of  hydrocarbon  fuel,  MJ/kg  O2  ( 1 3. 1  MJ/kg) 
(A//^         ^  =  heat  of  combustion  of  carbon  monoxide,  MJ/kg  O2  (17.6  MJ/kg) 
=  mass  flow  rate  in  exhaust  duct,  kg/s 

a  =  combustion  products  expansion  factor  (1.10) 

M,  =  molecular  weight  of  gas  i,  g/mole  (28.967  for  air) 

3.1.3        Net  Heat  Release  Rate 

Comparison  of  the  HRR  curves  with  the  direct  igniter  contribution  subtracted  allowed  more  appropriate 
comparison  of  the  various  test  conditions.  This  net  HRR,  H  ,  was  computed  using  the  formula: 

H"{t)^H{t)-H,{t)  ,  t<te. 

...  (3-3) 
H"{t)^H{t)-H,^(t)        ,  t  >tex 

where  H(t)  is  the  total  heat  release  measured  by  the  calorimeter,  H^(t)is  the  heat  released  by  the  liquid 
spray  burner,  and  (t) ,  the  background  heat  release  rate,  and  tex  is  the  time  the  spray  burner  is  turned 
off  The  middle  two  tenns  are  explained  as  follows: 

•  The  liquid  burner  had  a  approximately  2  MW  HRR  while  operating.  The  HRR  of  the  burner 
varied  by  as  much  as  20  percent  by  the  time  the  burner  was  shut  off.  The  expression 

for //y  (t)  varied  linearly  in  time  from  the  initial  burner  HRR  value  to  the  final  value. 

•  The  HRR  is  proportional  to  the  change  in  the  O2  volume  fraction  in  the  stack.  A  drift  in  the 
oxygen  analyzer  will  appear  as  a  HRR  even  though  there  is  no  fire.  Such  drifts  were  small 
and  were  represented  by  a  function  for (/)  that  varied  linearly  in  time  from  the  apparent 
heat  release  before  the  test,  which  was  very  small,  to  the  value  at  the  end  of  the  test  when  the 
fire  was  out. 

Figures  3-8  and  3-9  show  the  plots  of  net  HRR  versus  time. 
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3.1.4        Net  Heat  Released 

The  net  heat  released  by  the  workstation,  H,  is  obtained  by  integrating  the  net  HRR  over  the  time  of  the 
bum. 


Using  a  standard  plotting  routine,  the  mass  loss  during  each  test  was  plotted  as  shown  in  Figs.  3-10  and 
3-11.  The  smooth  S-shaped  curves  are  occasionally  interrupted  by  abrupt  steps  up  or  down,  the  result  of 
something  falling  off  the  weighing  platform  (usually  one  of  the  privacy  panels)  and  perhaps  being 
replaced. 

The  mass  loss  rate  (MLR)  at  time  t,  Ml(t),  was  computed  as  the  slope  from  the  best  linear  fit  to  the  data 
based  on  a  fixed  number  of  points  on  each  side  of  the  datum  point  of  interest.  The  best  balance  between 
noise  suppression  and  degraded  time  response  was  obtained  for  a  linear  fit  of  1 5  points  on  each  side  of  the 
datum  points.  Using  5  points  or  even  10  points  would  produce  significant  noise  artifacts.  The  use  of 
15  points  rather  than  a  smaller  number  of  points  had  less  than  a  2  percent  effect  on  the  peak  height  or  the 
time  of  the  peak  MLR.  Other  methods  for  computing  MLR,  including  the  direct  calculation  of  Am/At  for 
time  intervals  of  10  s,  20  s,  and  30  s  and  the  use  of  a  numerical  differentiation  algorithm,  were  tried,  but 
they  produced  a  slightly  noisier  fit  for  the  same  averaging  interval. 

The  MLR  histories  were  then  compared  with  the  net  heat  release  rate  curves,  as  indicated  in  Figs.  3-12 
through  3-16.  The  abrupt  steps  in  the  mass  loss  arising  from  objects  falling  off  the  load  cell  appear  as 
large  discontinuities  in  the  MLR  in  these  figures.  For  ease  of  viewing  these  have  been  removed  from  all 
but  Fig.  3-12,  where  they  have  been  left  for  demonstration  purposes. 

For  a  single,  homogeneous  material,  the  MLR  would  follow  the  HRR,  with  the  proportionality  being  the 
heat  of  combustion.  Here  the  workstation  includes  a  mix  of  materials  with  differing  heats  of  combustion 
and  differing,  time-dependent,  individual  MLR  s  that  sum  to  the  whole  and  the  correlation  between  the 
two  need  not  be  rigidly  fixed.  Nonetheless,  the  agreement  is  surprisingly  good,  reflecting  the  concurrent 
burning  of  multiple  cubicle  components  after  the  initial  buildup. 


(3^) 


where  /„  is  the  time  the  spray  burner  is  turned  off  and  t)  is  the  end  of  the  test. 


3.1.5 


Mass  Loss  and  Mass  Loss  Rate 
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Mass  Loss  Rate  (kg/s) 
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h.tass  Loss  Rate  (kg/s) 
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Mass  Loss  Rate  (kg/s) 
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Mass  Loss  Rate  (kg/s) 


in 

o 


CO 

o 


CD 
CD 


IT) 


o 
o 


o 
o 


o 

CD 
CD 


CD 

o 


o 
IE 

!_ 

E 


£ 
P 


(A 
(0 


c 

0)  ^ 
(A  0) 

ni  3 

0) 

TO  W 
0)  0) 

lo 

c  w 
o  a> 

II 
=  S 

o 

0) 

o 
m 

(A 

50  O 

I 

(O 

0) 
l_ 

3  . 

O) 


OD 


38 


NISTNCSTAR  1-5C,  WTC  Investigation 


Test  Results  and  Interpretation 


3.1.6        Effective  Heat  of  Combustion 

Two  methods  were  used  for  determining  the  effective  heat  of  combustion.  The  first  was  to  compute  the 
total  heat  release  and  the  total  mass  loss  over  an  interval  of  time  and  then  take  the  ratio.  These  results  are 
included  in  Table  3-1 .  The  second  method  was  to  compute  the  ratio  of  the  net  HRR  to  the  MLR  at  each 
second  of  time.  Before  computing  the  ratio,  a  10  s  smoothing  of  the  data  was  applied.  For  each  datum, 
the  five  points  prior  and  the  five  points  after  were  averaged  with  the  datum.  As  expected,  this  differential 
approach  varied  widely  with  time  but  did  provide  insight  on  the  time  dependence  of  the  heat  of 
combustion. 

Figure  3-17  shows  the  results  are  noisy  due  to  taking  the  ratio  of  two  quantities,  each  with  appreciable 
noise.  It  also  shows  some  time  variation,  starting  somewhat  below  20  MJ/kg  and  increasing  gradually  to 
roughly  25  MJ/kg.  The  discontinuity  around  1500  s  is  the  result  of  a  panel  falling  off  the  weighing 
platform.  The  heat  of  combustion  was  separately  computed  based  on  the  net  heat  released  during  the 
intervals  over  which  the  mass  decreased  0  kg  to  50  kg  and  100  kg  to  150  kg.  The  heat  of  combustion 
increased  from  18.4  MJ/kg  to  22.6  MJ/kg  over  these  intervals.  Part  of  the  reason  for  the  low  initial  value 
for  the  plotted  heat  of  combustion  is  a  slight  offset  between  the  MLR  and  HRR.  Early  in  the  fire 
development  there  may  be  pyrolysis  gases  leaving  the  workstation  without  combusting.  Late  in  the  fire, 
when  the  effective  heat  of  combustion  has  risen,  the  only  things  burning  are  charred  cellulosic  materials 
(wood  and  paper).  Such  a  rise  is  seen  only  at  the  end  of  the  cone  calorimeter  tests  with  these  materials. 


400  800  1200  1600 

Time  from  ignition  (s) 


Figure  3-17.  Heat  of  combustion  vs.  time  for  Test  2  (no  jet  A,  no  tiles). 
3.2         FIRE  GROWTH  BEHAVIOR 

The  following  text  presents  observations  of  the  progress  of  the  fires  in  the  set  of  tests.  The  purpose  is  to 
characterize  the  phenomena  that  controlled  the  spread  of  flame  on  workstation  components  and  the  spread 
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from  component  to  component.  Some  of  these  phenomena  are  not  included,  or  even  capable  of  being 
included,  in  a  computational  model  of  current  sophistication.  In  those  cases,  an  effort  was  made  to 
estimate  whether  the  phenomenon  was  likely  to  be  a  significant  contributor  to  the  shape  or  magnitude  of 
the  HRR  curve  over  time.  These  estimates  were  then  considered  when  evaluating  the  viability  of  FDS  for 
simulating  the  thermal  environment  from  the  combustion  of  fuel  packages  such  as  those  tested  here. 

3.2.1  Role  of  the  Ignition  Burner 

The  HRR  curves  all  began  with  an  abrupt  step  up  to  approximately  2  MW  when  the  heptane  burner  was 
ignited.  There  was  an  immediate  bathing  of  the  cubicle  surfaces  with  radiant  energy  from  this  flame,  as 
indicated  by  a  floor  flux  gage  reading  of  5  kW/m"^  in  Test  1 .  This  irradiance  reaching  workstation 
surfaces  from  this  source  was  greatest  early  in  the  test,  before  the  workstation  reached  its  peak  burning 
rate.  The  later  falloff  was  due  to  (a)  the  workstation  surfaces  eventually  collapsing  downward  away  from 
the  igniter  flames  and  (b)  the  workstation  flames  being  nearly  optically  thick  (and  thus  opaque  to  the 
igniter  radiation)  near  the  peak  burning  stage. 

Generally,  the  ignition  burner  was  turned  off  well  past  the  peak  HRR  was  reached.  By  then,  the 
combustion  of  the  collapsed  pile  of  burning  debris  was  minimally  dependent  on  the  extra  radiation.  In 
Test  4,  the  burner  was  turned  off  near  the  peak  HRR  because  the  hood  capacity  could  be  exceeded.  Some 
insight  into  the  effect  of  turning  off  the  igniter  early  can  be  gained  from  examining  the  net  HRR  curves 
from  Test  2  (generic  workstation)  and  Test  4  (WTC  workstation)  (Fig.  3-14).  In  Test  2,  the  burner  was 
kept  on  for  at  least  2000  s;  in  Test  4,  the  burner  was  turned  off  at  approximately  575  s,  just  after  the  HRR 
peak.  The  burner  did  not  actually  turn  off  abruptly,  and  a  sharp  decrease  in  the  Test  4  HRR  can  be  seen  in 
the  2  min  following  the  act  of  turning  off  the  fuel.  The  HRR  curve  for  Test  4  drops  significantly  below 
the  curve  for  Test  2.  This  difference  in  area  between  the  two  curves  was  larger  than  can  be  explained  by 
the  small  differences  in  combustible  mass  and  material  composition  between  the  two  workstations.  It 
more  likely  reflected  the  loss  of  the  boosting  radiation  from  the  burner. 

3.2.2  Growth  to  the  HRR  Peak  without  Jet  Fuel  Present 

In  Tests  1  through  4,  there  was  a  fairly  reproducible  progression  of  ignition  events  leading  up  to  the  HRR 
peak.  These  are  typified  by  the  sequence  of  photos  from  Test  2,  starting  with  the  source  fire  at  time  zero 
(Fig.  3-7a),  through  the  first  ignition  of  an  item  within  the  cubicle  (Fig.  3-7b),  some  secondary  ignitions 
(Fig.  3-7c),  reaching  the  peak  HRR  of  the  collapsing  chair  (Fig.  3-7d),  then  the  fire  peak  itself 
(Fig.  3-7e)  and,  finally,  a  late  stage  with  some  of  the  surfaces  collapsed  (Fig.  3-7f)- 

The  igniter  flame  bathed  the  entire  workstation  in  radiant  heat.  FDS  simulations  suggest  that  the  peak 
fluxes  on  the  side  of  the  desk  top  toward  the  igniter  were  of  the  order  of  30  kW/m"  or  more.  This  was 
well  above  the  minimum  flux  for  piloted  ignition"  of  the  various  exposed  surfaces  there  (Table  A-2). 
However,  the  visual  evidence  indicated  that  the  igniter  flame,  despite  its  size,  was  not  able  to  fimction  as 
a  pilot  for  any  of  the  workstation  components.  The  igniter  was  too  removed  from  the  irradiated  fuels  to 
be  effective.  Rather,  non-piloted  ignition,  or  autoignition,  or  pyrolysis  vapors  appeared  to  be  the 
initiating  event.  This  is  despite  the  fact  that  the  minimum  radiant  flux  for  autoignition  of  solid  materials 


"  Piloted  ignition  means  that  there  is,  in  the  gas  phase  very  near  the  igniting  object,  some  heat  source  sufficient  to  initiate  a  flame 
in  the  gases  evolving  from  a  fuel  surface  as  soon  as  those  gases  form  a  flammable  mixture  with  ambient  air. 
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tends  to  be  much  higher  than  that  for  piloted  ignition.  Boonmee  and  Quintiere  (2001),  for  example,  found 
that  a  bare  wood  surface  required  a  radiant  flux  of  about  40  kW/nr  for  non-piloted  ignition.  Other 
materials  would  be  more  susceptible  and  within  the  range  of  the  irradianccs  imposed  by  the  spray  burner. 

•  While  the  fabric  might  be  expected  to  ignite  first  since  a  portion  of  it  was  initially  within 
approximately  5  cm  (2  in.)  of  the  igniter,  it  melted  and  curled  downward  toward  the  rear  edge 
of  the  desk  top,  igniting  only  later. 

•  The  top  of  the  computer  monitor  shell,  which  gasified  extensively  only  20  cm  (8  in.)  from  the 
spray  burner  flames,  did  not  appear  to  ever  have  been  directly  ignited  by  those  flames. 

•  Instead,  the  progression  of  flaming  was  initiated  by  the  top  sheets  on  the  paper  stacks,  mid- 
length  along  the  side  nearest  the  spray  burner.  These  sheets  are  able  to  autoignite  by 
charring,  curling  upward  to  escape  the  heat  sinking  and  oxygen  blocking  effects  of  the  lower 
sheets,  glowing,  and  then  transitioning  into  flaming  combustion.  The  paper  was  the  only 
material  present  in  the  workstations  that  exhibited  such  an  autoignition  process.  Because  the 
paper  was  present  in  many  places,  it  was  highly  instrumental  in  spreading  the  flames  over  the 
area  at  and  above  the  desktop. 

The  flames  on  the  paper  stacks  then  piloted  the  ignition  of  the  computer  monitor  shell,  the  desktop,  and 
the  top  of  the  chair  back.  Later,  paper  ignition  brought  flames  to  the  two  other  sides  of  the  desk  area  in 
the  workstation.  Once  ignited  by  the  paper,  the  chair  fire  spread  from  the  top  of  the  back  to  the  arms  and 
seat.  The  fabric  on  the  chairs  used  in  the  basic  workstation  was  a  char  fonner,  possibly  cotton-based. 
Such  a  fabric  might  autoignite  via  a  glowing  char  oxidation  reaction.  However,  the  paper  seemed  to  pilot 
the  chair  top  ignition  in  cases  where  that  process  could  be  observed.  Invariably,  some  flaming  material 
would  then  fall  from  the  burning  arm  or  seat  onto  the  carpet,  providing  a  local  pilot  there.  The 
subsequent  spread  on  the  carpet  was  delayed  until  the  chair  was  fiilly  involved,  along  with  the  remainder 
of  the  upward-facing  surfaces  of  the  desktop. 

The  inward-facing  surfaces  of  the  wall  panels  were  not  important  here  since  the  only  combustible  external 
component  on  the  panels  of  the  generic  workstation  was  the  thermoplastic  fabric.  This  fabric  consistently 
melted  and  rolled  downward  into  a  mass  having  much  less  surface  area  than  the  original  fabric  prior  to 
igniting.  While  it  undoubtedly  was  burning  at  the  HRR  peak,  its  drastically  reduced  surface  area  would 
make  its  contribution  negligible.  The  fabric  on  the  actual  WTC  workstation  behaved  similarly.  This  type 
of  behavior  is  not  captured  in  cone  calorimeter  tests. 

The  behavior  of  the  chair  was  critical  to  the  rapid  engulfment  of  the  lower  compartment  in  flames.  As 
noted  previously,  the  chair  was  largely  composed  of  thermoplastic  materials.  When  only  the  upholstered 
surfaces  of  its  seat  and  back  were  burning  it  retained  its  original  shape  and  little  of  its  heat  reached  the 
lower  compartment  area.  One  comer  of  the  seat  was  deliberately  placed  approximately  15  cm  (6  in.) 
under  the  desk  surface,  under  the  area  on  which  the  computer  keyboard  rested.  This  was  to  promote  some 
heat  reaching  this  space  early  in  the  fire;  however,  the  flames  contacted  the  relatively  ignition  resistant 
underside  of  the  desk.  Only  as  the  support  shell  of  the  chair  began  to  melt  and  flow  to  the  floor  did  the 
flaming  chair  approach  its  peak  burning  rate.  Because  the  chair  had  a  partial  steel  skeleton,  the  burning 
process,  at  its  maximum,  was  three  dimensional,  extending  from  desk  level  to  floor  level  and  involving 
the  carpet  under  the  chair,  as  well.  It  radiated  strongly  into  the  space  beneath  the  desk  as  well  as  directing 
an  increased  portion  of  its  fire  plume  into  this  space.  This  quickly  ignited  the  carpet  and  underside  of  the 
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desk  (plus  the  wastebasket  and  the  combustible  face  of  the  computer  processor  enclosure)  leading  directly 
to  the  rapid  upward  step  in  HRR  immediately  prior  to  the  peak  seen  in  Figs.  3-1,  3-2,  3-3,  and  3-6.  This 
burning  behavior  of  the  chair  cannot  be  inferred  from  the  cone  calorimeter  data  obtained  on  the  chair 
materials  (Appendix  A).  Empirical  treatment  of  its  HRR  process  is  necessary. 

The  progress  of  the  ignited  workstation  toward  the  peak  HRR  behaved  as  if  composed  of  two  zones.  The 
first  was  the  entire  volume  at  and  above  the  top  of  the  desk  surface,  capped  on  top  by  the  false  ceiling 
above  the  workstation.  Typically,  this  zone  became  fully  involved  first  in  the  sense  that  all  combustible 
surfaces  ignited.  The  second  zone  was  that  volume  at  and  below  the  underside  surface  of  the  desk.  It 
reached  a  point  of  local  flashover  when  subjected  to  the  radiant  heat  flux  from  the  upper  zone  and  hot 
ceiling  layer  but  especially  the  heat  from  the  burning  chair.  When  all  flammable  surfaces  in  both  zones 
were  ignited,  the  HRR  quickly  peaked. 

Madrzykowski  (1996,  1998)  also  noted  this  correspondence  of  the  peak  HRR  and  the  simultaneous 
burning  of  all  the  accessible  combustible  surfaces  in  the  workstation  interior.  He  attempted  to  combine 
this  type  of  observation  with  cone  calorimeter  HRR  measurements  and  the  known  areas  of  the  different 
fuels  that  were  burning  to  predict  the  peak  HRR  value  with  limited  success.  The  difficulty  in  doing  this 
was  that  these  various  surfaces  saw  a  range  of  heat  fluxes  and  thus  there  was  no  firm  basis  for  applying 
the  cone  calorimeter  HRR  data,  which  are  sensitive  to  the  irradiance. 

3.2.3        Effect  of  Jet  Fuel  on  the  Growth  to  the  HRR  Peak 

The  presence  of  Jet  A  might  be  expected  to  short  circuit  this  whole  sequence  and  lead  to  an  almost 
explosive  engulfment  of  the  entire  workstation.  Indeed,  the  time  to  the  peak  HRR  was  shortened 
considerably.  However,  there  was  a  delay  in  the  ignition  of  the  cubicle  combustibles. 

The  flash  point  of  Jet  A  is  in  the  range  of  43°  C  to  66°  C;  thus,  there  was  no  flammable  mixture  of  vapor 
above  it  at  the  beginning  of  the  test  when  the  cubicle  was  at  ambient  temperature.  As  the  test  proceeded, 
the  strong  radiant  flux  from  the  spray  burner  quickly  overcame  this  deficit  in  temperature,  in  spite  of  the 
cooling  effect  that  the  vaporization  process  provided,  and  the  videos  showed  an  increasingly  dense 
aerosol  rising  from  the  various  wetted  surfaces  near  the  spray  burner.  However,  this  aerosol  did  not 
ignite;  evidently  turbulent  mixing  of  the  vapor  plume  above  the  surfaces  quickly  diluted  it  to  below  the 
flammable  limit.  Rather,  in  Test  5,  a  random  piece  of  flaming  material  rose  from  the  burner  area  and 
drifted  down  on  top  of  the  desk,  immediately  igniting  the  Jet  A.  A  flaming  piece  of  debris  also  arose 
from  the  burner  in  Test  6,  but  it  landed  on  top  of  the  bookcase  where  there  was  no  Jet  A,  and  it  had  no 
effect.  In  Test  6,  the  tops  of  the  paper  stacks  mid-length  along  the  side  next  to  the  spray  burner  finally 
dried  out  and  began  to  char  and  then  transitioned  into  flaming,  igniting  the  Jet  A.  Thus,  in  this  case,  the 
presence  of  Jet  A  and  the  evaporative  cooling  effect  it  provided  pushed  this  paper  auto-ignition 
mechanism  back  from  the  40  s  to  70  s  seen  in  "dry"  tests  to  150  s.  Presumably,  more  Jet  A  could  have 
delayed  this  mechanism  longer.  Of  course,  the  presence  of  any  other  type  of  pilot  flame  close  to  the 
wetted  surfaces  could  have  greatly  shortened  the  ignition  delay. 

This  was  a  strong  indication  of  the  difficulty  of  predicting  ignition  in  many  room  fire  growth  situations. 
Even  a  fuel  as  volatile  as  Jet  A  yielded  only  a  relatively  small  spatial  zone  in  which  a  flammable  vapor 
mixture  existed,  and  it  was  difficult  for  nearby  flames  to  reach  that  zone  and  ignite  that  mixture  before  air 
turbulence  diluted  it.  In  the  absence  of  a  char  that  got  much  hotter  than  the  fuel  surface  and  thereby 
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autoignited  the  fuel,  ignition  of  the  fuel  vapor  might  be  delayed  indefinitely.  Use  of  the  cone  calorimeter 
to  obtain  ignition  delay  time  data  for  use  in  a  fire  growth  model  could  thus  be  seriously  misleading  in 
some  cases  since  the  pilot  is  placed  directly  in  the  pyrolyzatc  plume  in  order  to  produce  ignition  at  the 
earliest  possible  time. 

Once  the  Jet  A  ignited  locally,  the  flames  spread  rapidly  but  not  continuously  along  the  desktop. 
Presumably,  the  initial  evaporation  period  left  some  dry  spots  that  interrupted  the  spread.  The  same 
turbulent  vapor  plume  dilution  process  discussed  above  prevented  flames  from  jumping  downward  to  ihc 
carpet  from  the  flaming  desk  surface,  even  though  the  carpet  was  emitting  a  visible  aerosol.  The  carpet 
ignited  when,  as  in  the  "dry"  tests,  the  chair  fire  dropped  a  flaming  pilot  to  the  floor.  Here,  however,  the 
subsequent  flame  spread  on  the  carpet  was  rapid.  None  of  these  steps  was  distinguishable  in  the  HRR 
curves  of  Figs.  3-5  and  3-6;  the  HRR  curves  jumped  rapidly  to  their  peak  values  soon  after  the  Jet  A 
ignited  since  the  time  span  over  which  the  separate  steps  occurred  was  compressed.  Here  again,  the  peak 
corresponded  to  all  accessible  combustible  surfaces  burning  simultaneously. 

3.2.4        Post-Peak  Behavior 

The  rather  rapid  decay  in  HRR  after  all  of  the  peaks  reflects  several  factors. 

•  The  chair  fire  rapidly  collapsed  to  a  pool  fire  on  the  floor,  and  the  reduced  burning  surface 
area  meant  a  reduced  HRR. 

•  The  various  paper  piles  developed  a  thick  ash  layer  that  drove  down  their  burning  rate. 

•  Similarly,  char  formation  on  the  desk  surfaces  drove  down  their  burning  rates. 

•  The  carpet  began  to  bum  out. 

As  the  various  cubicle  components  became  physically  weakened  by  the  fire,  their  geometric  arrangement 
changed.  The  computer  monitor  shell  was  the  first  item  to  change  shape,  with  the  top  of  the  shell  melted 
and  draped  down  over  the  CRT  even  before  ignition.  However,  its  peak  HRR  was  probably  less  than 
200  kW.  Next  came  the  chair  melt  down,  as  mentioned  above.  The  front  of  the  bookcase,  resting  on  top 
of  that  unit,  typically  fell  next.  The  fates  of  the  bookcase  itself  and  its  1 3  reams  of  paper  were  much  more 
difficult  to  discern  through  the  flames.  It  sometimes  appeared  to  come  down  at  the  same  time  as  the 
bookcase  front  panel  and  sometimes  it  came  down  distinctly  later.  The  desk  surface  bowed  progressively 
as  it  charred  through  and  then  it  collapsed,  with  separate  sections  doing  so  at  differing  times.  The  initial 
desk  collapse  did  little  to  change  the  available  burning  surface  area  and  thus  probably  did  not  greatly 
affect  its  burning  rate.  However,  ultimately  what  was  left  was  a  complex  rubble  pile.  Similarly,  the  wall 
panels  collapsed  at  random  times,  both  inward  and  outward,  typically  rather  late  in  the  fire. 

These  changes  led  to  departures  from  the  burning  rates  based  on  the  original  surface  areas  coupled  with 
cone  calorimeter  data.  The  further  past  the  peak  of  the  HRR  curve,  the  less  the  burning  rate  was 
predictable  by  an  FDS  calculation  that  retained  the  original  cubicle  geometry.  Fortunately,  the  most 
serious  disagreements  (other  than  the  complex  chair  collapse)  should  not  begin  until  well  after  the  desk 
surfaces  collapse. 
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In  all  tests,  the  paper  stacks  both  in  the  file  cabinets  and  on  the  desktop  above  the  file  cabinets  were 
comparatively  undisturbed  by  the  collapse  of  the  workstation.  The  paper  was  one  of  the  longer  flaming 
items,  since  it  formed  a  thick  layer  of  insulating  residue  that  slowed  the  involvement  of  the  deeper  layers. 
The  paper  stacks  on  the  desktop  appeared  to  roughly  double  in  thickness  as  a  result  of  this  residue 
formation.  In  cases  noted  below,  measurements  were  made  of  the  residual  mass  of  these  stacks  after  the 
bulk  of  the  fire  had  been  extinguished.  Often  these  paper  stacks  continued  to  smolder  right  up  to  the  time 
they  were  weighed  (i.e.,  they  were  continuing  to  lose  weight  slowly). 

An  implication  of  this  behavior  is  that  the  tail  of  the  fire  would  be  expected  to  be  more  intense  and 
continue  for  a  longer  time  in  a  workstation  having  more  paper  than  that  used  here.  In  particular,  these 
same  workstations  might  have  had  a  longer  fire  tail  had  their  file  cabinets  been  full  of  paper.  The  greater 
mass  combined  with  the  more  restricted  air  supply  in  a  filled  file  cabinet  could  be  expected  to  slow  the 
burning  process.  Relative  to  the  generic  workstation,  the  WTC  workstation  had  a  file  cabinet  design  that 
allowed  more  air  access  to  the  filled  file  drawers,  once  the  drawer  facing  burned  off.  This  could  possibly 
have  supported  more  extensive  burning  of  paper. 

The  consumed  mass  fractions  of  paper  were  as  follows  (averaged  from  replicate  measurements  in  Tests  2, 
3,  and  4,  none  of  which  involved  Jet  A). 

Paper  on  desk:  mass  fraction  burned  =  0.52,  0.68 

Paper  in  file  cabinets:  mass  fraction  burned  =  0.31,  0.44,  0.30 

The  higher  mass  fraction  of  file  drawer  paper  consumption  occurred  with  the  WTC  workstation.  The 
paper  stacks  on  the  desk  were  those  in  the  southwest  comer  of  the  cubicle  and  involved  three  side-by-side 
horizontal  layers,  each  one  ream  thick.  In  all  cases  there  was  some  uncharred  paper  in  the  center  of  the 
stacks.  Since  paper  is  a  thermally  thin  material  with  a  large  surface  to  volume  ratio,  it  has  the  potential  to 
bum  very  rapidly  if  it  is  dispersed  so  as  to  expose  this  large  surface.  However,  no  such  dispersal 
happened  here.  Instead,  the  compactness  of  the  paper  stacks  led  to  an  opposite  extreme  of  behavior  and 
incomplete  combustion. 

3.3  PEAK  HRR  VALUES  AND  TIMES  TO  PEAK  VALUES 

Table  3-1  lists  the  peak  HRR  values  and  the  times  at  which  they  were  reached.  The  HRR  peak  is  given  in 
two  ways:  ( 1 )  as  the  absolute  highest  single  reading  recorded  (for  a  1  s  interval)  and  (2)  as  the  average  of 
the  value  at  this  peak  plus  values  up  to  5  s  to  either  side  of  this  peak  (total  of  1 1  data  points).  The  latter 
compensated  for  both  noise  in  the  calorimeter  system  and  in  the  fire  itself.  Use  of  these  damped  values 
was  appropriate  since  the  principal  combustibles  in  a  serious  fire  have  thermal  inertias  that  cannot 
respond  to  small  HRR  fluctuations  that  are  under  10  s  in  duration.  The  damped  values  were  all  lower 
than  the  absolute  peak  values  (by  approximately  2  percent  to  4  percent)  due  to  the  sharpness  of  the  peaks. 
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A  few  points  emerged  from  Table  3-1  when  the  various  numbers  were  compared. 

•  The  half  workstation  yielded  a  peak  fire  intensity  that  was  about  two-thirds  that  of  the  lull 
workstation.  This  was  primarily  the  result  of  two  factors: 

-  The  same  chair  was  present  in  both  cases;  this  chair  had  an  estimated  HRR  peak  of 
Yi  MW  by  itself 

-  Twice  as  much  of  the  carpet  in  the  second  "half  of  the  workstation  was  covered  by  inert 
steel  file  cabinets,  precluding  its  participation  in  the  HRR  peak.  This  lowered  the  HRR 
contribution  from  the  second  half  of  the  workstation. 

The  computer  was  also  totally  within  the  half  workstation,  though  this  was  a  substantially  lesser  heat 
source  than  was  the  chair. 

•  Coverage  of  30  percent  of  the  upward  facing  horizontal  surfaces  by  inert  tiles  slightly 
increased  the  time  to  the  HRR  peak  and  decreased  the  peak  HRR  value  by  only 
approximately  14  percent.  These  tiles  effectively  prevented  direct  burning  of  the  surfaces  on 
which  they  lay.  As  described  above,  the  overall  HRR  peak  occurred  when  all  accessible 
flammable  surfaces  were  burning.  The  tiles  did  not  cover  any  fraction  of  such  surfaces  as  the 
underside  of  the  desk,  the  carpeting  below  the  desk,  or  the  underside  and  back  of  the  chair. 
Their  percent  reduction  in  HRR  was,  in  fact,  roughly  comparable  to  their  percent  coverage  of 
the  total  accessible  combustible  surface  area  of  the  workstation  (approximately  20  percent). 
Thus,  their  effect  on  the  peak  fire  size  was  roughly  linear  in  fractional  area  coverage.  Exact 
linearity  was  not  expected  since  differing  materials  have  differing  HRR  per  unit  area 
contributions.  The  response  might  have  become  strongly  non-linear  (in  the  direction  of 
greater  suppressive  effect)  in  a  more  extreme  case  where  nearly  all  of  the  upward  facing 
horizontal  surfaces  were  covered,  since  this  could  preclude  the  progressive  flame  spread  that 
gets  all  accessible  surfaces  involved. 

•  The  presence  of  4  L  of  Jet  A  led  to  a  dramatic  decrease  in  the  times  to  peak  HRR,  from 
approximately  10  min  to  approximately  3  min.  This  was  the  result  of  the  Jet  A  shortening  the 
time  to  involvement  of  all  accessible  combustible  surfaces.  The  peak  itself  was  boosted 
upward  by  only  approximately  5  percent,  presumably  because  the  Jet  A  helped  boost  the 
overall  fuel  gasification  rate  somewhat  while  adding  its  high  heat  of  combustion.  For 
comparison,  a  pool  of  pure  Jet  A  having  an  area  equal  to  that  of  the  plan  view  area  of  the 
cubicle  (5.9  m',  64  ft')  would  have  yielded  a  steady  heat  release  rate  of  13.1  MW  (Mudan  and 
Croce  1995).  Though  Jet  A  covered  a  comparable  area  here,  it  was  a  very  thin  layer,  either 
on  the  high-heat-sink  desk  surface  or  soaked  into  the  carpet;  this  kept  its  burning  rate  down. 
When  the  inert  tiles  were  also  present,  the  Jet  A  was  poured  across  their  top  surfaces, 
temporarily  rendering  these  surfaces  flammable.  Because  the  tiles  were  porous,  the  Jet  A 
burning  rate  on  them  was  reduced,  however,  and  the  tiles  still  managed  to  produce  an 
approximately  1 5  percent  reduction  in  the  peak  HRR  relative  to  the  case  with  Jet  A  and  no 
tiles. 
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Figure  3-18  shows  the  temporal  alignment  at  the  peak  of  various  measures  of  fire  intensity  in  Test  5. 
Each  of  the  measures  has  been  normahzed  by  its  maximum  value.  The  five  ineasures  are:  HRR,  MLR, 
carbon  dioxide  volume  fraction  in  the  exhaust  gases,  temperature  of  the  exhaust  gases,  and  radiant  flux 
seen  by  the  gage  viewing  the  entire  fire  from  the  south  side  of  the  cubicle.  The  radiometer  was  the  fastest 
responding  instrument  with  a  1  s  time  constant. 
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Figure  3-18.  A  comparison  of  five  measures  of  fire  intensity. 

The  five  measures  presented  a  consistent  representation  of  the  progress  of  the  fire,  given  their  individual 
measurement  characteristics  as  noted  in  the  following  bullets: 

•  HRR  and  CO:  concentration  rose  together  and  remained  close  throughout  the  test.  This  is 
reasonable  since  CO2  was  the  major  carbon-containing  product  of  combustion.  The  slight  dip 
in  the  post-peak  HRR  curve  relative  to  the  CO2  concentration  curve  is  within  the  uncertainty 
of  the  measurements.  However,  if  it  were  real,  it  probably  resulted  from  the  higher  H:C  ratio 
of  the  wood-based  fiiels  that  dominated  the  later  burning  compared  to  the  wood  and  plastics 
that  shared  the  earlier  burning.  The  combustion  of  the  hydrogen  atoms  to  fonn  water 
contributed  to  the  HRR  without  generating  CO2. 

•  These  two  peaks  lagged  the  radiant  flux  signal  peak  by  roughly  5  s,  though  the  noise  in  the 
latter  made  this  somewhat  ambiguous.  The  time  response  of  the  radiometer  was  about  0.2  s 
compared  to  about  a  10  s  response  time  of  the  calorimeter.  A  10  s  smoothing  of  the  peak 
radiometer  data  (5  data  points  before  and  after  the  point  being  smoothed)  was  used  to 
simulate  the  time  response  of  the  calorimeter.  The  ratio  of  the  smoothed  value  to  the  peak 
was  0.94  based  on  Tests  4,  5,  and  6  and  a  value  of  0.97  for  the  calorimeter.  This  result 
implied  that  the  calorimeter,  due  to  its  10s  time  response,  underestimated  the  peak  by  about 
3  percent.  In  principle,  if  the  radiant  flux  perfectly  reflected  the  total  HRR  of  the  fire  and  the 
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calorimeter  system  time  lags  were  fully  accounted  for,  the  HRR  peak  and  radiant  flux  peak 
should  have  coincided.  However,  perfect  coincidence  would  have  required  an  optically  thin 
fire  in  full  view  of  the  radiometer.  Here,  the  fire  was  approaching  optically  thick,  and  parts 
of  the  flames  were  blocked  from  radiometer  view  by  components  of  the  workstation.  Thus, 
in  practice,  there  is  no  guarantee  that  with  a  fire  of  this  size  and  spatial  complexity  the  radiant 
flux  to  an  external  point  completely  reflects  the  instantaneous  total  HRR  in  the  volume  it 
views.  It  was  also  possible  that  calorimeter  time  response  (in  addition  to  the  compensated 
system  flow  time  lag)  was  partly  responsible  for  the  apparent  time  lag  between  the  radiometer 
and  HRR  peaks.  In  the  case  of  the  exhaust  gas  temperature,  thermal  lags  in  the  duct,  due  to 
heat  loss  to  the  duct  walls,  were  the  source  of  the  slowed  response. 

•    The  rise  in  the  MLR  curve  preceded  the  HRR  curve  by  about  25  s  and  the  radiant  flux  curve 
by  about  15  s.  This  was  likely  a  result  of  mass  loss  by  evaporation  of  the  Jet  A  (and  other 
pyrolysis  processes)  that  yielded  a  net  mass  loss  in  excess  of  what  was  actually  burning. 

3.4  SIGNIFICANCE  OF  TEST  VARIABLES 

3.4.1        Effects  of  Surface  Coverage  and  Jet  Fuel 

Because  this  set  of  tests  included  a  full  factorial  design  for  the  effects  of  tile  coverage  and  Jet  A,  it  was 
possible  to  perform  a  sensitivity  analysis  of  how  the  various  burning  characteristics  of  the  workstation 
depend  on  these  two  factors  (Box,  Hunter,  and  Hunter  1978).  The  factorial  design  is  illustrated  in 
Fig.  3-19  with  a  box  for  which  the  lower  left  comer  corresponds  to  the  parameter  set  with  no  tiles  or 
Jet  A  (T-J-),  the  lower  right  to  tiles  and  no  Jet  A  (T+J-),  etc.  To  obtain  the  sensitivity  of  the  workstation 
fire  behavior  to  the  30  percent  tile  coverage  (of  upward  facing,  exposed  surfaces),  one  computes  the 
horizontal  differences  between  the  response  values  for  the  left  and  right  hand  sides.  The  response 
difference  values  shown  below  are  for  the  peak  HRR.  The  least  square  estimate  for  the  sensitivity  factor 
response  for  the  tile  coverage,  Pt,  is  identically  the  average  (-1.15)  of  these  two  differences  (-1.14  and 
-1.16)  normalized  by  the  average  (8.34)  of  the  data;  this  value  is  equal  to  -13.8  percent.  A  similar 
analysis  for  Jet  A  leads  to  Pj  =  4.9  percent.  The  interaction  sensitivity  response.  Pi,  is  equal  to  the 
difference  between  the  average  of  the  right  diagonal  terms  and  the  average  of  the  left  diagonal  terms 
divided  by  the  four  comer  average  and  is  equal  to  a  value  of  -0. 1  percent.  This  relatively  small  value  of 
the  interaction  sensitivity  indicates  that  the  tiles  and  Jet  A  are  nearly  independent  factors.  That  is,  the 
effect  of  the  tiles  on  the  burning  rate  of  the  cubicle  is  almost  independent  of  the  presence  or  absence  of 
the  Jet  A.  Note  that  the  sensitivity  factor  responses  calculated  here  are  for  the  specific  levels  of  change  in 
the  tile  coverage  or  Jet  A  factors  in  these  tests. 
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T-J+     M  ►  T+J++ 
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Figure  3-19.  Interaction  of  test  variables  in  determining  peak  HRR. 

Table  3-2  provides  a  complete  listing  of  the  sensitivity  factors  for  the  various  measured  quantities.  These 
are  the  percentage  changes  in  the  indicated  fire  quantities  resuhing  from  either  a  30  percent  coverage  of 
upward  facing  horizontal  surfaces  with  inert  tiles  or  from  adding  4  L  of  Jet  A  to  horizontal  surfaces  of  the 
generic  work  station.  The  interaction  sensitivity  measures  the  additional  effect  of  making  both  changes  at 
the  same  time. 


Tab 

e  3-2.  Sensitivity  analysis  factors. 

Quantity 

Mean 
Value 

pTi  Sensitivity' 
Factor  Response 
for  Tile  (%) 

Pj,  Sensitivity 
Factor  Response 
for  Jet  A  (%) 

Pi,  Interaction 
Sensitivity  Factor 
Response  (%) 

Instantaneous  peak  HRR  (MW) 

8.34 

-13.9 

5.0 

-0.1 

10  s  average  peak  HRR  (MW) 

8.11 

-14.4 

5.5 

0.2 

Time  to  instantaneous  peak 
HRR(s) 

370 

13.5 

102.7 

-2.7 

Net  instantaneous  peak  HRR  (MW) 

6.51 

-20.2 

8.2 

1.6 

Net  10  s  average  peak  HRR  (MW) 

6.28 

-21.3 

8.9 

1.9 

Peak  MLR  (kg/s) 

0.30 

-14.7 

9.7 

0.3 

Time  to  peak  MLR  (s) 

357 

7.0 

-104.9 

-1.4 

Net  heat  released  (GJ) 

3.88 

2.8 

-10.8 

0.8 

Total  mass  loss  (kg) 

206 

3.4 

-3.4 

-1.0 

Heat  of  combustion  (MJ/kg) 

18.8 

-0.8 

-7.7 

1.9 

It  is  seen  from  this  table  that  for  the  time  to  instantaneous  peak  HRR  response  and  for  the  time  to  peak 
MLR  response,  the  Jet  A  factor  had  the  largest  effect  (both  in  excess  of  100  percent ).  For  the  net 
instantaneous  peak  HRR  and  the  net  10  s  average  peak  HRR,  the  tile  factor  had  the  largest  effect 
(-20.2  percent  and  -21.3  percent,  respectively).  Also  it  is  noteworthy  that  there  was  less  than  a  5  percent 
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sensitivity  of  the  total  mass  loss  to  either  Jet  A  or  tiles.  The  importance  of  this  rcsuU  will  be  explained  in 
the  next  section. 

The  sensitivity  factors  were  based  on  experiments  with  30  percent  tile  coverage  of  upward  facing  surfaces 
and  4  L  of  Jet  A.  The  effect  of  smaller  coverage  and  a  smaller  amount  of  Jet  A  can  be  estimated  using  a 
linear  interpolation,  that  is,  a  15  percent  tile  coverage  would  result  in  sensitivity  factors  equal  to  half  the 
values  of  those  listed  in  Table  3-2.  This  can  only  be  extrapolated  to  larger  tile  coverage  or  Jet  A  quantity 
with  caution  because  it  is  empirical,  not  physical,  in  basis.  Furthermore,  it  applies  only  to  the  specific 
arrangement  and  materials  used  here. 

Another  important  observation  is  that  the  interaction  sensitivity  is  in  all  cases  less  than  3  percent  for  all 
10  parameters.  This  indicated  that  the  effect  of  the  coverage  factor  and  the  Jet  A  were  independent 
factors  in  determining  the  full  complement  of  flammabihty  parameters  in  Table  3-2.  This  finding  was 
exploited  in  assessing  the  uncertainty  regarding  the  repeatability  of  the  tests. 

3.4.2        Sensitivity  to  Loading  of  Jet  Fuel 

In  the  fires  in  the  twin  towers,  the  distribution  of  jet  fuel  was  likely  to  have  been  uneven.  It  is  also  likely 
that  the  burning  behavior  of  a  workstation  would  be  affected  by  the  initial  presence  of  differing  amounts 
of  jet  fuel.  This  study  examined  only  the  effect  of  the  addition  of  one  amount  of  jet  fuel  on  the  burning 
behavior  of  the  workstation.  This  level  was  comparable  to  the  uniform  coverage  of  the  five  affected 
floors  with  Jet  A  from  about  half  of  the  fijel  in  the  aircraft.  The  following  discussion  offers  some  insight 
into  the  effects  of  underestimating  or  overestimating  the  volume  of  jet  fiiel  present  when  modeling  the 
WTC  fires. 

•  The  presence  of  a  lower  volume  of jet  fuel  would  be  likely  to  have  had  a  less  than  linear  effect 
on  the  peak  HRR  value.  In  Test  5,  the  addition  of  4  L  ( 1 .06  gal)  of  Jet  A  uniformly  to  the 
surfaces  of  the  generic  workstation  resulted  in  a  6  percent  increase  in  the  peak  HRR  over  the 
value  from  Test  2.  The  addition  of,  for  example,  half  that  amount  would  have  resulted  in  an 
increase  of  less  than  3  percent  in  the  peak  HRR  value  because  some  amount  of  fuel  would 
have  been  consumed  before  the  peak  intensity  was  reached  and  that  amount  would  make  no 
contribution  to  the  peak  itself.  Since  the  effect  of  4  L  of  Jet  A  on  the  peak  HRR  was  small, 
the  HRR  results  would  not  have  been  likely  to  have  been  sensitive  to  lower  volumes  of  added 
jet  fuel. 

•  The  presence  of  a  larger  volume  of jet  fuel  would  be  likely  to  have  a  greater  than  linear  effect 
on  the  peak  HRR  value.  This  is  because  the  quantity  of  fuel  consumed  during  the  build-up  to 
the  peak  HRR  would  become  a  smaller  fraction  of  the  total  volume  added  and  the  volume  of 
fuel  remaining  for  pool-like  burning  at  the  peak  would  be  disproportionately  higher.  An 
upper  bound  to  the  sensitivity,  in  the  context  of  a  doubled  amount  of  Jet  A,  was  estimated  as 
follows: 

-    The  overall  HRR  was  approximated  as  the  sum  of  the  jet  fuel  "pool"  HRR  and  the  HRR 
from  the  burning  of  the  areas  not  covered  with  jet  fuel  (carpet  under  work  surface, 
underside  of  work  surface,  vertical  section  of  chair,  and  bookshelf).  It  was  assumed  that 
the  wetted  and  unwetted  surfaces  were  equal. 
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-  The  steady  state  HRR  of  a  deep  Jet  A  pool  with  area  equal  to  the  workstation  surface  area 
is  approximately  13  MW.  Here,  the  effective  layer  thickness  of  the  "poof  would  be 
approximately  1 .4  mm  and  for  such  a  thin  layer  the  actual  burning  would  be  transient 
with  a  significantly  reduced  HRR  compared  to  the  steady  state  value.  Gottuk  and  White 
(2002)  report  that  such  thin  fuel  layers  yield  only  about  1/5  the  steady  state  HRR  for 
several  types  of  surfaces.  By  analogy,  the  thin  Jet  A  fuel  on  the  desk  surfaces 
(approximately  62  percent  of  the  wetted  area)  would  have  a  steady  state  HRR  of  1.6  MW. 
The  remainder  of  the  "pool"  would  be  soaked  into  the  carpet  in  the  open  center  of  the 
work  station.  Ma  et  al.  (2004)  found  that  at  an  approximately  equivalent  loading  level  of 
hydrocarbon  fuels  on  nylon  caipeting,  the  HRR  is  essentially  unchanged  from  that  of  a 
deep  pool.  Thus,  this  area  would  contribute  about  4.9  MW  to  the  total  HRR. 

-  The  HRR  for  full  involvement  of  the  "dry"  surfaces  was  estimated  as  half  of  the  peak 
HRR  from  Test  2,  or  3  MW.  '  - 

-  The  sum  of  these  three  contributions  is  approximately  9.5  MW,  about  30  percent  higher 
than  the  peak  HRR  in  Test  5.  For  still  more  Jet  A,  the  layer  thickness  on  the  desk 
surfaces  will  not  get  appreciably  deeper  nor  will  the  HRR  from  such  wetted  surfaces. 
The  results  of  Ma  et  al.  (2004)  do  imply  that  the  HRR  from  the  central  carpet  area  could 
go  even  higher  with  higher  Jet  A  loadings.  This  could  have  occurred,  presumably,  very 
close  to  the  airplane  impact  area. 

•  The  presence  of  jet  fuel  had  only  a  small  effect  on  the  total  mass  loss  and  total  heat  released 
from  the  burning  workstation.  In  Test  5  (generic  workstation.  Jet  A)  the  mass  loss  of  the 
workstation  at  1 500  s,  the  time  at  which  the  ignition  burner  was  turned  off,  was  only 
approximately  10  percent  higher  than  in  Test  2  (generic  workstation,  no  jet  fuel).  The  total 
mass  loss  in  the  two  tests  was  indistinguishable.  The  total  heat  released  in  Test  2  was  slightly 
higher  than  in  Test  5.  As  discussed  below,  the  approximately  600  s  longer  exposure  to  the 
igniter  flame  and  the  800  s  longer  test  were  conducive  to  the  burning  of  more  of  the  paper.  It 
does  not  appear  that  varying  the  added  volume  of  jet  fuel  would  lead  to  a  significant  change 
in  this  small  effect. 

•  Increasing  the  quantity  of  jet  fuel  would  not  affect  further  the  time  to  full  involvement  of  the 
cubicle.  The  4  L  of  Jet  A  sharply  decreased  the  time  to  full  involvement  of  the  workstation 
from  530  s  to  less  than  160  s.  As  discussed  in  Sec.  3.2.3,  the  presence  of  Jet  A  and  the 
evaporative  cooling  effect  it  provided  actually  delayed  the  paper  auto-ignition  mechanism 
from  40  s  to  70  s  seen  in  "dry"  tests  to  1 50  s.  While  more  Jet  A  could  have  delayed  this 
mechanism  longer,  random  piloted  ignition  by  flying,  burning  debris  would  provide  an  upper 
limit  to  this  delay. 


3.4.3        Sensitivity  to  Fractional  Coverage  of  Surfaces  by  Inert  Material 

The  peak  HRR  was  decreased  roughly  in  proportion  to  the  20  percent  coverage  of  the  accessible 
combustible  surface  by  inert  tiles.  An  incremental  increase  in  coverage  would  be  expected  to  reduce  the 
peak  HRR  further  since  less  of  the  surface  area  would  have  been  accessible  for  burning.  Far  larger 
coverage  by  fallen  ceiling  tiles  was  not  likely  to  have  occurred.  Reduced-scale  tile-dropping  experiments 
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pointed  to  a  tendency  for  ceiling  tiles  to  accumulate  preferentially  on  the  open  central  area  of  the  carpet 
on  the  cubicle  floor,  covering  roughly  two  thirds  of  this  area.  Only  about  one-third  of  the  desk  surface 
became  covered  (Ohlemiller  2004).  Vertical  and  underside  combustible  surfaces  of  the  workstations 
remained  uncovered. 

Comparison  of  the  total  heat  released  in  the  tests  showed  no  effect  of  this  degree  of  coverage  by  inert 
material.  The  fire,  aided  by  thermal  radiation  from  the  hot  ceiling  layer,  was  able  to  sustain  the  burning 
albeit  at  a  slower  rate. 

Blanketing  of  the  cubicles  by  small  fragments  of  fallen  ceiling  tiles  and  shattered  wallboard  panels  could 
have  precluded  the  progressive  flame  spread  that  led  to  flame  involvement  of  all  accessible  surfaces. 
However,  there  were  clearly  locations  in  the  towers  where  intense  fires  existed  and  grew,  and  these  were 
the  most  likely  regions  in  which  the  thermal  environment  would  have  had  the  best  opportunity  to  weaken 
the  local  structural  components. 

3.4.4  Sensitivity  to  Mass  of  Paper  Products 

The  presence  of  paper  on  the  cubicle  surfaces  was  instrumental  in  the  ignition  of  the  heavier  workstation 
components.  Variation  in  the  mass  of  paper  is  not  likely  to  have  had  a  substantial  effect  on  this  role. 
However,  the  total  absence  of  exposed  paper  would  have  delayed  the  workstation  ignition  until  the  next 
most  susceptible  component  autoignited  or  random  flaming  debris  from  adjacent  fire  areas  served  as  a 
pilot. 

By  contrast,  incorrect  estimation  of  the  mass  of  paper  would  affect  the  total  heat  released.  Filling  the 
filing  cabinets  in  the  cubicle  could  have  added  150  kg  of  paper  and  increased  the  total  heat  released  from 
3.9  GJ  to  4.8  GJ  (based  on  an  upper  bound  estimate  of  30  percent  of  the  mass  of  this  additional  paper 
being  consumed).  There  would  not  likely  be  a  significant  increase  in  the  peak  HRR  due  to  the  poor  air 
supply  to  the  paper  in  the  file  cabinets.  However,  the  additional  heat  released  over  a  longer  period  of  time 
could  have  provided  a  longer  exposure  at  high  temperature  for  the  structural  components  of  the  building. 

3.4.5  Sensitivity  to  Workstation  Type 

The  overall  burning  behavior  of  the  two  types  of  workstations  was  quite  similar.  This  is  reasonable, 
given  the  general  similarity  in  layout,  combustible  mass,  and  composition  (Tables  2-1  through  2-4).  The 
following  two  factors  in  the  workstafion  composition  both  explain  the  modest  differences  in  burning 
behavior  and  could  lead  to  differences  from  or  similarities  to  other  workstation  types. 

•    The  net  heat  peak  HRR  for  the  WTC  workstation  was  about  10  percent  larger  than  that  for  the 
generic  workstation  (7.46  MW  versus  6.73  MW,  Table  3-1).  This  is  thought  to  result  fi-om 
the  flammable,  charring  plastic  surface  on  the  drawer  fronts  that  increased  the  exposed 
flammable  area  by  about  10  percent.  It  was  apparent  from  the  video  that  the  plastic  surface 
of  the  WTC  workstation  filing  cabinets  was  burning  during  the  test.  The  larger  mass  of  the 
chair  for  the  WTC  workstation  may  have  also  contributed  to  the  increased  HRR.  For  a 
typical  HRR  of  30  MJ/kg  for  these  particular  plastic  materials,  the  burning  of  1  kg  over  30  s 
would  be  equivalent  to  an  addifional  1  MW  of  heat.  So,  a  relatively  small  additional  mass  of 
these  materials  could  have  had  a  significant  effect  on  the  HRR. 
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•    The  WTC  workstation  had  a  file  cabinet  design  that  allowed  more  air  access  to  the  top  area  of 
filled  file  drawers  after  the  drawer  facing  burned  off  This  could  possibly  have  supported 
more  extensive  burning  of  paper  in  a  filled  file  cabinet  than  would  be  allowed  by  the  generic 
file  cabinets.  As  discussed  in  Sec.  4.2,  a  higher  mass  fraction  of  file  drawer  paper 
consumption  occurred  with  the  WTC  workstation. 

3.5  UNCERTAINTY  ANALYSIS 

3.5.1        Total  Heat  Release  Rate 

The  HRR  expressed  in  equation  (3-1)  is  a  function  of  gas  concentrations,  the  duct  flow  rate,  and  a 
number  of  constant  parameters  related  to  the  combustion  properties  of  the  fuel.  Bryant  et  al.  (2003) 
expressed  all  of  the  gas  analyzer  outputs  and  sensor  outputs  in  terms  of  their  measured  voltages  to 
eliminate  correlations.  One  example  of  correlated  measurements  is  that  of  the  ambient  oxygen 
measurement  (made  before  a  test)  and  the  oxygen  measurement  during  the  test,  both  of  which  are  made 
by  the  same  instalment.  The  following  equation  indicates  that  a  single  concentration  measurement 
requires  5  voltages  and  one  span  gas  concentration: 


V  -V 

^  O,   ~  y  y  02.Span  ^O..Zero) 

''^02. Span  ~     On, Zero 


Ot  .press , Zero 

V 

Y        O  2.  press  J 


(3-5) 


where: 


=  voltage  measured  by  oxygen  sensor  at  time,  t 
Zero  ^  voltage  measured  by  oxygen  sensor  for  the  zero  gas 
Vq^  5;,uh  ^  voltage  measured  by  oxygen  sensor  for  the  span  gas 

5^,^,,  =  oxygen  volume  fraction  of  the  span  gas,  bottled  air  (0.2095) 
X Zero  ~  oxygen  volume  fraction  of  the  zero  gas,  nitrogen 

X =  oxygen  volume  fraction  computed  with  respect  to  voltage  measurement  at  time,  t 

p^.^^^  =  voltage  measured  by  absolute  manometer  to  monitor  pressure  in  the  oxygen  analyzer  at 
time,  t  , 

Zero  ^  voltage  measured  by  absolute  manometer  to  monitor  pressure  in  the  oxygen 
analyzer  during  the  zeroing  of  the  instrument. 

The  modified  version  of  equation  3-1  included  78  input  voltages/parameters.  About  28  of  these  were  gas 
analyzer  voltages/calibration  constants,  and  30  were  related  to  measurement  of  the  velocity  profile  at 
eight  locations  using  pressure  transducers  to  record  the  pressure  difference  in  the  bidirectional  probes. 
The  remainder  mainly  involved  the  various  constants  appearing  in  the  HRR  calculation.  The  equation  for 
q  is  written  symbolically  as: 
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q  =  q{x^,x,,x^,  ,Xj^)  (3-6) 

By  using  voltages,  the  input  quantities  are  uncorrclated,  and  the  relative  total  type  B  standard  uncertainty, 
i/fi,  is  given  by: 


=-^Z^' ("(^;))'  (3-7) 

In  equation  (3-7),  ufxj  is  the  relative  standard  uncertainty  for  each  input,  and  Si  is  the  associated 
dimensionless  sensitivity  coefficient  given  by: 

dy  1 

^1=-   (3-8) 

dx.  y 

From  the  complete  analysis  based  on  the  uncertainties  for  the  78  quantities,  the  resulting  value  of  ub  is 
0.061. 

Section  3.3  reported  that,  based  on  fast  response  measurements  with  a  radiometer,  the  calorimeter 
underestimated  the  instantaneous  peak  heat  release  rate  by  about  3  percent.  Per  the  discussion  in  that 
section,  the  radiometer  output  was  not  directly  proportional  to  heat  release,  so  this  correction  term  itself 
had  a  large  uncertainty.  For  this  reason,  this  effect  was  treated  as  a  Type  B  uncertainty  with  i/(pk)  =  ± 
0.03  rather  than  correcting  the  readings  by  3  percent.  Including  this  term  in  the  total  Type  B  uncertainr\' 
resuhed  in  0.068. 

A  simplified  approximate  analysis  was  made  for  these  large  heat  release  rates,  because  the  uncertainty 
associated  with  the  gas  measurements  was  minor  compared  to  the  three  major  sources.  The  dominant 
relative  standard  uncertainties  were  the  bi-directional  velocity  probe  constant,  u(C).  the  HRR  per  kg  of 
oxygen  consumed,  i/(EHC),  and  the  cross  sectional  area  of  the  duct,  ii(A);  the  corresponding  values  were 
0.050,  0.027,  and  0.010.  The  HRR  was  approximately  proportional  to  the  product  of  these  three 
quantities: 

q«(EHC){C)(A)  (3-9) 

The  product  (C)(A)  is  part  of  the  expression  for  the  mass  flow  in  the  exhaust,  rh^ .  There  is  a  correction 
term  to  the  product  (EHC)^  in  equation  (3-9)  with  a  value  typically  less  than  1  percent  of  (EHC)^zJ  for 
the  conditions  of  these  workstation  bums.  The  value  of  ub  was  computed  from  equation  (3-6)  with 
sensitivity  factors  equal  to  1 . 

=  [{uiC)y-  +u{A)f  +uiEHC)f  \  "=0.058  (3-10) 

Adding  in  quadrature  the  value  for  w(pk)  gave  ub  =  0.065.  This  simplified  analysis  accounted  for  most  of 
the  type  B  uncertainty. 

The  Type  A  relative  standard  uncertainty  for  the  calorimeter,  Ua  (Cal.),  was  obtained  by  perfonning  repeat 
measurements  with  the  calorimeter  of  a  validation  gas  burner.  The  measured  value  was  0.017.  This 
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repeatability  included  both  that  of  the  calorimeter  and  that  of  the  burner.  The  relative  combined 
uncertainty  was  obtained  by  taking  the  root-sum-of-squares  for  the  Type  A  and  Type  B  uncertainties  with 
a  resuh  of  0.063  for  the  smoothed  peak  and  a  value  of  0.070  for  the  instantaneous  peak. 

There  were  two  issues  regarding  fire  measurement  repeatability.  One  was  the  repeatability  of  the 
calorimeter.  This  quantity  was  computed  above.  The  second  was  the  repeatability  of  the  workstation 
bums,  for  which  there  were  no  repeat  measurements.  This  precluded  the  ideal  direct  assessment  of 
repeatability.  However,  there  was  an  indirect  method  of  estimating  the  repeatability  that  entailed  a 
plausible  assumption  and  some  equivalence  calculations. 

The  rationale  was  as  follows:  The  repeatability  quantity  of  interest  to  estimate  was  Grepeat^  the 
(population)  standard  deviation— a  value  which  reflects  the  natural  spread  of  repeated  measurements.  If 
there  were  n  such  repeated  measurements,  it  would  be  natural  to  estimate  arepeat  by 


<7. 


epeal 


(  V 

V  /  ) 


-lI/2 


(3-11) 


In  the  absence  of  repeats,  however,  consider  the  following:  The  data  analysis  in  Table  3-2  shows  that  an 
estimated  interaction  effect  Pi  (for  any  of  the  responses)  was  quite  small,  of  the  order  of  .02  at  its  worst. 
If  it  is  assumed  that  the  population  interaction  effect  fii  was  in  fact  zero,  one  could  then  pose  the  question 
as  to  why  the  sample  estimate  fii  of  the  interaction  effect  was  not  also  identically  zero.  The  answer  to  this 
question  was  that  random  variation  caused  the  estimates  of  Pi  to  differ  from  effect  "truth"  fij 


Such  variation  in  fii  was  characterized  and  related  to  CFn-ped  as  follows: 

.2 


cj'p^  =  Var(p,  )  =  Var(y,  -  v  J  =  Var(y,  )  +  Var(y,  ) 
=  Var(y)/2  +  Var(y)/2  =  Var(y)  =  cr^^^, 


(3-12) 


Thus,  at  the  population  level  with  this  zero  interaction  assumption,  Crepeat  was  identically  ^/],  .  Further 
discussion  of  this  approach  can  be  found  in  Box  et  al.  (1978). 

Since  with  two  factors,  there  is  only  a  single  two-term  interaction,  the  estimate  of  cr^^  becomes  simply 


l(A-o)  /I 


1/2 


|A  1-0.02 


(3-13) 


Thus,  the  estimate  of  a, 


repeat 


0.02. 


In  summary,  this  estimate  of  repeatability  was  an  indirect  surrogate  of  the  preferred  replicate-based 
estimate;  on  the  other  hand,  in  the  absence  of  replication,  this  served  as  a  rigorous  estimate  of 
repeatability,  contingent  on  the  assumption  of  the  population  interaction  being  identically  zero,  and  the 
existence  of  only  a  single  interaction  on  which  to  base  our  estimate. 

The  relative  combined  uncertainty  of  the  calorimeter  for  the  workstation  tests  was  computed  as  the 
quadrature  sum  of  the  relative  combined  uncertainty  of  the  calorimeter,  0.063,  and  the  relative  type  A 
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uncertainty,  which  is  the  same  as  o^epeut,  0.02,  with  a  resuh  of  0.066.  This  was  the  result  for  the  smoothed 
peak.  For  the  instantaneous  peak,  the  value  increased  to  0.073. 

Finally,  the  relative  expanded  uncertainty  was  determined  from  the  combined  uncertainty  by  multiplying 
by  a  coverage  factor  of  2  giving  Ur  =  ±  0. 1 32  for  the  smoothed  peak  and  U  =  ±  0. 1 46  for  the 
instantaneous  peak. 

3.5.2  Net  Heat  Release  Rate 

Following  the  same  calculation  path,  the  relative  repeatability  standard  uncertainty  for  the  net  HRR  is 
0.02.  The  expanded  uncertainty,  U,  is  slightly  increased  to  ±0.136  for  the  smoothed  peak  and  ±  0. 1 50  for 
the  instantaneous  peak  because  the  difference  of  total  HRR  and  the  burner  release  rate  is  computed. 

3.5.3  Net  Heat  Released 

Two  additional  Type  B  uncertainties  for  the  total  HRR  are  the  uncertainty  of  the  spray  burner  (igniter) 
output  and  the  uncertainty  in  the  drift  of  the  baseline  resulting  from  a  drift  in  the  oxygen  meter.  The 
burner  Type  B  relative  standard  uncertainty,  0.035,  is  estimated  as  1/3  of  the  variation  of  the  burner 
output  during  its  operation  times  the  time  operating  divided  by  the  total  heat  released.  The  baseline  drift 
Type  B  standard  uncertainty,  0.017,  is  estimated  as  the  rate  of  drift  over  the  last  400  s  of  operation  times 
the  length  of  the  bum  after  the  liquid  burner  is  extinguished  divided  by  the  total  HRR.  The  total  Type  B 
relative  uncertainty  is  computed  as  a  quadrature  sum  of  the  spray  burner,  baseline  drift,  and  the  Type  B 
contribution  from  the  calorimeter  (0.061 )  with  a  result  of  0.072.  The  calorimeter  is  able  to  integrate  peaks 
too  narrow  to  be  quantified  so  that  the  uncertainty  associated  with  the  peak  height  is  not  included  in  the 
analysis.  The  relative  expanded  uncertainty  of  0.149  is  obtained  from  the  quadrature  sum  of  the  total 
Type  B  relative  uncertainty  and  the  relative  repeatability  uncertainty,  0.02,  multiplied  by  the  coverage 
factor  of  two. 

3.5.4  Total  Mass  Loss 

A  calibration  weight  of  22.7  kg  (50  lb)  was  added  to  the  load  cell  before  each  test,  and  the  relative 
standard  uncertainty  for  calibration,  Ucai.  is  estimated  as  0.005.  This  result  included  the  effects  of  both 
repeatability  and  the  calibration  uncertainty  of  the  22.7  kg  calibration  mass.  The  finest  increment  in  the 
output  for  a  single  load  cell  was  0.045  kg.  For  the  combined  system  of  four  load  cells,  it  was  obser\'ed 
that  the  background  reading  changed  in  discrete  increments  of  about  0.1  kg.  This  limited  resolution  was  a 
major  contributor  to  the  value  of  Ucai.  An  additional  source  of  uncertainty  during  the  fire  tests  was  objects 
falling  off  the  load  cell.  These  were  discrete  events  as  is  apparent  in  Fig.  3-9  at  1500  s  and  2000  s.  There 
was  an  uncertainty  in  measuring  this  step  change  in  mass  that  resulted  in  an  estimated  relative  standard 
uncertainty  of  approximately  0.030  for  the  total  mass  change,  which  was  typically  about  40  kg  relative  to 
a  total  mass  loss  of  approximately  200  kg.  The  relative  standard  uncertainty  for  falling  objects,  Ufau, 
expressed  as  a  fraction  of  the  total  weight  was  0.006.  The  relative  combined  uncertainty  was  obtained  by 
adding  Ucai  and  the  relative  repeatability  uncertainty,  0.02  in  quadrature  with  a  resuh  of  0.022.  The 
relative  expanded  uncertainty,  U,  with  a  coverage  factor  of  2  was  0.044. 

There  is  another  issue  in  terms  of  repeatability  related  to  the  test  protocol.  The  time  the  ignition  burner 
operated  and  the  total  time  of  the  test  were  shorter  for  the  tests  with  Jet  A  because  the  time  to  the  peak 
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HRR  was  about  350  s  shorter.  This  confounded  the  factor  analysis,  since,  in  effect,  two  additional 
parameters,  the  operation  time  of  the  ignition  burner  and  the  bum  time  of  the  workstation,  were 
inadvertently  changed.  In  Table  3-3,  the  mass  loss  values  for  the  four  tests  of  the  generic  workstation  are 
given  for  four  time  intervals.  The  shortest,  1500  s,  corresponds  to  the  one  case  where  all  four  tests  have 
the  same  environment — igniter  burner  on  for  that  entire  time,  hi  this  case,  the  mass  loss  is  larger  with  the 
presence  of  Jet  A.  On  the  other  hand,  if  one  compares  the  results  at  the  end  of  the  bum,  then  there  is  more 
mass  loss  without  Jet  A.  This  is  likely  a  result  of  the  longer  bum  time  and  longer  exposure  to  the  ignition 
flame,  causing  additional  mass  loss  from  the  paper  and  wood.  Also  there  is  a  reversal  in  the  effect  of  the 
tile  addition  from  reducing  the  total  heat  release  over  1 500  s  and  increasing  it  by  the  end  of  the  bum.  The 
results  for  2100  s  and  2500  s  are  intermediate  between  the  shortest  and  longest  times. 


Table  3-3.  Mass  loss  (kg)  vs.  time  interval. 


Time  Interval 

Test 

0  to  1500 

0  to  2100 

0  to  2500 

0  to  3200 

2 

163.5 

186.8 

195.7 

205.0 

3 

160.3 

188.2 

197.7 

213.1 

5 

181.3 

193.0 

200.2 

200.2 

6 

172.2 

196.3 

205.3 

205.3 

3.5.5  Peak  Mass  Loss  Rate 

The  MLR  ranged  from  a  peak  value  of  0.25  to  0.4  kg/s  and  then  decreased  to  a  value  as  small  as  0.01  kg/s 
at  long  times.  The  relative  uncertainty  in  the  peak  for  1  s  resolution  was  on  the  order  of  0.025  or  larger 
because  of  the  typical  0.1  kg  discretization  limit  of  the  load  cell.  The  MLR  was  smoothed  by  computing 
a  linear  fit  through  30  points.  This  reduced  this  Type  B  relative  uncertainty  by  roughly  the  square  root  of 
the  number  of  points  to  a  value  of  0.046.  The  rate  of  rise  of  the  MLR  near  its  peak  was  less  steep  than 
that  of  the  radiometer  output  and  calorimeter  (Fig.  3-16),  and  the  uncertainty  from  this  effect  was  treated 
in  a  similar  fashion  to  that  of  uncertainty  for  the  peak  HRR.  This  second  Type  B  uncertainty  was 
estimated  as  0.04  compared  to  the  value  of  0.03  for  the  calorimeter.  The  combined  relative  uncertainty 
was  obtained  by  combining  the  two  type  B  uncertainties  plus  the  0.02  repeatability  in  quadrature,  with  a 
resulting  value  of  0.064.  The  relative  expanded  uncertainty  was  obtained  as  twice  this  value,  0.128. 

3.5.6  Effective  Heat  of  Combustion 

The  uncertainty  in  the  effective  heat  of  combustion  was  computed  as  the  quadrature  sum  of  the  relative 
combined  uncertainty  of  the  net  heat  released  and  the  mass  loss.  The  resulting  value  was  multiplied  by  a 
coverage  factor  of  two  resulting  in  an  expanded  relative  uncertainty  equal  to  ±  0.150. 

3.5.7  Time  to  Peak  Heat  Release  Rate  and  Peak  Mass  Loss  Rate 

The  relative  standard  uncertainty  for  repeatability  was  0.020,  as  it  is  for  the  other  quantities.  The  major 
Type  B  component  for  the  peak  heat  release  rate  was  the  estimated  5  s  standard  uncertainty  associated 
with  correcting  the  nominal  30  s  time  delay  in  the  calorimeter  for  the  time  for  the  gas  to  flow  from  the 
duct  to  the  gas  analyzers  and  the  time  response  of  the  analyzers.  The  5  s  corresponded  to  a  relative 
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standard  uncertainty  of  0.010  for  the  peak  times  for  tests  not  involving  Jet  A  and  a  value  of  0.03  tests 
involving  Jet  A.  Combining  the  two  uncertainties  in  quadrature  and  multiplying  by  a  coverage  factor  of 
2  gave  relative  expanded  uncertainties  for  the  time  to  peak  of  0.022  for  the  tests  without  Jet  A  and  0.036 
for  the  tests  with  Jet  A. 

3.5.8        Burning  Time  and  Ignition  Time 

The  major  source  of  uncertainty  in  the  FWHH  for  the  net  HRR  was  determining  ti/2+,  the  time  after  the 
peak  at  which  the  net  heat  release  has  decreased  by  a  factor  of  two.  The  net  HRR  had  oscillatory 
behavior  in  this  region,  and  there  were  two  or  more  times  corresponding  to  the  half-height  of  the  net  heat 
release  weight.  The  average  of  the  shortest  and  longest  times  was  used  in  computing  a  smoothed  estimate 
of  t|/2+  and  is  denoted  t]/2+  (avg).  The  actual  value  of  ti/2+  was  within  the  interval  defined  by  the  shortest 
and  longest  times.  Assuming  the  values  in  this  interval  were  equally  probable,  the  standard  deviation 
about  t|/2+  was  equal  to  half  the  interval  divided  by  V3.  This  value  was  divided  by  the  FWHH  and 
multiplied  by  a  coverage  factor  of  two  to  obtain  the  expanded  uncertainty.  The  much  larger  uncertainty 
for  Test  3  compared  to  the  others  was  a  result  of  the  nearly  flat  net  heat  release  versus  time  for  an 
extended  period. 

The  uncertainty  in  the  total  mass  loss  was  one  of  the  two  major  sources  of  uncertainty  in  the  time, 
t(75  %),  at  which  75  percent  of  the  total  mass  loss  occurs.  This  uncertainty  was  estimated  as  one-half  the 
time  interval  determined  by  75  percent  of  the  total  mass  loss  plus  or  minus  the  mass  loss  combined 
uncertainty.  For  example,  for  a  total  mass  loss  of  205  kg,  the  time  interval  would  be  from  the  time  at 
which  the  mass  was  152.5  kg  to  the  time  it  was  155.0  kg.  For  Tests  2  and  3,  the  average  relative  standard 
uncertainty  was  0.01 1,  and  for  Tests  5  and  6  it  was  0.019.  Combining  these  values  in  quadrature  with  the 
relative  standard  uncertainty  for  repeatabihty  of  0.02  and  multiplying  by  a  coverage  factor  of  2  leads  to  an 
expanded  uncertainty  for  t(75  %)  equal  to  0.046  for  Tests  2  and  3  and  0.055  for  Tests  5  and  6. 

One  source  of  uncertainty  in  t,„  was  estimating  the  ignition  time  from  the  video.  This  standard  uncertainty 
was  estimated  as  ±  1  s.  The  standard  uncertainty  associated  with  repeatability  was  computed  as  the 
standard  deviation  of  the  four  tests  without  Jet  A  (Tests  1  through  4)  and  was  equal  to  ±  12  s.  The 
expanded  uncertainty  based  on  a  coverage  factor  of  2  was  ±  24  s.  The  uncertainty  is  expressed  in  terms 
of  relative  expanded  uncertainty  in  Table  3^. 
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Table  3-4.  Uncertainty  values  for  test  parameters. 


Parameter 

Uncertainty  (k=2) 

Peak  heat  release  rate"" 

14.6  %/13.2  % 

Time  lo  peak  for  both  HRR  and  MLR  ^ 

4.4  %/7.2  % 

Net  peak  HRR' 

15.0%/13.6% 

Peak  MLR 

12.8  % 

Net  heat  released  by  the  workstation 

14.9  % 

Total  mass  loss 

4.4  % 

Effective  heat  of  combustion 

15.0% 

FWHH  for  net  HRR^ 

11. 6  %,  43.6  %,  10.8  %,  12.8  % 

t(75  %)' 

4.6  %,  5.5  % 

62  %,  36  %,  43  %,  48  %,  27  %,  21  % 

a.  The  first  number  is  the  uncertainty  for  the  instantaneous  peak  and  the  second  is  the  10  s  average  about  the  peak. 

b.  The  first  number  is  the  uncertainty  without  Jet  A  and  the  second  number  is  with  Jet  A. 

c.  The  first  number  is  the  uncertainty  for  the  instantaneous  peak  and  the  second  is  the  1 0  s  average  about  the  peak. 

d.  The  four  values  are  the  uncertainties  for  Tests  2,  3,  5,  and  6. 

e.  The  first  value  is  the  uncertainty  for  Tests  2  and  3.  and  the  second  is  for  Tests  5  and  6. 
f  The  six  values  of  uncertainties  are  for  Test  1  through  Test  6. 
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Simulations  of  Workstation  Fire  Tests 


4.1  FIRE  dynamics  SIMULATOR  COMBUSTION  MODELING 

The  purposes  of  these  experiments  were  to  obtain  understanding  of  how  the  office  workstations  may  have 
burned  in  a  fire  and  to  provide  data  for  use  any  improvements  in  the  Fire  Dynamics  Simulator  (FDS) 
combustion  algorithm  and  its  inputs  needed  to  provide  a  sufficient  approximation  of  the  burning  of 
combustibles  as  complex  as  those  that  comprise  an  office  workstation.  A  general  description  of  the 
implementation  is  presented  here.  A  more  detailed  version  is  presented  in  NIST  NCSTAR  1-5F. 

The  existent  version  of  FDS  computed  the  radiant  flux  to  the  surfaces  of  the  workstation.  Given  the 
thermophysical  properties  of  each  combustible  material,  including  the  ignition  temperature,  heat  of 
vaporization/pyrolysis,  and  the  thermal  conductivity,  FDS  computed  when  each  ignited  and  its 
subsequent  burning  rate.  The  effluent  was  added  to  the  thermal  and  radiative  environment.,  and  the 
computation  proceeded  with  successive  iterations  until  the  fuel  was  fully  consumed. 

Section  2.1  describes  in  detail  the  various  materials  making  up  the  workstations  and  Appendix  A 
describes  the  cone  calorimeter  measurements  of  the  properties  needed  for  the  FDS  calculations.  The 
measurements  to  characterize  the  burning  characteristics  were  made  at  three  fluxes  for  the  carpet,  desk 
(wood),  computer  monitor,  chair,  privacy  panel,  and  stacked  paper.  Additional  piloted  ignition  time 
measurements  at  other  heat  fluxes  were  made  on  these  six  materials  to  infer  their  ignitability  parameters. 

Direct  use  of  all  of  these  data  did  not  result  in  acceptable  simulation  of  the  workstation  fire  growth 
behavior.  One  significant  contributor  to  this  disagreement  was  that  only  piloted  ignition  had  been 
measured,  and  as  noted  in  Sec.  3.2.2,  the  paper  throughout  the  cubicle  appeared  to  have  been  ignited  from 
irradiance  alone  (i.e.,  it  exhibited  non-piloted  or  spontaneous  ignition).  In  addition,  the  chair  fire,  whose 
complex  behavior  was  a  key  determinant  of  the  timing  and,  to  a  lesser  extent,  the  magnitude  of  the  heat 
release  peak  from  the  work  station,  could  not  be  predicted  from  the  cone  calorimeter  results  alone. 

Thus,  only  the  carpet,  desk,  and  privacy  panel  data  from  the  cone  calorimeter  were  used  direcfly. 

•  The  carpet  and  privacy  panel  were  modeled  as  thermoplastics,  that  is,  the  burning  rate  was 
assumed  to  be  proportional  to  the  heat  flux  from  the  surrounding  environment. 

•  The  desk  was  modeled  as  a  charring  solid,  in  which  a  pyrolysis  front  propagated  through  the 
material,  leaving  a  layer  of  char  behind  that  insulated  the  material  and  reduced  the  burning 
rate.  Details  of  the  pyrolysis  models  can  be  found  in  the  FDS  Technical  Reference  Guide 
(McGrattan  2004).  More  details  on  the  specific  FDS  model  inputs  used  are  given  in 
Appendix  A. 

•  The  chair,  computer,  paper,  and  other  miscellaneous  items  within  the  workstation  were 
modeled  as  a  single  item  by  lumping  their  masses  together  into  large,  homogeneous  "boxes" 
and  distributing  these  throughout  the  workstation. 
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•  Paper  cartons  containing  various  amounts  of  plastic  have  been  used  as  surrogates  for  a  wide 
range  of  retail  commodities  in  other  fire  analyses  and  simulations.  A  particular  distribution 
of  plastic  in  a  carton  teiined  "Group  A  Plastic"  has  been  shown  to  be  fairly  representative  of 

;  fires  fueled  by  a  mixture  of  cellulosic  and  thermoplastic  materials.  Because  the  Group  A 

Plastic  had  been  used  in  numerous  FDS  simulations,  the  contents  of  the  office  workstations 
was  modeled  with  a  fuel  similar  to  Group  A  Plastic. 

•  Blind  predictions  of  the  single  open  workstation  bums  were  made  using  material  properties 
obtained  previously  during  a  sprinkler/roof  vent  study  of  large-scale  rack  storage  fires 
(McGrattan,  Hamins,  and  Forney  2000).  These  properties  were  then  adjusted  to  match  the 
results  of  the  current  experiments.  Thus,  the  results  from  the  single  workstation  bums  served 
to  calibrate  the  model. 

•  The  nature  of  the  boxes  that  provided  a  reasonable  simulation  of  the  cubicle  fires  was  as 
follows: 

-  Three  boxes,  each  0.5  m  (20  in.)  square  by  1 .0  m  (39  in.)  high. 

-  The  ignition  temperature  was  370  °C,  and  the  heat  release  rate  (HRR)  per  unit  area  was 
.  450kW/m-. 

-  The  HRR  of  each  box  underwent  a  100  s,  linear  ramp  from  the  point  of  ignition. 

-  The  heat  release  of  the  boxes  decreased  in  time  resulting  from  a  decrease  in  surface  area 
as  the  boxes  burned. 

-  Unlike  the  desk,  partition,  and  carpet,  the  boxes  were  simply  given  a  burning  rate  rather 
than  a  heat  of  vaporization,  meaning  that  the  boxes  would  bum  at  the  given  rate 
regardless  of  the  heat  flux  upon  them.  The  reason  is  that  the  complex  details  of  the 
ignition  and  fire  spread  on  a  chair,  followed  by  the  ignition  of  the  carpet  via  flaming  drips 
from  the  chair  were  not  sufficiently  understood  to  be  incorporated  into  FDS. 

The  obstmctions  (tiles)  were  modeled  by  placing  inert  objects  on  the  desk  and  floor,  covering  roughly 
30  percent  of  the  surface  area. 

The  addition  of  Jet  A  was  modeled  as  a  single  nozzle  spraying  liquid  fuel  droplets.  Four  liters  of  fuel  are 
released  in  2  s  and  combust  immediately. 

4.2  FDS  SIMULATIONS 

Figures  4-1  through  4-4  compare  the  predicted  HRR  over  the  duration  of  four  of  the  tests  with  the 
measured  equivalent.  Since  Test  1  was  used  only  for  scoping  the  experiments,  no  simulation  was 
performed.  The  FDS  simulation  for  Test  2  applies  implicitly  to  Test  4,  since  the  only  (small)  difference 
in  the  input  data  was  the  combustible  mass. 

The  workstation  bums  served  to  check  that  the  assumptions  made  in  modeling  the  workstation  with  the 
simplified  fuel  packages  would  produce  reasonable  results.  Shown  in  Fig.  4-1  is  a  comparison  of  the 
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HRR  from  Test  2  (no  jet  A,  no  tiles)  and  that  from  the  FDS  simulation.  The  initial  small  peak  near  t=0  in 
the  simulation  was  from  the  ignition  of  the  fabric  of  the  workstation  panels,  an  event  that  did  not  occur  in 
the  test.  The  peak  heat  release  occurred  at  full  involvement  of  everything  except  one  of  the  boxes  and 
perhaps  a  portion  of  the  underside  of  the  desk.  The  value  of  the  HRR  peak  is  slightly  larger  ( 1 0  percent ) 
than  the  experimental  value.  Also  the  peak  structure  includes  a  double  peak.  The  decrease  in  the  heat 
release  past  the  peak  is  a  resuh  of  the  decreasing  surface  area  of  the  boxes  and  of  the  formation  of  a  char 
layer  on  the  work  surface.  The  peak  HRR  in  the  simulation  occurred  sooner  than  in  the  experiment,  in 
the  experiment,  the  time  to  peak  HRR  was  strongly  influenced  by  the  melting  of  the  chair  plastic  onto  the 
carpet.  This  level  of  detail  was  not  captured  in  the  numerical  model,  especially  given  the  fact  that  the 
chair  had  been  lumped  together  with  various  other  combustible  items. 
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Figure  4-1.  Comparison  of  measured  and  computed  HRR  in  Test  2 
(generic  workstation,  no  tiles,  no  Jet  A). 
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Figure  4-3.  Comparison  of  measured  and  computed  HRR  in  Test  5 
(generic  workstation,  no  tiles,  with  Jet  A). 
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Figure  4-4.  Comparison  of  measured  and  computed  HRR  in  Test  6 
(generic  workstation,  with  tiles,  with  Jet  A). 

The  predicted  decrease  in  the  peak  heat  release  rate  resuhing  from  tiles  (Figs.  4-2  and  4-4)  and  the 
increase  in  the  peak  from  the  Jet  A  (Figs.  4-3  and  4-4)  are  similar  in  magnitude  to  the  experimentally 
measured  results.  The  simulations  did  predict  a  relatively  small  increase  in  the  time  to  the  peak  with  the 
addition  of  the  tile,  but  they  did  not  predict  the  more  than  a  factor  of  three  reduction  in  the  time  to  the 
peak  for  the  addition  of  Jet  A.  FDS  consumed  the  Jet  A  immediately  with  relatively  small  effect  on  the 
growth  of  the  heat  release  rate.  For  the  actual  fire,  the  ignition  of  the  Jet  A  did  not  occur  until  about 
120  s.  Then  the  ignition  and  rapid  burning  of  the  Jet  A  caused  the  entire  cubicle  quickly  to  become  fijlly 
involved,  first  above  the  desk  surface  and  then  below  via  dripping  from  the  chair  to  the  carpet. 

One  qualitative  difference  between  experiment  and  simulation  is  that  the  predicted  peak  structure  is 
similar  to  a  square  wave,  while  the  measured  HRR  decreases  in  a  more  nearly  linear  fashion  at  longer 
times.  The  overall  burning  time  is  similar  for  the  experiments  and  the  simulations. 

4.3  NUMBER  OF  WORKSTATIONS  BURNING 

To  put  into  context  the  importance  of  the  deviations  of  these  simulations  from  the  experimental  results,  it 
was  useful  to  have  an  estimate  of  the  number  of  workstations  involved  in  the  actual  World  Trade  Center 
(WTC)  fires. 

The  estimated  peak  HRR  of  the  fires  in  WTC  1  was  of  the  order  of  1  GW  to  1 .5  GW  (NIST 
NCSTAR  1-5F),  with  fires  being  observed  on  the  92nd  through  1 00th  floors  (NIST  NCSTAR  1-5 A).  The 
peak  HRR  of  the  WTC  workstation  identical  in  design  to  those  used  in  the  primary  impact  floors  of 
WTC  1  was  approximately  7  MW.  These  numbers  suggest  that  the  WTC  1  peak  HRR  was  the  equivalent 
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of  about  200  workstations  at  their  peak  HRR  values.  However,  it  was  likely  that  at  any  moment,  many  of 
these  workstations  were  not  generating  their  peak  heat  output  because: 

•  Not  all  would  have  been  at  the  same  point  in  their  combustion  history  simultaneously, 

•  Many  were  at  least  partially  covered  with  debris  from  fallen  ceiling  tiles  and  wall 
fragments,  and 

•  Not  all  would  have  free  access  to  an  unrestricted  supply  of  air. 

Thus,  a  greater  number  of  workstations  were  burning  as  the  overall  fire  peaked.  Assuming  that  the 
number  was  doubled  and  that  they  were  spread  evenly  over  the  nine  floors,  of  the  order  of 
40  workstations  would  have  been  at  some  stage  of  burning  on  any  given  floor.  There  were  on  the  order  of 
200  workstations  or  equivalents  on  each  floor  (NIST  NCSTAR  1-7),  so  about  one-fifth  of  the 
workstations  on  each  floor  would  have  been  burning  during  the  most  intense  stages  of  the  fire.  The  fire 
intensity  on  each  floor  (and  thus  the  number  of  work  stations  burning)  would  not  be  expected  to  be  equal 
but  the  argument  still  holds  that,  in  general,  numerous  work  stations  were  burning  on  each  floor  as  the 
fires  reached  peak  intensity. 

With  a  large  number  of  workstations  burning  simultaneously,  with  each  workstation  burning  for 
approximately  20  min,  and  with  the  uncertainties  in  the  burning  conditions  in  the  interior  of  the  buildings, 
it  is  unlikely  that  the  differences  in  the  predicted  HRR  curves  and  ignition  delay  times  would  have  a 
dominant  effect  on  the  prediction  of  the  overall  fire  growth  and  magnitude. 

Two  further  factors  help  alleviate  any  concern  about  the  timing  differences  between  the  experimental  and 
simulated  HRR  curves  seen  here.  The  next  stage  of  the  experimental  program  involved  the  burning  of 
three  workstations  within  a  large  enclosure  (NIST  NCSTAR  1-5E).  The  ignition  of  the  first  workstation 
was  by  the  same  means  as  here  (2  MW  spray  burner),  but  the  involvement  of  the  other  two  workstations 
was  the  result  of  a  fire  spread  process  after  the  first  station  ignited.  FDS,  using  the  approach  inferred 
from  the  tests  described  in  this  report,  successfully  described  the  overall  HRR  behavior  of  the  three 
stations.  For  the  three  stations  and  also  for  the  full  WTC  floor  simulations  (NIST  NCSTAR  1-5F),  the 
limiting  factor  in  determining  the  overall  HRR  evolution  over  time  was  found  to  be  the  rate  of  air  supply 
to  the  fire.  Therefore,  the  results  were  relatively  insensitive  to  the  types  of  timing  discrepancies  seen  in 
the  open  bums  here. 

Thus,  overall,  the  chosen  set  of  parameters  and  approximated  component  burning  descriptions  gave  an 
adequate  description  of  the  actual  workstation  heat  release  rate  behavior  and  its  dependence  on  the  factors 
varied  in  the  tests  described  here  (inert  tile  coverage  and  presence  of  Jet  A).  That  description  was  found 
to  be  robust  when  applied  to  the  larger,  air-supply  limited  fires  most  pertinent  to  the  WTC. 
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To  prepare  the  National  Institute  of  Standards  and  Technology  (NIST)  Fire  Dynamics  Simulator  (FDS) 
for  use  in  reconstructing  the  World  Trade  Center  (WTC)  fires,  it  was  necessary  to  establish  its  accuracy  in 
simulating  the  combustion  of  workstations,  the  principal  combustibles  in  the  WTC  buildings.  A  series  of 
six  fire  tests  of  individual  cubicles  was  conducted  to  provide  data  for  use  in  developing  an  FDS 
combustion  algorithm  capable  of  replicating  the  burning  of  such  complex  combustibles.  The  tests 
exammed  the  effect  of  var}'ing  the  t^'pe  of  workstation  and  inclusion  of  the  effects  of  inert  debris  from 
fallen  ceiling  tiles  or  fractured  walls  and  Jet  A  from  the  airplane.  The  principal  measurements  were  of  the 
heat  release  and  mass  loss  rates,  the  total  heat  release  and  mass  loss,  and  the  ignition  delay  time,  as  well 
as  the  uncertainty  in  each.  A  detailed  description  of  the  progress  of  each  fire  was  also  prepared.  The  test 
results  showed  the  following: 

•  Five  measures  of  fire  intensity  (heat  release  rate  [HRR].  mass  loss  rate,  carbon  dioxide 
volume  fraction  in  the  exhaust  gases,  temperature  of  the  exhaust  gases  and  radiant  flux  seen 
by  the  gage  viewing  the  entire  fire  from  the  south  side  of  the  cubicle)  all  presented  a 
consistent  representation  of  the  progress  of  the  fire,  given  their  individual  measurement 
characteristics. 

•  The  peak  HRR  for  a  workstation  similar  to  those  in  the  WTC  was  approximately  7.0  MW 
with  a  relative  expanded  uncertainty  of  ±  15  percent  and  with  the  T}^e  A  (repeatability) 
relative  standard  uncertainty  of  ±  2  percent.  The  time  to  the  peak  HRR  was  approximately 
540  s  with  an  expanded  relati\'e  uncertainty  of  4  percent. 

•  In  all  tests,  the  peak  HRR  corresponded  to  all  accessible  combustible  surfaces  burning 
simultaneously. 

•  Approximately  75  percent  of  the  heat  release  and  mass  loss  occurred  over  a  period  of 
approximately  20  min,  insensitive  to  all  the  test  variables. 

•  In  all  tests,  the  post-peak  behavior  included  both  the  burning  out  of  the  various  combustibles 
and  the  collapsing  of  three-dimensional  components  mto  more  compact  piles  that  decreased 
their  burning  rates. 

•  The  overall  burning  behavior  was  not  sensitive  to  workstation  design,  with  the  only  (small) 
effect  being  due  to  the  difference  in  mass  and  the  exposed  flammable  area. 

•  The  peak  HRR  of  the  workstation  was  diminished  in  rough  proportion  to  the  fraction  of 
combustible  surface  covered  by  inert  materials.  It  was  estimated  that  a  moderate  increase  in 
the  surface  co\'erage  would  ha\  e  decreased  the  peak  HRR  proportionately  and  that  a  very 
large  increase  in  coverage  (possible  in  heavily  rubbled  areas)  would  have  had  a  more  than 
proportionate  effect. 
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•  The  presence  of  4  L  of  Jet  A  dramatically  decreased  the  times  to  peak  HRR,  from 
approximately  10  min  to  approximately  3  min.  Increasing  the  quantity  of  jet  fuel  was  not 
expected  to  affect  further  the  time  to  full  involvement  of  the  cubicle.  The  presence  of  jet  fuel 

.  ,    ;       had  only  a  small  effect  on  the  total  mass  loss,  the  total  heat  released,  and  the  peak  HRR  from 
the  burning  workstation.  The  peak  HRR  value  was  estimated  to  be  insensitive  to  a  lesser 
amount  of  jet  fuel  present,  but  sensitive  to  a  larger  amount. 

•  The  tiles  and  Jet  A  were  nearly  independent  factors:  the  effect  of  the  tiles  on  the  burning  rate 
of  the  cubicle  was  almost  independent  of  the  presence  or  absence  of  the  Jet  A. 

•  Varying  the  mass  of  paper  in  the  cubicle  would  have  affected  the  duration  of  the  tailing  off  of 
the  workstation  fires.  .  > 

FDS  simulations  of  the  fire  tests  in  which  all  the  combustibles  were  volatilized  using  thermophysical 
properties  derived  using  a  cone  calorimeter  were  unsatisfactory.  This  was  due  to  fire  phenomena  that 
could  not  be  replicated  in  any  apparatus  that  used  a  small  test  specimen:  changes  in  the  shape  (e.g.,  by 
collapse)  and  surface  area  of  combustibles  during  the  fire  test,  three-dimensional  burning  behavior, 
especially  of  the  chair;  fabric  peeling  on  the  privacy  panels;  and  anticipation  of  whether  a  combustible 
would  ignite  from  direct  flame  impingement  or  by  radiation-driven  autoignition. 

In  a  modified  FDS  representation  of  the  workstation  combustion,  the  carpet,  desk,  and  privacy  panel  data 
from  the  cone  calorimeter  were  used  directly.  All  the  other  items  within  the  workstation  were  modeled  as 
a  single  entity  by  lumping  their  masses  together  into  three  homogeneous  "boxes"  and  distributing  these 
throughout  the  workstation.  The  combustion  properties  were  adjusted  to  match  the  results  of  the 
experiments.  Thus,  the  results  from  the  single  workstation  bums  served  to  calibrate  the  model. 

The  modified  FDS  then  computed  HRR  versus  time  curves  whose  shapes  and  magnitudes  were  similar  in 
form  and  magnitude  to  the  tests.  The  predicted  decrease  in  the  HRR  peak  and  time  to  that  peak  resulfing 
from  the  presence  of  inert  tiles  and  the  increase  in  the  peak  from  the  Jet  A  were  similar  in  magnitude  to 
the  measured  results.  There  were  differences  in  the  predicted  times  to  the  HRR  peaks.  It  was  anticipated 
that  these  effects  could  be  averaged  when  multiple  workstations  were  at  various  stages  of  their  burning. 
An  estimate  was  that  of  the  order  of  40  workstations  were  burning  on  each  floor  at  the  height  of  the  fires 
in  WTC  1 .  In  fact,  the  workstation  description  derived  here  was  successful  in  describing  the  overall  HRR 
from  the  spread  of  fire  over  three  workstations  in  a  single  compartment  (tests  conducted  in  another  part  of 
this  study).  There,  as  well  as  in  simulations  of  full  WTC  floor  burning,  the  factor  controlling  the  overall 
HRR  from  the  fire  was  the  rate  of  air  supply,  not  timing  of  the  HRR  peak  from  a  single  workstation. 

Thus,  the  chosen  set  of  parameters  and  approximated  component  burning  descriptions  gave  a  reasonable 
description  of  the  actual  workstation  heat  release  rate  behavior  for  application  to  compartment  fires 
involving  multiple  workstations. 
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Determination  of  Flammability  Properties  of  Office 
Workstation  Materials 


A.1       materials  examined 

The  Fire  Dynamics  Simulator  (FDS)  model  required  that  the  exposed  surface  of  each  significant 
flammable  object  be  described  by  an  appropriate  set  of  parameters  that  determine  its  case  of  ignition  and 
the  rate  at  which  it  evolves  heat  once  ignited.  For  these  workstations,  it  was  possible  to  break  the  most 
significant  surfaces  down  to  the  six  materials  or  composites  described  in  Table  A-1. 


Ta 

ble  A-1.  Descriptions  of  workstation  materials. 

Material 

Description 

Privacy  panel 

Fabric  cover  (220  g/m",  6.5  oz/yd",  unidentified  thermoplastic  fiber)  over  6  mm  of 
fiberglass  over  approximately  1  mm  perforated  steel  sheet  over  20  mm  layer  of 
fiberglass  with  an  aluminized  backing  paper.  Structure  repeated  past  this  mid-depth 
point  as  mirror  image.  Total  panel  thickness  was  54  mm."" 

Desk  surface 

Approximately  1  mm  thick  laminate  layer  glued  on  top  and  bottom  surfaces  of  a 
28  mm  thick  layer  of  medium  density  fiberboard  (mean  density  =  820  kg/m"). 

Carpet  tile 

Approximately  4  mm  deep  nylon  face  pile  bonded  to  3  mm  thick  backing  (lowest  layer 
was  approximately  1 .5  mm  thick,  dense  foam  rubber). 

Chair  seat 

Approximately  1  mm  thick  fabric  layer  (unidentified  charring  fiber)  over  40  mm  of 
polyurethane  foam  over  4  rmn  thick  plastic  shell  (probably  polyolefin)  (Chair 
assembly  also  included  2.3  kg  nylon  leg/foot  stand). 

Computer  monitor  shell 

ABS  shell  nominally  2.8  mm  thick.  17  in.  CRT  monitor  contained  little  other 
flammable  material. 

Paper 

75  g/m'  white  copier  paper  (30  percent  recycled  fiber).  Used  in  varying  layer 
thicknesses  in  full-scale  tests,  but  characterized  here  as  a  horizontal  layer 
approximately  25  mm  thick. 

a.  Note  the  minimal  amount  of  flammable  material.  Fiberglass  typically  contains  only  a  few  percent  of  an  organic  binder.  Each 
panel  did  have  a  wood  frame  around  its  periphery,  but  this  was  deemed  to  be  a  minimal  contributor  to  any  thing  but  the  tail 
end  of  the  burning  of  this  workstation. 


The  materials  used  for  the  characterization  tests  of  these  first  two  items  were  taken  from  the  generic 
workstation.  The  computer  had  a  separate  keyboard  and  processor  enclosure  made  of  differing  materials, 
but  these  minor  components  were  assumed  to  behave  the  same  as  the  monitor  shell. 

In  the  full-scale  tests  of  the  workstation,  the  paper  was  arrayed  in  both  horizontal  stacks  and  in  vertical 
groups  of  sheets  within  document  boxes.  However,  it  was  characterized  only  in  the  form  of  a  horizontal 
stack  (of  approximately  250  sheets). 
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A.2         TEST  PROTOCOL 

Square  specimens,  10  cm  (4  in.  )  square,  of  the  workstation  materials  were  cut  from  spare  parts.  A  spare 
chair  seat  was  similarly  cut  into  the  same  size  pieces,  being  careful  to  select  flat  portions.  The  computer 
monitor  shell  was  able  to  supply  sufficient  samples  only  by  restricting  the  size  to  75  cm  squares,  so  the 
heat  release  rate  (HRR)  per  unit  area  data  for  these  was  adjusted  to  compensate  for  this  smaller  size. 

All  samples  were  tested  in  a  cone  calorimeter  (ASTM  E  1354)  in  a  horizontal  orientation  for  both  ignition 
delay  time  and  HRR  as  a  function  of  incident  radiant  heat  flux.  An  electrical  spark  was  used  as  an  igniter 
for  the  evolved  gases  (piloted  ignition).  Triplicate  tests  were  conducted  at  heat  fluxes  of  20  kW/m", 
45  kW/m-,  and  70  kW/m".  In  the  case  of  the  desk  surface,  which  would  not  ignite  at  20  kW/m",  the  lower 
flux  tests  were  conducted  at  25  kW/m~.  The  results  of  the  exposures  at  these  three  fluxes  were  used  in 
fitting  a  linear  (least  squares)  variation  of  HRR  with  incident  flux.  Rather  than  use  the  peak  HRR  value, 
however,  the  line  was  fitted  to  points  derived  from  some  average  of  the  HRR  data  versus  time  at  a  given 
flux,  as  discussed  later  in  this  appendix.  The  ignition  data  at  these  fluxes  were  supplemented  by  farther 
duplicate  tests  in  which  the  flux  was  lowered  until  the  ignition  delay  time  exceeded  1 200  s  (taken  as  non- 
ignition).  The  borderline  between  ignition  and  non-ignition  was  determined  to  within  approximately 
±1.5  kW/m".  The  ignition  delay  time  versus  flux  was  then  used  to  infer  an  ignition  temperature  and 
effective  thermal  inertia  for  each  material  by  matching  the  experimental  curves  against  the  predictions  of 
a  one-dimensional,  numerical  heat  conduction  model. 

The  resulting  data  are  compiled  in  Table  A-2. 
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A.3         HEAT  RELEASE  RATE  BEHAVIOR 

The  measured  HRR  behavior  of  the  six  material  combinations  at  the  tested  fluxes  are  shown  in  Figs.  A-1 
tlirough  A-6.  None  of  the  materials  yielded  a  simple  constant  HRR. 

Neither  the  desk  surface  nor  the  paper  was  as  simple  as  pure  wood  in  physical  behavior.  Each  gave  a 
strong  initial  release  of  heat  followed  by  a  rapid  decay  as  the  char  layer  built  up  and  insulated  the  virgin 
material  below.  The  top  layer  of  material  on  the  desk  surface  (decorative  laminate)  separated  and 
fractured  irregularly,  forming  a  somewhat  erratic  radiation  shield.  The  horizontal  stack  of  copier  paper 
was  even  more  irregular  as  successive  sheets  curled  and  charred,  forming  a  thick  insulative  layer.  The 
bulk  of  the  subsequent  burning  was  on  the  sample  periphery. 

The  chair  seat  sample  behavior  was  complex  because  the  sample  itself  was  complex  (three  material 
layers).  The  initial  peak  was  due  to  the  fabric  burning,  which  left  a  char  layer  in  place  over  the  melting 
foam.  As  the  foam  melt  burned  out  it  ignited  the  more  massive  plastic  support  layer  beneath  it,  which 
burned  vigorously,  yielding  the  second,  larger  heat  release  rate  peak.^ 

The  carpet  sample  consisted  of  initially  distinct  layers  of  nylon  fibers  over  a  rubber  base.  The  nylon 
melted  and  then  burned,  igniting  the  lower  rubber  layer  in  a  continuum  that  yielded  only  a  single  HRR 
peak. 

The  ABS  polymer  of  the  computer  monitor  shell  consistently  gave  a  late  HRR  peak  that  is  normally 
associated  with  a  boiling  polymer  melt  layer.  No  such  melt  layer  was  evident  here,  and  the  overall 
behavior  was  suggestive  of  a  material  which  bums  more  vigorously  after  its  flame  retardant  content  has 
been  driven  off.  (Its  calculated  heat  of  combustion  went  up  strongly  as  the  heat  release  rate  peaked.) 


This  thermoplastic  layer  collapsed  toward  the  floor  when  the  full-scale  chair  was  burned,  causing  difficult-to-predict  changes 
in  HRR.  This  collapse  ignited  the  nylon  base  of  the  chair  (and  surrounding  carpet),  resulting  in  a  pool  fire  on  the  floor. 
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Figure  A-1.  Workstation  desk  surface  heat  release  rate  vs.  time. 
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Figure  A-2.  Workstation  chair  seat/back  heat  release  vs.  time. 
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Figure  A-3.  Workstation  carpet  heat  release  rate  vs.  time. 
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Figure  A-4.  Paper  stack  heat  release  rate  vs.  time. 
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Figure  A-5.  Computer  monitor  shell  heat  release  rate  vs.  time. 
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Figure  A-6.  Workstation  privacy  panel  heat  release  rate  vs.  time. 
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In  spite  of  these  complex  behaviors,  the  samples  were  treated,  for  modeling  purposes,  as  if  they  were 
exhibiting  a  simple  square-wave  heat  release  rate  behavior,  i.e.,  the  heat  release  rate  is  taken  to  jump  at 
ignition  to  some  constant  value  and  remain  at  that  value  until  the  material  bums  out.  This  was  a  usual 
approach  to  complex  materials  and  has  been  used  with  reasonable  success  in  describing  fire  growth 
(Dillon,  Kim,  and  Quintiere  1998).  It  was  not  feasible  to  model  the  complexities  of  most  real  material 
composites,  not  least  because  each  new  composite  exhibits  unique  behavior.'*  However,  by  replacing  a 
strong  initial  HRR  spike,  for  example,  with  a  flat  plateau  at  some  lower  (for  example,  overall  test 
average)  heat  release  rate,  one  impacted  any  prediction  for  how  flames  will  spread  on  a  material, 
particularly  in  an  upward  direction.  There  was  thus  an  incentive  to  develop  more  than  one  description  for 
some  of  the  materials  and  choose  between  them  depending  on  how  that  material  is  involved  in  the  fire 
growth  process.  Note  that  since  the  available  total  energy  per  unit  area  in  a  material  was  fixed,  there  was 
a  trade-off  between  the  chosen  magnitude  of  the  constant  HRR  and  the  burning  duration. 

Here,  for  some  materials,  calculated  two  types  of  square  wave  description  were  calculated:  one  which 
emphasized  the  initial  spike  in  HRR  and  a  second  which  emphasized  either  the  overall  test  average  HRR 
or  some  compromise  between  these  two  extremes.  In  principle,  a  fire  growth  model  could  switch  from 
the  first  to  the  second  of  these  two  descriptions  at  some  pre-determined  extent  of  burning  (fractional 
weight  loss  or  fractional  release  of  the  total  available  heat),  if  this  were  deemed  desirable  to  more  closely 
capture  the  true  time-dependence.  The  calculated  values  are  presented  below. 

Once  an  average  HRR  is  defined,  the  available  nine  values  from  the  cone  calorimeter  tests  (three  rephcate 
tests  at  three  heat  fluxes)  can  be  used  to  define  a  least  squares  straight  line  fit  to  predict  the  HRR  (as 
approximated  by  a  square  wave)  at  any  external  radiant  heat  flux  (but  see  discussion  below).  The  slope  of 
this  line  is  a  measure  of  the  effective  heat  of  gasification  of  the  material.  In  principle,  the  heat  of 
gasification  is  the  amount  of  heat  required  to  convert  unit  mass  of  a  material  at  room  temperature  into 
flammable  gases  during  steady  burning.  Thus,  from  a  steady-state  heat  balance  at  the  burning  surface  of 
the  material: 

g^,,,-mL  (A-1) 

where: 

rii  is  the  flux  of  gases  from  the  surface  of  a  burning  material  and  L  is  the  heat  of  gasification^ .  It  can  be 
shown  from  Equation  A-4  below  that 

L  =  AH^  I  {Slope  of  HRR  vs  q^^„.^,  )  (A-2) 


The  desk  surface  heat  release  rate  behavior  happens  to  closely  follow  the  predictions  of  a  published  model  for  wood  behavior 
(Ritchie  et  al.  1997),  in  spite  of  the  detaching  top  laminate  layer.  This  model  thus  may  be  used  to  give  a  more  complete  picture 
of  the  time-dependent  heat  release  rate. 
^  In  steady  burning  the  heat  of  gasification  is  the  sum  of  the  sensible  heat  required  to  heat  unit  mass  of  a  material  from  room 
temperature  to  its  burning  surface  temperature  plus  the  amount  of  heat  required  to  convert  that  mass  into  gases  at  the  burning 
surface  temperature. 
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where: 

A//^,  is  the  heat  of  combustion  of  the  sample,  and  the  slope  referred  to  is  that  of  a  plot  of  some 
appropriate  measure  of  heat  release  rate  versus  the  external  flux  absorbed  by  the  sample. 

Also,  by  definition, 

and  q  ^!^,„,^.  is  the  heat  flux  to  the  surface  from  the  flame,  ^„„,,,/  is  the  absorbed  amount  of  external 
radiation  (from  the  cone  calorimeter  heater)'',  and         is  the  radiation  lost  from  the  hot  surface  to  the 
surroundings.  Equation  (A-1)  says  that  the  net  heat  flux  into  the  condensed  phase  goes  into  generating 
the  gaseous  mass  flux  m  at  a  cost  of  L  kJ/g  of  gasified  material.  When  applied  here,  this  heat  balance 
glossed  over  non-steady  behavior  that  is  inherently  present  in  complex,  real  materials  due  to  such  things 
as  char  accumulation,  an  external  heat  flux  that  changes  on  a  time  scale  faster  than  the  condensed  phase 
themial  wave  equilibration  time,  a  layered  composition,  or  a  finite  sample  thickness.  Those  effects 
(except  for  a  time-varying  external  flux)  were  present  in  the  cone  calorimeter  tests  but  they  were  being 
averaged  into  the  apparent  value  of  L.  Thus,  as  used  here,  L  was  the  apparent  or  effective  amount  of  heat 
required  to  sustain  the  gasification  process  over  a  duration  of  burning  that  corresponds  to  the  defined 
average  HRR  (and  total  heat  available)^.  This  would  carry  over  to  a  full-scale  fire  with  reasonable 
accuracy  if  the  factors  affecting  the  averaged  transient  effects  were  essentially  the  same  on  both  scales. 
Here  this  meant  that  the  physical  behavior  of  the  sample  had  to  be  the  same  on  both  scales,  and  the  rear 
surface  boundary  conditions  also  had  to  be  nominally  the  same.  There  was  a  reasonable  expectation  that 
the  physical  behavior  would  be  the  same  on  both  scales  -  up  to  a  point.  That  point  corresponded  to 
physical  collapse  of  the  various  workstation  structures,  the  chair,  the  computer  housing,  etc.  These 
various  collapses  would  occur  at  differing  times  in  the  course  of  a  full-scale  workstation  fire.  The  rear 
surface  boundary  conditions  were  essentially  adiabatic  in  the  cone  calorimeter  tests.  They  were  not  so  far 
from  adiabatic  in  the  full-scale  as  to  lead  one  to  expect  appreciable  deviations  here.  An  exception  might 
have  occurred  when  the  underside  of  the  desk  surface  began  to  bum  and  did  so  long  enough  for  its 
associated  thermal  wave  to  interact  with  the  thermal  wave  associated  with  the  top  surface  (roughly  1000  s 
after  local  ignition  of  the  second  side). 

The  FDS  model  actually  required  the  rate  of  mass  loss  from  the  local  surface  of  a  burning  object  resulting 
from  a  given  total  incident  heat  flux,  plus  the  heat  of  combustion  of  the  evolved  gases.  The  model 
calculated  this  total  incident  flux,  used  it  to  obtain  the  mass  evolution  rate  from  the  surface,  then 
distributed  the  combustion  heat  from  that  mass  in  the  gas  phase  in  accord  with  the  buoyant,  turbulent  flow 
field  the  model  calculates.^  In  calculating  the  total  incident  heat  flux,  FDS  did  not  distinguish  between 
"external  radiation"  and  "flame  heat  flux"  in  the  manner  shown  in  Equation  (A-3).  Thus,  the  FDS  model 
could  not  directly  use  plots  of  average  heat  release  rate  versus  external  radiant  flux  since  the  abscissa 
there  was  the  external  radiant  heat  flux  from  the  cone  heater,  not  total  incident  heat  flux  to  the  surface. 


*  The  external  radiant  flux  is  assumed  to  pass  undiminished  through  the  flame  since  this  layer  is  thin. 
^  The  burning  duration  is  (Total  Available  Heat/ Average  HRR). 

^  The  model  assumes  infinitely  fast  gas  phase  oxidation  kinetics  so  that  the  combustion  takes  place  on  the  surface  where  fuel 
vapors  and  oxygen  meet  in  stoichiometric  proportions.  That  surface  location  and  shape  is  constantly  changing  due  to  the 
turbulence. 
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The  cone  samples  respond  to  the  combined  external  radiant  flux  and  the  flame  feedback  flux  so  their  total 
incident  heat  flux  is  greater  than  the  indicated  abscissa  values  on  such  plots'. 

Heat  release  rate  is  simply  the  product  of  mass  loss  rate  and  heat  of  combustion  of  the  evolved  gases. 
Using  this  fact  plus  Equations  (A-1)  and  (A-3),  it  can  be  shown  that,  in  principle,  the  zero  flux  intercept 
on  a  plot  of  average  HRR  versus  external  flux  is  a  measure  of  the  flame  feedback  flux  which  the  samples 
in  the  cone  calorimeter  see.  Thus: 

Here,  the  first  term  on  the  right  is  the  intercept  of  such  a  plot;  if  the  other  quantities  in  this  term  were 
known  from  other  measurements  or  estimates^  then  g       can  be  calculated.  In  practice,  given  the  scatter 
in  the  data  points  that  define  a  fitted  straight  line,  plus  the  approximations  that  this  model  encompasses, 
use  of  this  intercept  is  a  very  poor  way  in  which  to  infer  the  flame  feedback  flux.  It  is  much  better  to  rely 
on  measured  values  of  heat  fluxes  from  flames  to  surfaces  in  bench  scale  experiments  such  as  this 
(Quintiere  and  Rhodes  1994;  Ohlemiller  1999;  Kashiwagi  2003).  These  values  tend  to  be  in  the  range 
from  25  kW/m"  to  40  kW/m".  In  the  present  work  a  value  of  33  kW/m'  was  assumed  to  be  an  adequate 
estimate  for  all  of  the  materials.  This  imposed  an  added  point  through  which  a  straight  line  fit  to  the  cone 
data  had  to  pass  (since  that  line  was  governed  by  Equation  (A^).  The  slope  of  the  least  squares  line  was 
thus  altered  to  force  it  through  this  intercept  value.  Since  the  heat  of  gasification,  L  ,  appears  in  both 
ternis  of  equation  3  above,  NIST  proceeded  iteratively  to  get  a  value  for  it  when  fitting  a  line  to  the 
combination  of  the  nine  measured  HRR  data  points  plus  the  forced  point  at  an  external  flux  of  zero. 

The  plots  in  Figs.  A-7  through  A- 12  are  a  result  of  the  above  process.  The  lowest  point  on  the  left  in 
each  plot  corresponds  to  the  imposed  flame  flux  of  33  kW/m'.  The  fitted  line  was  forced  through  this 
point  but  gave  a  least  squares  fit  to  the  remaining  points  which  were  averaged  HRR  values  from  the  cone 
calorimeter  experiments.  The  abscissa  is  not  the  external  flux  measured  in  the  cone  calorimeter  tests  but 
rather  an  estimate  of  the  total  flux  seen  by  the  sample  surface.  Thus  the  abscissa  is  simply  the  sum  of  the 
cone  calorimeter  flux  plus  the  estimated  33  kW/m'  flame  flux.  This  is  more  in  keeping  with  the  heat  flux 
the  FDS  model  would  use,  though,  in  fact,  the  model  used  the  gasification  heat  (not  the  fitted  line  in  the 
Figures)  to  get  mass  loss  rate  from  total  incident  heat  flux  via  Equation  (A-1),  as  described  above,  letting 
the  model  distribute  the  resulting  heat  release  in  the  gas  phase.  In  effect,  however,  the  model  followed 
the  HRR  values  predicted  by  the  fitted  lines  in  the  Figures.  Thus,  it  predicted  heat  release  rates  in  accord 
with  the  cone  calorimeter  experiments. 


^  The  other  significant  estimate  here  is  that  the  surface  re-radiation  loss  is  the  same  for  all  of  the  materials  at  a  value  of 
lOkW/m^.  This  correspond  to  a  surface  temperature  of  375  °C.  '~ 
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Figure  A-7a.  Heat  flux  dependence  of  averaged  HRR;  workstation 
desk  surface;  HRR  averaged  within  80  percent  of  peak  HRR. 
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Figure  A-7b.  Heat  flux  dependence  of  averaged  HRR;  workstation 
desk  surface;  HRR  averaged  over  entire  burn  duration. 
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Figure  A-8a.  Heat  flux  dependence  of  averaged  HRR;  chair  seat 
assembly;  HRR  averaged  within  80  percent  of  peak  HRR. 
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Figure  A-8b.  Heat  flux  dependence  of  averaged  HRR;  chair  seat 
assembly;  HRR  averaged  over  entire  burn  duration. 
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Figure  A-9.  Heat  flux  dependence  of  averaged  HRR;  HRR  averaged 
within  80  percent  of  peak  HRR 
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Figure  A-10.  Heat  flux  dependence  of  averaged  HRR;  workstation 
privacy  panel;  HRR  averaged  within  80  percent  or  peak  HRR. 
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Figure  A-11a.  Heat  flux  dependence  of  averaged  HRR;  paper  stack; 
HRR  averaged  within  80  percent  of  peak  HRR. 
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Figure  A-1 1  b.  Heat  flux  dependence  of  averaged  HRR;  paper  stack; 
HRR  averaged  within  60  percent  of  peak  HRR. 
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Figure  A-12a.  Heat  flux  dependence  of  averaged  HRR;  computer  monitor  shell;  HRR 

averaged  within  80  percent  of  peak  HRR. 
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Figure  A-12b.  Heat  flux  dependence  of  averaged  HRR;  computer  monitor  shell;  HRR 

averaged  within  50  percent  of  peak  HRR. 
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The  values  of  the  heats  of  gasification  in  Table  A-2  were  sometimes  calculated  in  two  ways  to  emphasize 
either  early  or  overall  behavior  (as  described  above)  for  the  six  material  combinations  characterized  here. 
The  first  two  parameters  pertain  to  ignitability  and  are  discussed  below.  The  heat  of  combustion  was,  as 
noted  above,  obtained  from  the  cone  calorimeter  tests  as  the  ratio  of  the  HRR  to  the  mass  loss  rate  from 
the  sample.  These  are  both  instantaneous  values,  and  the  heat  of  combustion  may  vary  with  time.  Only 
the  sample  mass  was  actually  measured  as  a  function  of  time.  The  mass  loss  rate  was  obtained  by 
differentiating  mass  versus  time,  which  could  yield  quite  a  noisy  result.  Thus,  the  heat  of  combustion 
data  from  the  cone  calorimeter  were  filtered  and  then  averaged,  typically  over  a  range  from  10  percent  to 
90  percent  weight  loss  from  the  sample.  Each  heat  of  gasification  was  obtained  from  the  slope  of  a  fitted 
line  in  Figs.  A-7  through  A- 12  and  the  average  heat  of  combustion.  The  total  available  heat  per  unit 
sample  area,  which  can  vary  with  heat  flux,  is  the  time  integral  of  the  HRR  curve  and  is  reported  out  by 
the  cone  calorimeter  software.  The  burning  times  reported  in  the  Table  are  derived  from  the  total 
available  heat  and  the  average  HRR,  as  described  above.  They  were  not  used  directly  by  the  FDS  model 
and  are  shown  just  to  indicate  their  magnitude  (and  their  variation  among  the  six  materials).  Note  that  the 
fabric  cover  on  the  privacy  panels  had  a  bum  time  of  less  than  1  min  while  the  desk  surface  had  a  bum 
time  that  approached  30  min.''*  For  the  upper  values  of  heat  of  gasification  and  of  buming  time  in  the 
Table,  the  HRR  curves  were  averaged  around  the  peak  values;  typically  the  averaged  values  were  those 
within  80  percent  of  the  absolute  peak  value,  a  technique  used  with  some  success  in  (Dillon,  Kim,  and 
Quintiere  1998).  This  averaging  helped  smooth  some  of  the  noise  inherent  in  a  peak,  particularly  from  a 
physically  complex  composite.  The  lower  values  were  typically  close  to  the  overall  test  average  HRR. 

As  discussed  in  Sec.  3.2.2,  some  of  the  materials,  notably  the  fabric  on  the  privacy  panels  and  the  overall 
chair,  behaved  in  a  substantially  different  way  than  in  the  above  cone  calorimeter  tests.  This  was  a  result 
of  their  thermoplastic  namre.  As  a  result,  it  was  ultimately  necessary  to  treat  the  buming  of  these  in  a 
more  heuristic  manner. 

A.4         IGNITION  BEHAVIOR 

Piloted  ignition  delay  time  was  measured  at  several  incident  radiant  heat  fluxes,  as  described  above.  The 
two  ignition  parameters  reported  in  Table  A-2  were  derived  by  matching  the  experimental  ignition  time 
data  with  a  one-dimensional  numerical  model  of  heat  conduction  in  a  homogeneous  solid  of  finite 
thickness  with  constant  thermal  properties.  The  front  boundary  condition  includes  an  incident  radiant 
flux  (here  taken  to  be  absorbed  at  the  surface),  a  natural  convective  loss,  and  a  gray-body  re-radiative 
loss.  The  rear  boundary  condition  was  taken  to  be  adiabatic  since  the  cone  calorimeter  samples  were 
insulated  on  their  back  side.  The  front  surface  convective  loss  was  estimated  for  a  horizontal  cone 
calorimeter  sample  using  a  standard  heat  transfer  correlation. 

The  ignition  temperature  was  first  estimated  by  repeated  runs  of  the  model  near  the  minimum  flux  for 
ignition.  The  thermal  inertia  was  then  estimated  by  matching  the  measured  ignition  delay  time  at 
20  kW/m"^  (somewhat  higher  for  the  desk  surface).  This  low  flux  matching  was  chosen  to  enhance  the 
likelihood  that  the  FDS  model  would  correctly  predict  the  early  fire  growth  on  the  workstation.  The 
thermal  inertia  thus  obtained  was  an  effective  value  that  absorbs  any  deviations  between  the  inert  heat 
conduction  model  and  the  actual  behavior  of  the  sample.  Thus,  for  example,  the  endothermic  heat  of 
melting  of  the  nylon  fibers  on  the  face  of  the  carpet  was  implicit  in  the  inferred  value  of  thermal  inertia 


These  are  square  wave  bum  times.  The  real  materials  tend  to  tail  off  gradually  and  thus  to  bum  weakly  for  a  longer  time. 
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for  the  carpet.  Typically,  thermal  inertia  values  inferred  in  this  manner  were  appreciably  larger  than  the 
product  of  handbook  values  for  thermal  conductivity,  density  and  heat  capacity.  Nevertheless,  ihcy  were 
quite  capable  of  reproducing  the  observed  ignition  behavior.  Figures  A- 13  through  A- 18  show  both  the 
experimental  ignition  behavior  measured  in  the  cone  calorimeter  and  the  calculated  behavior  based  on  ilie 
ignition  temperature  and  effective  thermal  inertia  given  in  Table  A-2. 

The  ignition  behavior  shown  in  these  figures  is  all  based  on  the  use  of  a  pilot  (in  this  case  an  electric 
spark)  in  the  gas  phase  just  a  few  centimeters  above  the  gasifying  surface.  Thus,  as  soon  as  an  ignitable 
mixture  was  present  in  this  space  immediately  above  the  sample  surface,  it  did  ignite.  The  full-scale 
experiments  on  actual  workstations  had  no  such  guarantee  of  immediate  ignition.  As  discussed  in  the 
above  section  on  the  workstation  fire  tests,  this  led  to  behavior  that  was  dependent  on  the  appearance  of 
pilot  flames,  which  were  often  provided  by  paper  piles  that  could  auto-ignite. 
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Figure  A-13.  Workstation  desk  surface;  ignition  delay  time  vs.  incident  heat  flux. 
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Figure  A-14.  Chair  seat  assembly;  ignition  delay  time  vs.  incident  heat  flux. 
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Figure  A-15.  Carpet;  ignition  delay  time  vs.  incident  heat  flux. 
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Figure  A-16.  Paper  stack;  ignition  delay  time  vs.  incident  flux. 
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Figure  A-17.  Computer  monitor  shell;  ignition  delay  time  vs.  incident  heat  flux. 
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Figure  A-18.  Workstation  privacy  panel;  ignition  delay  time  vs.  incident  heat  flux. 
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